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PREFACE 


This  document  serves  dual  purposes:  (l)  provides  a  status  report  on  the  Air 
Force's  basic  research  program  on  rocket  propulsion  and  (2)  is  the  program  for 
the  contractors'  meeting  on  rocket  research. 

Most  of  the  abstracts  follow  a  specific  format.  The  text  begins  with  a  short 
statement  of  relevant  scientific  questions  addressed  by  the  research,  followed 
by  an  explanation  of  the  scientific  approach.  A  statement  of  the  uniqueness  of 
each  approach  was  solicited  from  the  investigators.  The  major  portion  of  the 
text  describes  the  results  obtained  during  the  last  twelve  months.  The 
abstracts  should  describe  two  figures:  Figure  1  illustrates  the  main  (or  a 
representative)  feature  of  the  scientific  approach  and  Figure  2  presents  a 
primary  accomplishment. 

Hard  copies  of  the  vugraph  material  and  collateral  information  are  in  file 
folders  (one  for  each  presentation)  on  a  table  at  the  rear  of  the  meeting  room. 

Since  the  research  on  combustion  instability  and  advanced  ingredients  have  had 
adequate  national  forums  in  recent  months,  that  research  will  be  presented  in 
the  March  1985  Rocket  Research  Meeting.  The  abstracts  of  that  research  are 
included  in  this  booklet.  The  research  on  advanced  diagnostics  of  reacting  flow 
is  the  subject  of  a  separate  meeting  and  abstract  booklet. 

A  primary  objective  of  the  meeting  and  this  booklet  is  to  encourage  the 
participants  to  interpret  the  technological  barriers  and  to  consider  new 
research  approaches.  Since  a  25  to  30  percent  annual  turn-over  is  built  into 
the  program,  each  year  opportunities  exist  for  new  research  approaches  and  for 
new  principal  investigators.  Several  of  the  presentations  provide  introductions 
to  some  of  the  barriers  and  technological  challenges.  Prospective  principal 
investigators  should  not  feel  constrained  by  these  presentations  and  are 
encouraged  to  look  beyond  the  identified  items.  The  location  of  the  meeting 
(i.e.,  near  organizations  interested  in  the  research)  promotes  interchanges 
among  the  investigators  and  those  responsible  for  Air  Force  programs. 
Accordingly,  many  of  the  participants  will  be  able  to  provide  specific 
information  on  Air  Force  requirements.  Questions  concerning  the  meeting  can  be 
directed  to  either: 

Leonard  H  Caveny  or  Robert  A  Biggers 

AFOSR/NA  AFRPL/XRX 

Bolling  AFB  Edwards  AFB,  CA  93334 

Washington,  DC  20332  Phn:  (805)  277-5206 

Phn:  (202)  767-4937  Autovon;  550-5206 

Autovon:  297-4937 
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0800-1200 

0830 


AGENDA  SUMMARY 

1984  AFOSR/AFRPl  ROCKET  RESEARCH  MEETING 
March  12  -  Monday 


Registration 

Announcements 

Materials  and  Structural  Mechanics 
Combustion 

Combustion  (continued) 

Exhaust  Plumes 

Discussion  Symposia  (Snacks  and  Cash  Bar)  -  Essex  House 


March  13  -  Tuesday 

Registration 
Metallized  Propellants 
Announcements 

WELCOME  by  Don  A  Hart,  Director.  AFRPL 
Overview  -  Air  Force  Rocket  Propulsion, 

Richard  R  Weiss.  Chief  Scientist,  AFRPL 
AFOSR  Interests  In  Rocket  Propulsion,  Leonard  H  Caveny,  AFOSR 
AFRPL  Research  Interests  In  Propulsion,  Robert  A  Diggers,  AFRPL 
Shuttle  Based  Experiments,  William  Oran  and  Thomas  Labus,  NASA 
Thermodynamics  and  Kinetics 
Administrative  Meeting  for  AFOSR  Contractors 
Discussion  Symposia  (Snacks  and  Cash  Bar)  -  Essex  House 

March  14  -  Wednesday 
Late  Registration 
Beamed  Energy  In  Flowing  Media 
Electromagnetic  Acceleration  of  Plasmas 
Discussion  Symposia  (Snacks  and  Cash  Bar)  -  Essex  House 


I  .  .  J 


March  15  -  Thursday 
Energy  Conversion 

NASA  Sponsored  Research  on  Nonconventlonal  Propulsion, 
Thomas  H  Cochran  and  David  C  Byers,  NASA 
Working  Sessions 
Adjourn  Meeting 
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1984  AFOSR/AF RPL  ROCKET  RESEARCH  MEETING 
Essex  House,  Lancaster,  CA 
MONDAY  (AM) 

1?  March  1984 

0730  Registration 

0820  Announcements 

Session  Chairman:  Ourwood  I  Thrasher,  AFRPI/MKP6 
Topic:  Materials  and  Structural  Mechanics 
TIME  NUM. 

0830  1  CARBON-CARBON  PROCESSING  VARIABLES.  W  P  Hoffman,  AFRPL/MKBN 

0900  2  ST A I ISTICAL  NATURE  OF  CRACK  GROWTH  IN  SOLIO  PROPELLANTS.  Chi  T 

(Jimmie)  Liu,  AFRPL/MKPB 

0930  3  PROPELLANT  NONLINEAR  CONSTITUTIVE  THEORY  EXTENSTION.  Gene 

Francis,  Chemical  Systems  Division,  United  Technologies, 
Sunnyvale,  CA 

100U  4  PROPFIIANT  AGING  RESEARCH.  OougTas  B  Olsen,  Aerochem  Research 

Laboratories,  Inc,  Princeton,  NJ  (New  Start) 

TOTS  Break 


Session  Chairman:  Jay  N  Levine,  AFRPL/OYC 
Topic:  Combustion 

1030  5  MECHANISMS  FOR  ACOUSTIC  SUPPRESSION.  Merrill  W  Beckstead,  R  1 

Raun,  and  P  C  Bralthwalte,  Brigham  Young  University,  Provo,  UT 

1100  b  COMBUSTION  MECHANISMS.  David  P  Weaver  and  J  T  Edwards,  AFRPl 2DYP 
David  H  Campbell  and  Susan  C  Hulslzer,  University  of  Dayton 
Research  Institute 

1130  /  IHFRMOPHYSICAL  PR0PER1Y  DETERMINATIONS  USING  1 RANSIENT 

TECHNIQUES.  Raymond  E  Taylor,  Purdue  University, 

West  Lafayette,  IN 


1200 


lunch  (Reconvene  at  1300) 
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HONDA Y  (PH) 
12  HARCH  1984 

Session  Chairman:  Jay  N  Levine,  AFRPL/OYC 
Topic:  Combust on  (continued) 

Tine  NUH. 


% 

1300 

8 

CHEHICAL  KINETICS  OF  NITRAHINE  PROPELLANT  COHBUSTION.  Helvyn  C 
Branch,  University  of  Colorado,  Boulder,  CO  (New  Start) 

1330 

9 

RAPID  SCAN  FTIR  SPECTROSCOPY  OF  NIIRAHINES.  Thomas  B.  Brill, 
Richard  J  Karpowlcz,  and  Yoshlo  Oyuml,  University  of  Delaware, 
Newark,  DE 

1400 

10 

SOL 10  PROPELLANT  TRANSIENT  RESPONSES.  Robert  L  Gllck,  Purdue 
University,  West  Lafayette,  IN 

11 

*FUEL  RICH  PROPELLANT  COHBUSTION  BEHIND  A  BLUFF  BODY.  Warren  C 
Strahle  and  Jechlel  (Jeff)  I  Jagoda,  Georgia  Institute  of 
Technology,  Atlanta,  GA 

1430 

Break 

Session  Chairman:  David  P  Weaver,  AFRPL/DYP 


Topic: 

Exhaust  Plumes 

1500 

12 

SUPRESSION  OF  AFTERBURNING  IN  SOLID  ROCKET  PLUHES  BY  POTASSIUH 
SALTS.  Eugene  Hiller,  University  of  Nevada,  Reno,  NV  (New  Start) 

1530 

13 

AFTERBURNING  SUPPRESSION  KINETICS.  Jay  D  Eversole, 
Dayton  Research  Institute 

University  of 

1600 

14 

COHBUSTION  KINETICS  OF  HETAL  OXIDE  AND  HALIOE  RADICALS.  Arthur 
Fontljn,  Rensselaer  Polytechnic  Institute,  Troy,  NY 

1630 

ADJOURN  SESSION 

1630 

Discussion  Symposia  (Snacks,  soft  drinks,  cash  bar) 

-  Essex  House 

‘Abstract  only,  no  presentation. 


TUESDAY  (AH) 

13  March  1984 

0730  Registration 

Session  Chairman:  Wayne  E  Roe,  AFRPL/DYC 
Topic:  Metallzed  Propellants 
TIME  HUH. 

0600  15  WIDE  DISTRIBUTION  PROPELLANTS.  Robert  A  Frederick,  John  C  Matson, 

and  John  R  Osborn,  Purdue  University,  West  Lafayette,  IN 

0030  16  BEHAVIOR  OF  ALUMINUM  IN  SOL  10  PROPELLANT  COMBUSTION.  Edward  W 

Price  and  Robert  K  Slgman,  Georgia  Institute  of  Technology, 
Atlanta,  GA 

0900  1/  AERODYNAMIC  BREAKUP  -  METAL  DROPLETS.  James  E  Craig,  Spectron 

Development  Laboratory,  Costa  Mesa,  CA 

0930  IH  AGGLOMERATE  BREAKUP  CONTROL.  Kevin  K  Nack,  AE RPL/OYC 

0945  Break 

1030  19  PUI  SEIl-LASER  CINEMATOGRAPHY  OF  PROPELLANT  DEFLAGRATION.  Roger 

Becker,  Unlv  of  Dayton  Research  Institute. 

1100  20  HOLOGRAPHIC  DATA  RETREV1AL.  Oavld  Netzer,  U  S  Naval  Postgraduate 

School,  Monterey,  CA 

21  *FUEI  -RICH  SOLID  PROPELLANT  BORON  COMBUSTION  (ANALYSIS). 

Merrill  King,  Atlantic  Research  Corp,  Alexandria,  VA 

22  *FUEI  -RICH  SOLID  PROPELLANT  BORON  COMBUSTION  (EXPERIMENT). 

James  J  Komar  and  Ronald  S  Fry,  Atlantic  Research  Corp, 

Alexandria,  VA 


1130  Lunch  (Reconvene  at  1300) 


‘Abstract  only,  no  presentation. 


1984  AFOSR/Af RPL  CHEMICAL  ROCKET  RESEARCH  MEETING 
Essex  House,  Lancaster,  CA 
TUESDAY  (PH) 

13  March  1984 


Session 

Chairman:  Leonard  H  Caveny,  AFOSR/NA 

TINE  NUN. 
1300 

Announcements 

1305 

WELCOME : 

Don  A  Hart,  Director,  AFRPL/CC 

1315 

OVERVIEW: 

Scientist 

AIR  FORCE  ROCKET  PROPULSION.  Richard 
.  AFRPL 

R.  Weiss, 

Chief 

1415 

23 

OVERVIEW: 

AFOSR/NA. 

AFOSR  INTEREST  IN  ROCKET  PROPULSION. 

Leonard  H 

Caveny 

1445 

OVERVIEW: 
A  Blggers 

AFRPL  RESEARCH  INTERESTS  IN  ROCKET  PROPULSION. 

,  AFRPL/XRX 

Rober  t 

1500 

Break 

1530 

Invited: 

SHUTTLE  BASED  COMBUSTION  EXPERIMENTS. 

William  Oran, 

Topic: 

1545 


1615 


NASA  Washington,  DC  and  Thomas  tabus,  NASA  Lewis  Research  Center, 
Cleveland,  OH 


Thermodynamics  and  Kinetics 


24  CRITICAL  EVALUATION  AND  COMPILATION  OF  THE  THERMODYNAMIC 

PROPERTIES  OF  HIGH  TEMPERATURE  SPECIES.  Malcom  W  Chase,  The  Dow 
Chemical  Company,  Midland,  MI 


25  CRITICAL  EVALUATION  OF  HIGH  TEMPERATURE  CHEMICAL  KINETIC  DATA. 
Norman  Cohen  and  Karl  Westberg,  Aerospace  Corporation,  Los 
Angeles,  CA 


1645 

1645-1715 


ADJOURN  SESSION 

ADMINISIRATIVE  MEETING  FOR  AFOSR  CONTRACTORS  ONLY 


1645 


Discussion  Symposia  (Snacks,  soft  drinks,  cash  bar)  -  Essex  House 


WEDNESDAY  (Am 
14  March  1984 

NONCONVENTIONAL  PROPULSION 

0730  Registration 

Session  Chairman:  Clark  Hawk,  AFRPL/LK 
Topic:  Beamed  Energy  to  Flowing  Media 

TINE  NUH. 

0800  26  LASER  THERMAL  PROPULSION.  Dennis  Keefer,  University  of  Tennessee 

Space  Institute.  Tullahoma,  TN 

0830  27  LASER- INOUCED  GAS  BREAKDOWN  IN  THE  VISIBLE/UV.  David  I  Rosen  and 

Guy  H  Weyl,  Physical  Sciences  Inc,  Andover,  HA 

0900  28  HIGH  TEMPERATURE  MOLECULAR  ABSORBERS  FOR  CW  LASER  PROPULSION. 

David  0  Rosen,  David  0  Ham,  and  Lauren  N  Cowles,  Physical  Sciences 
Inc,  Andover,  HA 


Break 


1000  29  PROPAGATION  OF  LASER  SUPPORTED  PLASMAS  IN  FLOWING  MEDIA.  Charles 

L  Merkle  and  Thomas  M  York,  The  Pennsylvania  State  University, 
University  Park,  PA 

30  HIGH  POWER  Nd-GLASS  LASER  INSTRUMENT  FOR  ADVANCED  PROPULSION  AND 
DIAGNOSTICS.  Charles  L  Merkle  and  Thomas  M  York,  The  Pennsylvania 
State  University,  University  Park,  PA 

1030  31  PLASMA  INITIATION  MECHANISMS,  ABSORPTION  OR  CW  LASER  ENERGY  AND 

MIXING  FLUIO  MECHANICS.  Herman  Krler,  Jyotl  Mazumder,  and  Ronald 
J  Glumb,  University  of  Illinois,  Urbana-Champalgn,  IL 

1100  32  Invited:  EXPERIMENTAL  STUOY  OF  LASER  SUPPORTED  HYDROGEN  PLASMAS. 

T  Dwayne  McCay,  0  M  VanZandt,  and  R  H  Eskridge,  NASA/Marshall 


1130  33  AFRPL  SOLAR-THERMAL  ROCKET  ACTIVITIES.  Curtis  C  Selph,  AFRPL/LKCS 

1200  LUNCH  (Reconvene  at  1300) 


Xi  - 


WEDNESDAY  (PH) 
UHarch  1984 

Session  Chairman:  Robert  J  Vondra,  AFRPL/LKCJ 


loplc:  Electromagnetic  Acceleration  of  Plasmas 


TINE  NUM 

1300  34  BASIC  PROCESSES  IN  PLASMA  PROPULSION.  Herbert  0  Schrade  and 

Michael  Klrschner,  Instltut  fur  Raumfahrtantrlebe  der  Unlversltat 
Stuttgart 


1330  35  MPO  THRUSTER  PERFORMANCE.  R  Joseph  Cassady,  Purdue  University, 

West  Lafayette,  IN 

Robert  J  Vondra  and  W  M  Schmidt,  AFRPL,  Edwards  AFB,  CA 


1400  36  MAGNETOPLASMAOYNAMIC  THRUSTER  EROSION  RESEARCH.  David  Q  King,  Jet 

Propulsion  Laboratory,  Pasadena,  CA 

1430  37  EROSION  AND  0NSE1  MECHANISMS  IN  MAGNETOPLASMADYNAMIC  THRUSTERS. 

John  L  Lawless,  Carnegle-Mellon  University,  Pittsburgh,  PA 

1500  Break 


1530  38  PERF0RMANCE-LIMITA1 ING  FACTORS  FOR  HPD  THRUSTERS.  Manuel 

Martlnez-Sanchez,  Massachusetts  Institute  of  Technology, 
Cambridge,  MA 

1600  39  PULSED  INDUCTION  ENERGY  TRANSFER.  Peter  Mongeau  and  Henry  Kolm, 

Electromagnetic  Launch  Research,  Inc,  Cambridge,  MA, 

1630  40  UNIFIEO  STUDY  OF  PLASMA-SURFACE  INTERACTIONS  FOR  SPACE  POWER  AND 

PROPULSION.  Peter  Turchl,  R  and  D  Associates  Washington 
Laboratory,  Alexandria,  VA 

1700  ADVANCED  DIAGNOSTICS  FOR  PLASMA  FLOWS.  Ronald  K  Hanson,  Stanford 

University,  Stanford,  CA 

42  *N0NC0NVENTI0NAL  ELECTROMAGNETIC  PROPULSION  CONCEPTS.  Michael  M 
Mlccl,  Pennsylvania  State  University,  University  Park,  PA 

43  *CLOSED-SPACE  TEMPERATURE  KNUDSEN  FLOW.  John  B  McVey,  Rasor 
Associates,  Inc,  Sunnyvale,  CA 

1730  Discussion  Symposia  (Snacks,  soft  drinks,  cash  bar)  -  Essex  House 


‘Abstract  only,  no  presentation. 
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1030  41  AN1 1  PROTON  ANNIHILATION  PROPULSION.  Robert  L  Forward,  Foward 
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1300  WORKING  SESSIONS: 

DAI  A  SKIS  FOR  MPD  PHYSICAL  MODEL  VALIDATION.  (Discussion  leader  to 
be  determined) 

INSIRUMENTATION  TO  ESTABLISH  PROPERTIES  OF  LASER  SUPPORTED 
PLASMAS.  (Discussion  leader  to  be  determined) 

1500  Adjourn  Meeting 


'Abstract  only,  no  presentation. 
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CARBON-CARBON  PROCESSING  VARIABLES  INVESTIGATION 

W.  P.  HOFFMAN 

AIR  FORCE  ROCKET  PROPULSION  LABORATORY 
AFRPL/MKBN  EDWARDS  AFB,  CA  93523 


The  chemistry  and  physics  of  carbon-carbon  composites  is  not  well  understood  even 
though  they  are  being  evaluated  as  nozzle  materials  in  a  number  of  Air  Force 
programs.  This  lack  of  understanding  includes  the  surface  ablation  kinetics  of  these 
materials.  There  is  conjecture  about  the  specie  (s)  principally  responsible  for 
recession.  However,  there  is  no  definitive  data  available  either  pertaining  to  the 
principle  specie  (s),  or  to  the  kinetics  of  the  heterogeneous  reactions  with  carbon 
under  rocket  nozzle  conditions.  Therefore,  the  nozzle  designer  has  no  firm  theoretical 
basis  on  which  to  predict  the  performance  of  his  nozzle  design.  He  is  forced  to  make 
predictions  based  on  data  from  actual  motor  firings.  The  nozzle  and  conditions  used  in 
these  firings  may  or  may  not  represent  the  composite  material  and  operating 
conditions  that  the  designer  has  chosen.  The  objective  of  this  program  is  to  fill  this 
gap  in  our  knowledge  of  surface  reaction  kinetics  and  their  relation  to  carbon-carbon 
material  construction  and  processing.  This  is  being  acomplished  by  (I)  extensive 
surface  characterization  of  the  carbon-carbon  materials  and  (2)  by  measuring  the 
kinetics  of  the  heterogeneous  reactions  between  carbon-carbon  materials  and  the 
various  species  present  in  the  rocket  motor  exhaust.  This  will  give  the  designer  a  firm 
basis  on  which  to  predict  the  performance  of  his  nozzle. 

During  the  past  year,  characterization  studies  of  two  commercially  prepared  carbon 
fibers  (T-300  and  VSB-32)  that  had  been  heat  treated  to  2700oC  was  continued.  Work 
has  centered  on  gaining  better  understanding  of  fiber  structure.  In  addition,  much 
effort  has  been  expended  in  refining  the  techniques  developed  to  measure  the  total 
surface  area  (TSA),  active  surface  area  (ASA),  and  porosity  of  these  fibers  as  well  as 
composites  manufactured  from  them.  Very  precise  techniques  are  required  since  the 
parameters  measured  are  typically  very  small  for  adsorption  techniures.  For  instance, 
the  TSA  of  these  fibers  is  usually  0.3m2/g  and  the  ASA  is  typically  five  percent  of 
this  voiue.  One  of  the  major  modifications  completed  this  year  was  in  the  automated 
data  acquisition  system.  This  modified  the  signal  to  noise  (S/N)  ratio  by  an  order  of 
magnitude  which  revealed  information  previously  obscurred.  This  improvement  in  the 
S/N  ratio  made  in  possible  not  only  to  confirm  the  microporosity  of  these  fibers  but 
also  to  quantify  it.  Although  these  fibers  have  been  found  to  have  similiar  TSA's,  they 
have  significantly  different  ASA's  and  pore  size  distributions.  The  intriging  fact 
however  is  that  although  there  are  differences  between  the  fibers,  these  differences 
are  not  as  significant  as  one  would  expect  considering  the  diversity  of  the  precursor. 
For  the  two  fibers  studied,  the  ratio  of  ASA's  in  an  oxygen  atmosphere  was  similar  to 
the  ratio  of  recession  rates  of  billets  manufactured  from  these  fibers. 
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STATISTICAL  NATURE  OF  CRACK  GROWTH  IN  SOLID  PROPELLANTS 

C.  T.  LIU 

Air  Force  Rocket  Propulsion  Laboratory 
Edwards  AFB,  CA  93523  (805)  277-5502 

The  high  volume  concentration  of  fillers  of  different  sizes  in  a  composite  solid 
propel lapt  produces  a  heterogeneous  material.  The  heterogeneity  of  the  micro¬ 
structure  will  undoubtedly  affect  the  material  properties,  as  well  as  the  failure 
behavior  of  solid  propellants.  Thus,  depending  upon  the  local  microst^ucture  or  the 
local  stress  and  strength  near  the  crack  tip,  the  rate  of  crack  growth,  -m  may  show  a 
considerable  fluctuation  within  a  specimen,  or  it  may  show  considerably  larger  scatter 
among  specimens  even  though  the  testing  conditions  remain  identical.  In  addition  to 
the  intrinsic  scatter  of  due  to  the  heterogeneity  of  the  microstructure,  the 
methods  of  the  -r-  calculation  may  also  be  a  contributing  factor  to  the  scatter  of  the 
rate  of  crack  growth  due  to  the  discrete  nature  of  the  crack  growth  data.  Therefore, 
it  is  the  purpose  of  this  study  to  investigate  the  effect  of  the  methods  of  the  -re¬ 
calculation  on  the  variability  of  the  crack  growth  rate  and  the  statistical  character¬ 
istics  of  the  crack  growth  behavior  in  a  solid  propellant. 

The  crack  growth  behavior  of  a  highly  filled  composite  propellant  under  constant 
strain  rate  condition  was  studied  through  the  use  of  centrally  cracked  biaxial 
specimens.  In  this  study,  a  total  of  38  replicate  tests  were  conducted  under  identical 
conditions.  The  raw  data  (load,  time,  and  crack  length)  obtained  from  the  tests  were 
used  to  calculate  the  instantaneous  crack  growth  rate  and  the  stress  intensity  factor 
for  a  given  instantaneous  crack  length  and  the  associated  load.  In  calculating  the 
crack  growth  rate,  four  different  methods  were  used.  The  calculated  crack  growth 
rate  and  the  stress  intensity  factor  were  analyzed  statistically  to  determine  the 
functional  relationship  between  the  crack  growth.rate  and  the  stress  intensity  factor, 
as  well  as  the  statistical  distribution  function  of  a  function  of  the  stress  intensity 
factor.  A  detailed  discussion  on  the  methods  of  4^  calculation  and  the  statistical 
analysis  was  given  at  the  1983  Air  Force  Office  of  Scientific  Research/Air  Force 
Rocket  Propulsion  Laboratory  (AFOSR/AFRPL)  Chemical  Rocket  Research  Meeting 
and  subsequently  shown  in  the  meeting's  report.  Therefore,  those  subjects  will  not  be 
discussed  in  this  abstract.  However,  for  an  illustrative  purpose,  the  schematic 
representation  of  the  basic  approach  is  reproduced  and  shown  in  Figure  I. 

Based  on  the  analysis,  the  following  conclusions  are  drawn: 

I.  The  method  of  calculation  has  a  significant  effect  on  the  variability  of  the 


data 


da 


2.  The  secant  method  introduces  the  least  amount  of  error  into  the  g^data. 


da 


data. 


3.  Data  smoothing  action  introduces  a  large  error,  but  less  scatter,  into  the 


4.  For  the  range  of  K  |  consid^yd  in  this  study,  the  2-parameter  log-normal 
distribution  provides  a  best  fit  for  the  gy  data. 

5.  A  power  law  relationship  exists  between  the  crack  growth  rate,  4?,  and  the 
stress  intensity  factors,  and  mathematically,  it  can  be  written  as 

3T=  C2KI  C| 

6.  The  parameters,  C,  and  C2,  are  mutually  dependent,  and  all  the  crack  growth 
lines—plotted  as  In  a  vs  In  K|  --  teria  to  converge  to  a  small  region. 


Figure  1.  Schematic  Representation  of  Basic  Approach 
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(a)  Coefficient  of  Variation  of  77-  vs  Stress  Intensity  Factor 
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PROPELLANT  AGING  RESEARCH 
Douglas  B.  Olson 

AeroChem  Research  Laboratories,  Inc. 

P.O.  Box  12,  Princeton,  NJ  08542 

The  overall  objective  of  this  work  is  to  gain  a  fundamental  understanding 
of  the  chemical  reaction  mechanisms  and  kinetics  of  solid  rocket  propellant 
aging.  If  the  aging  chemistry  of  the  individual  propellant  ingredients  can  be 
accurately  modeled,  then  the  useful  lifetimes  of  current  and  improved  formula¬ 
tions  will  be  predictable.  The  approach  taken  (see  Fig.  1)  is  to  separately 
study  the  two  main  steps  in  the  propellant  aging  process:  the  decomposition 
of  oxidizers  and  plasticizers;  and  the  attack  of  the  highly  reactive  decompo¬ 
sition  products  on  the  binders,  stabilizers,  and  other  propellant  ingredients. 
Mechanisms  will  be  postulated  (using  the  measured  kinetic  parameters)  and  the 
predictions  of  this  chemical  aging  "model"  will  be  compared  with  available 
data  on  mechanical  property  changes  that  occur  from  long  term  aging  of  propel¬ 
lants.  Experiments  to  investigate  the  effects  of  minor  ingredients,  e.g., 
bonding  agents  or  combustion  modifiers,  will  be  performed  as  needed  to  improve 
the  predictive  accuracy  of  the  model. 

Measuring  oxidizer  and  plasticizer  decomposition  rate  coefficients  at 
near  ambient  temperature  as  in  this  work  requires  ultrasensitive  techniques. 
Gas  evolution  from  the  primary  decomposition  steps  of  several  oxidizers  and 
plasticizers  (RDX,  AP,  and  BTTN)  is  being  investigated  using  the  AeroChem  NOx/ 
NHj  +  03  chemiluminescence  analyzer  which  can  measure  real  time  decomposition 
rate  coefficients  as  small  as  about  10_la  s_I  (half-life  greater  than  10'* 
years),  corresponding  to,  e.g.,  an  NO  concentration  of  about  1  part  in  1010. 
Obtaining  data  such  as  these  at  realistic  propellant  storage  temperatures  is 
important  (see  Fig.  2)  because  it  has  been  found  that  rate  coefficients  ob¬ 
tained  at  higher  temperatures  (above  as  100°C)  sometimes  have  significantly 
different  activation  energies  from  those  at  lower  temperature,  and  give  gross¬ 
ly  incorrect  estimates  of  compound  stabilities  when  extrapolated  to  ambient 
temperature . 

The  reactivity  of  two  propellant  binders,  HTPB  and  PEG,  with  and  without 
stabilizer,  toward  several  gaseous  decomposition  products  from  the  oxidizers 
and  plasticizers  will  be  measured.  Various  techniques  will  be  used  in  these 
studies,  including  measurements  of  gas  evolution  rates  from  binder-oxidizer  or 
plasticizer  mixtures  (using  the  N0x/NH3  analyzer),  binder  oxidative  chemilumin¬ 
escence,  and  physical  property  changes  (using  several  methods  including  gel 
permeation  chromatography) . 

Experimental  measurements  on  RDX  decomposition  are  currently  underway 
using  the  N0X  analyzer.  Previous  unpublished  results  (labeled  AeroChem  in 
Fig.  2)  gave  an  activation  energy  of  about  88  kJ  mol-1  over  the  temperature 
range  60  to  130°C,  significantly  lower  than  the  higher  temperature  values. 

The  measurements  in  this  program  are  being  made  from  as  60°C  down  to  20°C  using 
an  improved  apparatus. 


•V*  -V.V.V.-.-A  . 
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PROPELLANT  AGING  RESEARCH 

PROBLEM 


APPROACH 

MEASURE  OXIDIZER/PLASTICIZER 
DECOMPOSITION  RATES 

RDX,  AP,  BTTN,  T  =  20-60°C 
Use:  •  NOx/NH3  +  O3  Analyzer 

MEASURE  BINDER  REACTIVITY  WITH 
DECOMPOSITION  GASES  (NO,  N02,  HCI04,  ...) 
Use:  e  NO*/NH3  +  O3  Analyzer 

e  Binder  Chemiluminescence 
e  Physical  Property  Changes 

DEVELOP  CHEMICAL  AGING  MODEL 
HTPB  and  PEG 

PREDICT  PROPELLANT  SERVICE  LIFETIMES? 


DECOMPOSITION  RATE  COEFFICIENT 
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T,  °C 


300  200  100  50  20 


103/T,  K~' 


RDX  and  BTTN  Decomposition  Data 

ADDITIONAL  DATA  WILL  BE  OBTAINED  FROM 
T  =  20  TO  60°C  FOR  USE  IN  THE 
PROPELLANT  AGING  MODEL. 
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MECHANISMS  FUR  ACOUSTIC  SUPPRESSION 

Principal  Invistigator:  Merrill  W.  Beckstead 
Coinvestigators:  R.  L.  Raun  and  P.  C.  Braithwaite 
Brigham  Young  University,  Provo,  Utah  86402 

Acoustic  suppressants  are  commonly  added  to  low  smoke  and  smokeless 
propellants  to  avoid  the  problem  of  combustion  instability.  However,  the 
basis  for  selecting  the  type,  size,  and  concentration  of  suppressant 
particle  for  optimum  performance  is  not  well  founded.  Suppressants 
apparently  work  by  one  or  more  of  three  mechanisms: 

(1)  energy  loss  due  to  viscous  dissipation  due  to  drag  forces, 

(2)  modification  of  the  propellant  combustion  response  function,  or 

(3)  energy  interchange  due  to  distributed  combustion. 

Viscous  dissipation  of  acoustic  energy  by  drag  forces  on  a  suspended 
particle  is  the  most  commonly  recognized  mechanism,  being  a  direct  analogy 
to  suspended  moisture  in  a  fog.  A  second  mechanism  whereby  an  additive  can 
effect  acoustics  occurs  at  the  burning  boundary  of  the  system.  Solid 
additives  within  the  propellant  can  nave  a  catalytic  influence  on  the 
propellant  combustion,  modifying  the  combustion  response.  The  best  known 
example  of  this  is  probably  AI0O3.  Al^  is  the  particulate  product  of 
aluminum  combustion  which  provides  the  viscous  drag  dissipation  previously 
discussed.  However,  A1203  is  also  a  weak  combustion  catalyst  in  composite 
propellants.  An  additive  that  influences  the  steady  burning  of  a  propellant 
can  also  be  expected  to  influence  the  transient  response  of  the 
propellant.  The  third  mechanism  considered  is  due  to  the  effect  of  a 
particle  burning  as  it  traverses  a  relatively  large  portion  of  the  system. 
As  it  does  so,  the  interchange  of  energy  between  the  burning  particle  and 
the  acoustic  environment  can  result  in  either  a  driving  or  damping 
contribution  to  the  acoustics  of  the  system.  This  third  mechanism  is  the 
primary  area  of  study  of  the  current  contract. 

Figure  1  is  a  graphic  illustration  of  the  approach,  some  of  the 
technical  advantages,  the  data  to  be  collected  and  the  data  analysis  to  be 
performed.  The  experimental  basis  for  the  technical  approach  is  the  Rijke 
tube,  a  gas  burner  with  a  paddle  to  control  (stop  and  start)  combustion 
oscillations,  over  frequencies  ranging  from  500  to  1500  Hz.  Advantages  of 
this  approach  include  allowing  a  separation  of  the  three  mechanisms 
mentioned  above  by  testing  an  additive  independent  of  the  propellant  burning 
surface,  avoiding  use  of  solid  propellants,  and  allowing  for  independent 
control  of  frequency,  0/F  ratio,  temperature  and  particle  addition. 

The  principle  objective  of  this  work  is  to  identify  and  develop  an 
understanding  of  the  mechanisms  whereby  acoustic  suppressants  modify  an 
acoustic  wave.  It  is  expected  that  several  mechanisms  will  be  shown  to  act 
together  to  dampen  acoustic  oscillations,  with  one  mechanism  perhaps 
dominating  under  certain  conditions  and  another  mechanism  dominating  at 
other  conditions. 

Figure  2  summarizes  the  activities  of  the  past  year. 
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ACCOMPLISHMENTS 


1.  LITERATURE  SURVEY  COMPLETED 

•  Rijke  Tube  Acoustics 

•  Particle  Damping  in  Solid  Propellants 


2.  MECHANISMS  FOR  PARTICULATE  ACOUSTIC  SUPPRESSION  IDENTIFIED 

•  Energy  loss  due  to  viscous  dissipation  of  drag  forces 

•  Modification  of  the  propellant  combustion  response  function 

•  Energy  interchange  due  to  distributed  combustion 


3.  BASIC  EQUATIONS  OUTLINED 

•  Conservation  Equations  for  Gas  Plus  Particles  Defined 

•  Neglect:  wall  friction,  viscous  dissipation,  kinetic  energy 
terms,  body  forces 

•  Assume:  one-dimensional  flow,  constant  transport  properties, 
constant  area,  small  amplitude  oscillations,  discontinuity  at 
flame 

•  Perturbation  form  of  equations  obtained 


4.  EXPERIMENTAL  PROGRESS  WITH  RIJKE  BURNER 

•  A  quantitative  particle  feeder  developed 

•  Improved  paddle  designed  and  constructed 

•  Test  procedures  established 
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How  Do  Suppressants  Work  in  Solid  Propellant  Rocket  Motors? 


Mechanisms  of  Acoustic  Suppressants: 


Viscous  dissipation  due  to  drag  forces 
Modification  of  propellant  response  function 
Effect  of  distributed  combustion 


CURRENT  PROGRAM: 


•  Separate  Effects  of  Different  Mechanisms  Using 
Unique  Experiments 

•  The  Modified  Rijke  burner 

•  Use  gas  (eliminates  solid  propellant) 
t  Evaluate  acoustics  w/o  particles 

•  Can  feed  particles  independently 

•  Vary  particle  type,  size  and  concentration 

•  Evaluate  reactive  or  inert  particles 
t  Use  paddle  to  control  acoustic  growth 

Particles 


Paddle 


•  Results 


•  Monitor  particle  combustion  with  high  speed  movies 

•  Evaluate  acoustic  growth  rates  with  and  without  particl 

•  Evaluate  mechanisms  of  reactive  versus  inert  particles 
t  Measure  effect  of  pertinent  parameter 
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COMBUSTION  MECHANISMS 

O.P.  Weaver,  D.H.  Campbell,  S.C.  Hulsizer , J.T.  Edwards 
Air  Force  Rocket  Propulsion  Laboratory 
Edwards  AFB,  California 


The  objective  of  all  propellant  development  programs  is  the 
ability  to  tailor  propellant  ballistic  properties  and  performance 
through  the  selective  modification  of  propellant  ingredients. 
While  the  influence  of  additives  on  burning  rate  and  combustion 
instabilities  has  been  well  documented,  as  yet  no  detailed 
understanding  of  the  specific  chemical  mechanisms  or  their 
influence  on  propellant  combustion  has  been  developed. 
Information  on  the  influence  of  such  detailed  chemical  kinetics 
changes  on  the  overall  propellant  burning  behavior  is  of 
fundamental  importance  in  establishing  the  influence  of 
composition  variables  on  macroscopic  properties  such  as 
performance  and  ballistics.  The  objective  of  this  research  is  to 
provide  detailed  chemical  kinetic  mechanisms  of  the  gas  phase 
reacting  species  at  or  near  the  surface  of  HTPB/AP  propellant 
samples  burning  at  pressures  typical  of  solid  rocket  motor 
chambers.  Additionally,  chemical  kinetic  and  statistical 
modelling  of  the  burning  process  is  to  be  coupled  to  these 
results  to  aid  in  the  interpretation  of  the  specific  chemistry 
of  the  reaction  process. 

On  a  microscopic  scale,  the  combustion  of  a  solid  rocket 
motor  propellant  is  seen  to  be  a  heterogeneous,  unstable  process 
exhibiting  multiple  flames,  local  wide  variation  of  gas  phase  0/F 
ratios,  and  ignition  delays  of  certain  particles.  However,  the 
aggregate  propellant  sample  burns  at  a  measurable,  steady  rate  on 
the  macroscopic  scale.  This  forms  the  basis  for  statistical 
approaches  to  combustion  modelling  and  requires  sufficient 
spatial  and  temporal  resolution  of  diagnostic  capability  as  to 
isolate  the  combustion  mechanisms  on  the  microscopic  scale  of  the 
intermittent  burning.  Further,  the  information  must  be  specific 
to  the  molecular  and  atomic  species  within  such  localized 
regions.  The  accumulation  of  such  localized  specific  data  will 
identify  the  competing  mechanisms  and  lead  to  an  understanding  of 
the  effect  on  burning  rate  of  chemistry  variation.  As  a  first 
step,  mechanisms  leading  to  the  localized  0/F  ratio  variations 
must  be  specified  and  the  affect  on  these  ratios  of  specific 
propellant  chemistry  detailed.  The  experimental  apparatus  used  in 
this  study  along  with  a  represen  tat i on  of  the  resulting  data  is 
presented  in  Figure  1.  Emphasis  at  this  point  is  on  coherent  anti 
Stokes  Raman  scattering,  laser-induced  fluorescence,  and  emission 
spectroscopic  surveys  of  HTPB/AP  propellants  in  the  pressure 
range  200-1000  PSI .  Typical  results  on  OH  are  shown  in  Figure  2. 
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FIGURE  2:  OH  EMISSION  FROM  I1TPB/AP  PRO  TP  1.1  ANT,  750  PSI 
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THEKiiOPHYSICAL  PROPERTY  DETERlIIh'ATIOI'S  UiXL.C 


vi;Au:;iLi.T  tlchuiuuel 


Principle  Investigator:  a.  i..  Tayior 
Thcri.opnysical  Propertied  Lesearea  Laboratory 
Purdue  University 
Lest  Lafayette,  Indiana 


Knowledge  of  the  theruophysicsl  properites  of  rochet  propellants,  especially  at 
higher  temperatures  where  significant  decomposition  occurs,  is  necessary  lor  under¬ 
standing  various  physical  processes  which  take  place  in  the  combustion  of  these 
propellants.  ilueh  effort  has  been  put  into  the  development  of  computer  codes  to 
model  various  aspects  of  rocket  propellant  behavior,  such  as  combustion  mechanics  and 
DSDT.  However,  there  is  a  lack  of  accurate  thermophysical  property  data  for  the  pro¬ 
pellants  being  used  in  these  models  in  addition  to  the  lack  of  understanding  of  the 
Lasic  heat  transfer  mechanisms.  In  this  program,  ueterminations  of  the  thormopim  si- 
cal  properties  of  AP,  RDX,  and  HhX  are  being  completed  in  an  attempt  to  quantitively 
understand  heat  transfer  through  the  propellants  under  transient  nesting  at  tempera¬ 
tures  near  combustion. 

Shown  in  Figure  1  is  the  general  approach  for  making  the  theruophysioai  property 
determinations  of  a  material  when  the  thermal  conductivity  is  difficult  to  measure 
directly  because  of  sample  instability.  Traditional  thermal  conductivity  measure¬ 
ments  require  large  samples  and  steady  state  thermal  gradients  for  accurate  results 
while  the  pulse  difrusivity  method  uses  small  samples  and  transient  temperature  rises 
of  the  order  of  a  few  degrees  C.  The  pulse  difrusivity  method  therefore  allows  the 
investigation  of  transient  heat  transfer  in  a  sample  while  the  sample  is  undergoing 
changes  which  more  closely  resembles  the  conditions  of  a  burning  propellant.  Two 
thermal  properties  are  being  determined:  (1)  specific  heat  using  a  Perkin-Elmer 
iioacl  DSC-2  Differential  Scanning  Calorimeter,  and  (2)  thermal  diffusivity  using  a 
laser  flash  apparatus.  All  measurements  are  made  under  digital  computer  control  for 
data  acquisition  and  analysis.  The  product  of  these  results  and  the  densities  of  the 
propellants  yields  the  thermal  conductivity  of  the  material. 

Specific  heat  results  have  been  obtained  for  all  three  materials.  The  effects 
of  partial  decomposition  on  the  specific  heat  have  also  been  studied.  These  results 
will  be  presented. 

The  rate  of  diffusion  of  heat  through  a  material  depends  upon  the  mean  free  path 
length  which  in  turn  is  a  function  of  the  crystal  structure  (including  imperfections) 
and  intergranular  contacts.  The  mean  free  path  lengths  within  the  crystal  is  signi¬ 
ficantly  influenced  by  imperfection  sites  due  to  phase  changes,  acuitives, 
dcconposition  and  lattice  distortion.  The  diffusion  path  length  for  pressed  powuers 
is  dependent  upon  several  variables:  the  grain  size,  the  intergrain  contact  resis¬ 
tance.  the  structural  integrity,  the  grain  composition  and  the  packing  density.  As  a 
result  the  thermal  diffusivity  of  a  powdered  propellant  is  highly  dependent  upon  the 
propellant's  composition,  particle  size  ana  its  structural  condition.  If,  for  exam¬ 
ple,  a  propellant  bed  has  been  exposed  to  acoustic  vibrations,  uicro-crackin^,  coulu 
occur  which  would  lower  the  thermal  diffusivity  of  the  propellant.  Therefore,  know¬ 
ing  a  propellant's  composition  is  not  sufficient  to  assign  an  accura  to  thermal 
diffusivity  value  to  it. 

Thu  approach  taken  has  been  to  determine  the  thermal  diffusivity  of  pure  crys¬ 
tals  and  powders  of  the  propellants  of  interest.  Single  crystal  and  pov;dcred  AP 
samples  have  been  measured  as  well  as  powdered  samples  of  HEX.  The  thermal  cifx'usiv- 
ity  along  two  principle  crystal  axes  for  AP  were  investigated.  The  uncorrecttc 
thermal  diffusivity  values  along  the  principle  axes  increase  as  tne  molecular'  Done 
lengths  decrease  as  expected.  These  results  are  plotted  in  Figure  2.  The  AP  crystal 
is  transparent  to  infared  light  and  thus  it  was  necessary  to  form  a  layered  composite 
by  placing  the  sample  between  two  opaque  materials.  The  data  presented  therefore 
require  the  application  of  a  new  heat  loss  correction  technique  that  is  presently 
under  development. 
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gure  2.  The  Uncorrected  Thermal  Diffusivity 
of  Two  Orientations  of  Single  Crystal  AP. 
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CHEMICAL  KINETICS  OF  NITRAMINE  PROPELLANT  OOMBUSTION 
Melvyn  C.  Branch 

Mechanical  Engineering  Department 
University  of  Colorado 
Boulder,  Colorado  80309 

Although  nitramine  solid  rocket  propellants  have  advantages  of  high 
specific  impulse  and  low  infrared  and  smoke  emissions  from  their  exhaust, 
better  control  of  the  burn  rate  is  needed.  Propellant  additives  have  been 
suggested  for  burn  rate  control  but  the  understanding  of  the  chemical  mechan¬ 
ism  of  burning  is  not  well  enough  developed  to  predict  the  effects  of  these 
burn  rate  modifiers.  The  objective  of  the  studies  reported  here  is  to  des¬ 
cribe  the  detailed  chemistry  of  reactions  of  gas  phase  species  which  are 
important  in  the  combustion  of  nitramine  propellants  and  additives.  Species 
composition,  temperature  and  burning  rate  are  being  measured  for  reacting  gas 
mixtures  representative  of  the  pyrolysis  products  of  nitramines  and  nitra- 
mines  with  additives.  The  reaction  is  followed  experimentally  by  precise, 
spatially  resolved  measurements  of  species  concentration  and  temperature 
profiles  above  a  one  dimensional,  laminar  flat  flame  burner.  Chemical  kin¬ 
etic  and  fluid  mechanic  modeling  of  the  flame  structure  is  also  being  used 
to  correlate  the  experimental  data  and  evaluate  critical  reaction  paths  and 
energetics  and  their  importance  in  burn  rate  control. 

Preliminary  chemical  kinetic  modeling  of  the  nitramine  system  are  one 
dimensional,  adiabatic  reaction  without  diffusion  of  heat  and  mass.  The 
reaction  mechanism  is  based  on  existing  studies  of  the  NH3-NO-O2  reaction 
system  (Miller,  Branch  and  Kee,  1981).  These  mechanisms  include  reactions 
of  the  major  nitrogen  species  NO,  NO2 »  N2O,  NH3  and  N2  and  reactions  in  the 
H2/O2  system.  The  reaction  mechanism  also  includes  reactions  of  CH2O,  00  and 
CO2  and  their  intermediates  taken  from  the  recent  review  of  Westbrook  and 
Dryer  (1983) .  The  mechanism  can  only  be  regarded  as  tentative  and  prelimin¬ 
ary,  but,  nonetheless,  several  interesting  results  are  derived  from  its  use. 

The  reaction  has  a  distinct  two  zone  structure.  The  rapid  first  stage 
is  reaction  between  CH2O ,  NO2  and  N2O  and  the  intermediates  H  and  OH.  These 
reactions  are  responsible  for  most  of  the  heat  release  of  the  reaction. 

These  first  stage  reactions  also  keep  the  radical  concentrations  (H,0H,0)  low 
so  that  second  stage  reactions  are  slow.  The  primary  reactions  in  the  second 
stage  are  those  reducing  NO.  These  reactions  are  also  exothermic  and 
increase  the  gas  temperature.  This  stage  of  the  reaction  does  not  begin, 
however,  until  the  CH2O  and  NO2  are  depleted  and  the  radicals  (primarily  H 
and  OH)  are  able  to  increase  in  concentration.  Although  reactions  between 
NH2  and  NO  also  deplete  the  NO  concentration  these  reactions  are  not  dominant 
because  NH2  is  present  only  in  low  concentration.  Attempts  to  provide 
sources  of  NH3  or  NH2  to  the  propellant  to  accelerate  this  reaction  process 
are  under  study  (Flanagan,  et  al,  1982).  Other  calculations  at  lower 
pressures  show  that  the  sensitivity  of  the  reaction  zone  thickness  and 
reaction  rate  to  pressure  is  similar  to  that  of  the  burning  propellant. 

These?  preliminary  calculations  suggest  the  feasibility  of  describing  the  two 
zone  combustion  and  the  effects  of  pressure  on  this  structure  by  a  compre¬ 
hensive  chemical  kinetic  mechanism. 
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SOLID  PHASE 


Figure  1.  Schematic  representation  of  the  reaction  zone  near  the  surface  (solid 
line)  of  a  burning  nitramine  propellant.  The  "fizz"  reaction  zone  in 
the  gas  phase,  the  gas  preparation  zone  and  gas  flame  zone  are 
distinguished.  This  study  is  concerned  with  describing  the  chemistry 
of  these  regions  dominated  by  gas  reactions. 


MOLE  FRACTION 
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Figure  2.  Products  and  product  temperature  from  the  one  dimensional 
adiabatic  reaction  of  the  products  of  decomposition  of  RDX 
at  40  atm.  Initial  mixture  is  30^0  +  NO2  +  N2O  +  A  N2 . 

A  rapid  first  stage  reaction  between  CH2O,  NO2  and  *■ 

N2O  is  shown  which  ends  at  about  0.01  sec.  The  second  stage 
reaction  occurs  with  the  reduction  of  NO  and  begins  at  about 
0.01  sec. 


TEMPERATURE(K) 
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RAPID  SCAN  FTIR  SPECTROSCOPY  OF  NITRAMINES 

Thomas  B.  Brill,  Richard  J.  Karpowicz,  Yoshio  Oyumi 
Department  of  Chemistry,  University  of  Delaware 
Newark,  Delaware  19716 

The  identity  and  structure  of  the  chemical  species  in  the  solid  and  liquid 
phase  of  nitramines  near  their  melting  points  are  ill-defined.  An  experimental 
attack  with  the  goal  of  defining  the  nature  of  condensed  phase  of  energetic 
materials  in  real-time  has  been  undertaken  in  our  laboratory.  Our  current 
effort  is  concentrated  on  developing  and  using  rapid  scan  FTIR  spectroscopy  for 
this  role.  The  results  have  been  very  promising.  For  the  first  time,  an 
answer  has  been  provided  for  one  of  the  longer-standing  questions  of  importance 
in  the  nitramine  field;  i.e. ,  what  is  the  chemical  nature  of  RDX  and  HMX  melts? 
Also  for  the  first  time,  the  molecular  structure  of  RDX  in  the  gas  phase  and 
the  comparative  structures  of  RDX  and  HMX  in  the  solid,  solution,  melt  and  gas 
phases  have  been  addressed.  An  empirical  characterization  of  AZMTTC  was 
accomplished  by  joined  study  using  rapid-scan  FTIR  and  X-ray  crystallography. 

The  "melt"  phase  of  RDX  and  HMX  was  spectroscopically  characterized.  With  a 
heating  rate  of  150  C  sec-^  the  spectra  reveal  that  freshly  melted  RDX  (<3  sec 
old)  is  comprised  of  at  least  95%  intact  RDX.  The  steady-state  concentrations 
of  the  Intermediate  degradation  species  in  the  condensed  phase  are  below  the 
level  detectable  by  FTIR.  The  predominant  liquid  phase  degradation  product 
appears  to  be  an  amide  such  as  N-hydroxymethyl-f ormamide .  HMX  decomposition 
occurs  to  some  extent  before  the  melt  forms,  but  the  fresh  melt  is  mostly 
composed  of  intact  HMX. 

The  molecular  structure  of  RDX  in  the  a  and  B  solid  phases,  in  solution,  in  the 
melt,  and  in  the  gas  phase  was  determined  by  FTIR.  No  information  was 
previously  available  on  the  B-RDX  solid,  the  melt,  and  the  gas  phase.  The 
molecule  has  essentially  relaxed  to  C3V  symmetry  in  all  phases  except  the  a-RDX 
solid  phase  where  it  has  Cj  symmetry.  The  similarity  of  the  molecular  struc¬ 
ture  of  HMX  in  the  melt  to  that  in  CD3CN  solution  was  established.  It  has 
mostly  the  chair-chair  conformation. 

The  azldomethyl  derivative  of  HMX,  C4HgN4(N02>3(CH2N3) ,  called  AZMTTC  prepared 
by  Frankel,  Woolery  and  Flanagan  (Rocketdyne)  was  examined  by  X-ray  crystal¬ 
lography  and  rapid-scan  FTIR.  The  decomposition  mechanism  was  empirically 
determined.  An  unusually  long  C-N3  bond  exists.  The  initial  gas  phase  decom¬ 
position  product  is  HN3  from  cleavage  of  this  bond,  but  HN3  is  rapidly  eclipsed 
in  concentration  by  N2O  and  CH2O  formed  from  the  nitramine  portion  of  the 
molecules.  The  IR  spectrum  of  the  solid  phase  during  degradation  supports  this 
observation  as  does  the  sequence  of  weight  loss  from  a  TGA  study  of  AZMTTC. 

The  decomposition  is  triggered  by  liberation  of  a  small  amount  of  HN3  owing  to 
the  weak  C-N3  bond.  At  this  point,  the  ring  rapidly  depolymerizes  producing 
N2O  and  CH2O.  The  solid  residue  that  remains  then  finally  degrades. 
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HgCdTe 


Approach 

Rapid  Scan  FTIR  Studies  of  Nitramines 
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SOI  II)  PROPf  l  L  ANT  TRANSIf  Ml  RFSPONSES 

Robert  L.  Glick  and  John  R.  Osborn 
Purdue  University 
W.  Lafayette,  IN  47907 

Study  addresses  nonsteady  interactions  between  flowfield  in  semi- 
enclosed  duct  with  deflagrating  walls  of  heterogeneous,  condensed  media  and 
that  deflagration  process.  Figure  1  illustrates  basic  physics  and  approach. 
Monotonically  increasing  flow  Reynolds  number  leads  to  laminar,  transitional, 
and  fully  developed  turbulent  flow  sequence.  However,  deflagration  produced 
normal  flow  and  turbulence  kinetic  energy  lead  to  unique  features:  inflow 
stabilized  turbulence  fronts  and  sequential  turbul ization  and  profile  tran¬ 
sition.  As  deflagration  process  is  dependent  upon  pressure  and  fine  scale 
turbulence  close  to  wall  and  flowfield  depends  upon  both  deflagration  induced 
flow  and  turbulent  kinetic  energy,  flowfield/deflagration  coupling  is  also 
unique.  Key  scientific  questions  are  roles  of  turbulence  front,  bursting 
turbulence,  and  deflagration  produced  turbulence  in  nonsteady  interaction 
and  surface  population  fluctuations  in  heterogeneous  deflagration.  Approach 
is  to  define  nonsteady  burning  of  heterogeneous  media  and  turbulent  kinetic 
energy  production  experimental ly  with  local,  optical  remote  sensing  technique; 
flowfield  to  be  defined  with  combination  of  analysis,  gross  profile  measures, 
and  cold  flow  measures  in  another  program;  coupling  to  be  defined  experi¬ 
mentally  at  macro  (transfer  function)  level  in  semi -enclosed  duct  apparatus 
forced  into  axial  mode  oscillations.  All  studies  are  unique  to  this  program; 
local  burning  rate  technique  will  open  way  to  greatly  improved  understanding 
of  same  thru  two  point  correlations. 

Key  accomplishments  were  concerned  with  the  semi-enclosed  duct  apparatus 
and  flow  field  analysis.  Magnetic  velocimeter  design  was  thoroughly  reviewed; 
field  coil  design  was  improved;  circuit  constraints  needed  for  good  bandwidth 
were  developed.  Similarity  solution  obtained  for  nonsteady,  two-dimensional 
flowfield  near  velocity  antinode  in  porous  and  semi-porous  ducts  showed  sig¬ 
nificant  normal  velocity  effect;  acoustic  boundary  layer  on  blown  surface  is 
thickened.  For  blowing  Reynolds  numbers  >  ICr  (typical)  acoustic  layer  fills 
duct  while  acoustic  velocities  become  independent  of  viscosity  and  lead  mean 
acoustic  velocity  roughly  n/2  near  blown  wall.  Figure  2  illustrates  magnitude 
of  acoustic  velocity  for  semi-porous  duct  (blowing  Reynolds  number  20,000; 
frequency  Reynolds  number  2x10$;  permeable  wall  at  nondimensional  distance  0). 
Note  damped  oscillatory  character  of  the  acoustic  velocity  field  that  fills 
the  duct  and  overshoot  near  wall.  Acoustic  layer  on  impermeable  wall  is 
that  expected  on  blown  surface  showing  clearly  that  utilization  of  this 
thickness  as  small  parameter  for  flowfield/deflagration  studies  is  bad  physics. 
Out  of  phase  acoustic  velocity  near  transpired  surface  implies  reverse  flow 
(inflected  velocity  profiles)  at  amplitudes  far  below  those  for  mean  flow 
reversal.  Therefore,  potential  for  energy  transport  from  mean  to  nonsteady 
flow  seem  probable;  may  be  mechanism  for  strong,  non-linear  wall  heat 
flux/flowfield  coupling  observed  by  Brown.  Analysis  of  near  wall  region 
shows  these  effects  can  be  expected  in  tangential  mode  oscillations  and 
that  pressure  coupling  can  modify  phenomena. 


Re  =  2x10 
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FUEL  RICH  PROPELLANT  COMBUSTION  BEHIND  A  BLUFF  BODY 


Warren  C.  Strahle 
Jechiel  I.  Jagoda 

School  of  Aerospace  Engineering 
Georgia  Institute  of  Technology 

i 

The  solid  fueled  ramjet,  as  with  other  ramjet  types,  requires  a  flame 
anchoring  region  at  the  head  of  the  combustion  region.  This  is  frequently 
provided  by  a  recirculation  zone  behind  a  backward  facing  step  as  shown  in 
Fig.  1.  The  scientific  question  being  posed  is  whether  or  not  this  flow  field  is 
predictable.  The  flow  field  is  enormously  complex,  containing  interactions 
between  many  processes,  and  many  features  of  the  flow  absolutely  require 
advanced  diagnostics  for  interrogation  of  the  field.  On  the  other  hand,  there 
are  regions  of  the  flow  field  for  which  conventional  intrusive  diagnostics  are 
superior  to  the  non-intrusive,  optical  ones,  so  a  variety  of  techniques  are 
being  used. 

The  analytical  approach  is  currently  based  upon  k-e  modelling 
methods.  This  is  being  carefully  monitored  as  the  experimental  results  come 
in  from  more  and  more  complex  flows.  A  modification  has  already  been 
made  in  the  treatment  of  pressure  velocity  correlation  which  improves  the 
non-bleed  cold  flow  reattachment  length  in  comparison  between  theory  and 
experiment.  A  unique  feature  of  the  program  is  the  union  of  coworkers  in 
both  theory  and  experiment  at  the  same  laboratory  to  ensure  valid 
modelling. 

• 

Current  work  has  proceeded  through  the  cold-flow,  no-bleed  phase. 
An  example  of  the  many  techniques  used  for  measurement,  and  their 
comparison  with  theory  is  shown  in  Fig.  2,  where  the  mean  velocity  profile 
near  the  step  is  shown.  Shown  are  our  measurements  and  also  those  obtained 
at  another  laboratory  with  a  highly  similar  configuration.  Experiments 
compared  well  with  the  theory.  More  importantly,  in  regions  of  overlap  of 
validity  of  the  various  measurement  methods  (near  the  freestream),  there  is 
very  good  agreement  between  methods;  in  the  region  where  use  of  the  laser 
is  absolutely  required  (regions  of  reversed  flow,  intermittent  reversed  flow 
and  near  zero  velocity)  the  ability  to  obtain  data  is  substantiated. 


TYPICAL  TECHNICAL  ISSUES 
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LOW  FREQUENCY  RANDOM  OSCILLATIONS  REPORTED  IN  LITERATURE 
Careful  facility  design  eliminates  this  oscillation, 
suggesting  t.tany  results  are  configuration  sensitive. 


MEAN  STREAMWISE  VELOCITY 
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lift  SUP  I  'KE.  b'i  I UIJ  (Ji-  .’il  lERHilKNIlIk.  1 N  ',01.  I  p  I  .in  It  I 
FLUMES  By  POTASSIUM  SAL  I S 

F  I  lUENL  H  11  I  U< 

MACK  AY  !3CHt  (III  01  MINES,  I  HEM  f,  ME  I  l  NO  '  U  DU  I  . 

UNIVERSITY  Of  NEVADA,  RENO.  ITS1  BOSS 

In  recent  /ears,  the  Services  have  i ncreasi nq J v  emphasized 
the  development  and  use  at  law  signature  tactical  solid  rocket 
motors.  Visible  primary  and  secondary  smoke  have  been  1 aroel v 
eliminated  by  the  removal  ot  ammonium  perchlorate  oxidizer  arid 
most  of  the  other  energy  and  ballistic  modifier  additives  iron. 
the  propellant  formulation  -  the  so-called  "minimum  smoke 
(mi n-smoke)  propellants.  The  exhaust  qases  from  min -smoke 
propellants  however  contain  significant.  concentrati  ons  of 
hydrogen  and  carbon  monoxide  which  when  mixed  with  ambient  air 
m  the  plume  react  to  water  and  i_  rbon  dioxide  producing 
visible  flash  and  increased  infrared  radiation.  Also,  some  of 
the  apparent  secondary  smoke  advantage  of  mi n-smoke  propellant? 
over  reduced  smoke  propellants  (ammonium  perchlorate  oxidizer 
with  low  aluminum  content)  is  lost  since  the  hydrogen  in  the 
plume  reacts  to  form  additional  water  which  .is  available  for 
potential  condensation  to  smoke.  The  research  i s  directed 
toward  preventing  or  at  least  inhibiting  the?  signatures  due  to 
aft er burning. 

It  is  known  that  potassium  salts  inhibit  the  reactions  of 
hydrogen  and  carbon  monoxide  to  water  and  carbon  dioxide 
respectively.  Potassium  salts  such  as  KNOs  and  KaSO*  have  been 
added  to  propellant  charges  at  a  level  of  1  3  wt  pet  to 

suppress  gun  muzzle  flash  and  rocket  plume  infrared  signature.. 
The  mechanism  by  which  the  potassium  salts  inhibit  afterburning 
is  controversial ,  but  it  probably  involves  K,  KUH  and  possibly 
K0S.  reacting  with  H  and  OH  radicals  to  break  the  chain 
reactions  controlling  the  combustion  of  hydrogen  and  carbon 
monoxide.  Since  experimental  evidence  indicates  that  only  small 
concentrations  of  the  I  ,  KOH  and  K0»  are  required,  the  amounts 
of  salt,  that  have  been  used  previously  in  guns  and  rockets  mav 
h:'vr.'  been  excessive.  Minimizing  t.he  amount  used  i s  important 
because  the  potassium  salts  increase  radar  signatures  and  the 
propensity  for  smoke  formation. 

The  effects  of  potassium,  KOH  and  KOa  on  the  afterburning 
reactions  will  be  studied  by  introducing  them  as  a  vapor  into  a 
flai  diffusion  flame  of  methane-02-Na  and  scanning  the  flame 
incremental  1 y  by  infrared  and  visible  spectroscopy.  It  should 
he  possible  to  define?  the  conditions  of  gas  composition  and 
v  t.riip<->r  lure  for  which  they  are  effective.  An  opposed— jet 
dr  I  I  ijm  on  bur  ne-1  an»J  gas  flow  system  is  already  being  used  for 
related  afterburning  research.  A  modified  vaporizer-burner  head 
ls  bei  rig  built  to  permit  the  vaporization  of  the  potassium  and 
the  salts  prior  to  entering  the  flame.  A  Beckman  IF 
riper t:  ophotometer  has  been  modified  for  detection  of  emission 
spectra  and  the  scanning  device  is  in  place. 


MIN-SMOKE  SOLID  ROCKET  PLUME  AFTERBURNING : 
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EXPERIMENTAL  APPROACH 
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AFTERBURING  SUPPRESSION  KINETICS 
JAY  0.  EVERSOLE 

University  of  Dayton/  Air  Fore*  Rocket  Propulsion 
Laboratory,  Edwards  AFB,  CA  93523 

Methods  of  reducing  or  suppressing  rocket  afterburning  have 
been  pursued  largely  on  an  empirical  basis.  This  paper 
describes  current  research  efforts  to  combine  experimental 
studies  of  flame  suppression  with  a  computational  model  in  order 
to  determine  conclusively  the  detailed  chemical  kinetic 
mechanisms  for  specific  suppressant  species. 

The  experimental  arrangement  consists  of  a  premixed  flat 
flame  burner  mounted  on  a  motorized  X-Y-Z  platform.  The  flame  is 
currently  run  at  atmospheric  pressure,  but  will  eventually  be 
enclosed  in  a  vacuum  chamber  to  allow  operation  at  reduced  pres¬ 
sures.  Flow  rates  of  the  hydrogen,  oxygen,  and  nitrogen  to  the 
burner  and  the  temperature  of  the  porous  plug  are  monitored 
during  operation.  An  attempt  has  been  made  to  acheive  laminar 
one-dimensional  flow  in  the  flame  zone.  Gas  temperature  and 
molecular  hydrogen  concentration  profiles  as  a  function  of  axial 
distance  are  determined  by  spontaneous  Raman  scattering.  Hydro¬ 
gen  bromide  has  been  introduced  as  a  suppressant  species  into  the 
premixed  gases,  and  efforts  are  being  made  to  introduce  potassium 
vapor  as  a  suppressant.  Figure  1  is  a  schematic  represen tai on 
of  the  experimental  approach.  Additional  diagnostic  measurements 
will  be  made  to  determine  intermediate  specie  concentrations. 

The  computational  model  incorporates  detailed  reaction  kine¬ 
tics  with  laminar  one-dimensional  flow.  Resulting  temperature 
and  specie  mole  fractions  are  plotted  as  a  function  of  axial 
distance.  The  rate  constants  for  a  basic  set  of  seventeen 
elementary  reactions  of  the  hydrogen/oxygen  flame  are  known  with 
reasonable  accuaracy  from  previous  work,  and  comparisons  of 
calculated  flame  structure  to  published  data  and  calculations 
have  resulted  in  good  agreement  (Fig.  2).  Therefore,  the  effect 
of  adding  a  particular  suppressant  species  to  the  flame  may  be 
characterized  by  adding  a  hypothesized  set  of  reactions  (and 
their  rates)  involving  that  species  to  the  baseline  set.  For 
hydrogen  bromide  a  five-reaction  set  has  been  postulated  ,  and 
for  potassium  a  two-reaction  set  has  been  considered  .  Compari¬ 
son  of  predicted  quantities  and  experimental  measurements  are 
then  made  for  both  the  baseline  and  suppressed  cases.  For  con¬ 
clusive  identification  of  a  kinetic  mechanism,  agreement  between 
the  calculated  and  measured  profiles  must  be  acheived  over  a  wide 
range  of  initial  conditions  (pressure,  dilution,  suppressant 
concentration,  and  stoichiometry)  at  the  same  time  that  any 
alternate  reactions  (or  rates)  do  not  generate  the  same  results. 


1.  Dixon-Lewis,  G. ,  Proc  Royal  Soc  London,  292.  45-99  (1979). 

2.  Dixon-Lewis,  G.,  Combustion  and  Flame,  2£,  1-14  (1979). 

3.  Jensen,  D.E.,  Jones, G. A.,  Chem.  Soc.  Faraday  I,  73.,  2377-85 
(1979). 
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Figure  1.  Schematic  drawing  of  H^/O*  flat  flame  burner  and 
suppressant  species  injector,  and  laser  diagnostics  arrangement . 
Spatial  maps  (primarily  Z-axis)  of  temperature  and  species 
concentrations  are  obtained  from  data  (Raman  scattering,  LIF)  by 
translating  burner. 
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Figure  2.  Shows  calculated  Z-axis  profiles  of  temperature 
(upper  curves)  and  major  specie  concentrations  (lower  curves)  for 
baseline  and  suppressed  H^/O^/Nx  flame.  For  the  baseline  case 
the  narrow  flame  zone  is  located  at  about  0.1  cm  by  the  steep 
gradients  of  the  profiles.  When  7%  mole  fraction  of  potassium  is 
added,  the  flame  zone  is  extended  by  over  an  order  of  magnitude 
as  shown  by  the  much  smaller  gradient  of  all  the  profiles. 
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COMBUSTION  KINETICS  OF  METAL  OXIDE  AND  HALIDE  RADICALS 

Arthur  Fontijn 

Department  of  Chemical  Engineering  and  Environmental  Engineering 
Rensselaer  Polytechnic  Institute,  Troy,  NY  12181 

Current  ability  to  influence  combustion  involving  metals  is  severely  hampered  by 
a  lack  of  understanding  and  knowledge  of  the  ways  and  manner  by  which  tem¬ 
perature  affects  the  rate  coefficients  of  individual  reactions  and  reaction 
channels.  While  the  simple  Arrhenius-type  equation  k(T)  =  AT1/2  exp(-E^/RT)  has 
over  limited  temperature  ranges  been  of  great  value,  when  applied  to  wide  tem¬ 
perature  ranges  It  is  often  not  obeyed.  Particularly  for  exothermic  and 
slightly  endothermic  reactions,  order  of  magnitude  errors  can  be  made  by  extrapo¬ 
lations  based  on  the  Arrhenius  equation1.  It  is  the  goal  of  the  present  program 
to  provide  an  insight  in  the  kinetic  behavior  of  metallic  radical  oxidation 
reactions,  as  influenced  by  temperature,  based  on  accurate  measurements. 

Figure  1  illustrates  the  HTFFR,  high-temperature  fast-flow  reactor,  method  we  are 
using  to  generate  the  metallic  radicals  and  measure  their  oxidation  kinetics. 

This  unique  tool  provides  accurate  measurements  on  isolated  elementary  reactions 
in  a  heat  bath.  With  traditional  high  temperature  techniques  such  isolation  is 
usually  impossible  to  achieve;  as  a  result,  data  on  any  given  reaction  depend  on 
the  knowledge  of  other  reactions  occurring  simultaneously,  leading  to  large 
uncertainties.  Basically,  the  apparatus  consists  of  two  vertical  high  tem¬ 
perature  reactors  in  series.  In  the  upstream  reactor  the  metallic  radical  is 
produced,  BC1  in  the  example.  In  the  downstream  reactor  the  relative  con¬ 
centration  of  these  radicals  is  measured  by  laser-induced  fluorescence,  as  a 
function  of  [ 02] ,  P,  T  and  t,  to  determine  k(T).  The  sensitivity  of  the 
fluorescence  technique2  allows  observations  on  very  low  radical  concentrations 
(109-10n  cm-3),  which  prevents  interference  from  nucleation  reactions. 

Figure  2  shows  the  type  of  results  obtained  in  HTFFR  studies.  The  Al/SO?  reac¬ 
tion,  which  is  endothermic,  adheres  closely  to  the  Arrhenius  equation.  At0/0?, 
the  only  metallic  radical  reaction  for  which  information  is  available,  has  a  ^ 
slight  negative  temperature  dependence  due  to  intermediate  complex  formation. 

The  Ai/C02  reaction  increases  with  temperature  much  faster  than  Arrhenius  would 
predict,  probably  due  to  the  thermal  equilibrium  increase  in  the  concentration 
of  C02  in  bending  vibrational  modes  which,  contrary  to  ground-state  linear  C02, 
can  react  with  A i  at  larger  distances  by  an  electron  jump  mechanism.  This  beha¬ 
vior  of  the  A10/02  and  A1 /C02  reactions  was  not  predictable  and  indeed  unex¬ 
pected.  By  obtaining  data  on  a  number  of  metallic  radical  oxidation  reactions, 
we  anticipate  to  provide  a  basis  for  understanding  and  hence  prediction. 


1.  A.  Fontijn  and  R.  Zellner  in  "Reactions  of  Small  Transient  Species.  Kinetics 
and  Energetics",  A.  Fontijn  and  M.A.A.  Clyne,  Eds.,  Academic  Press,  1983, 
Chap.  1. 

2.  A.  Fontijn,  The  Use  of  Laser-Induced  Fluorescence  for  Fundamental  Gas-Phase 
Kinetic  Measurements",  Progress  in  Astronautics  and  Aeronautics,  in  press. 
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HIGH  TEMPERATURE  COMBUSTION 
KINETICS  OF  METALS 

SCIENTIFIC  APPROACH 

GOALS 

•  UNDERSTAND  THE  MECHANISMS  OF  METALLIC  RADICAL  OXIDATION 
REACTIONS 

•  MEASURE  REACTION  RATE  COEFFICIENTS  FOR  REACTIONS  IMPORTANT 
TO  ADVANCED  PROPULSION  SYSTEMS- 

•  BC<  +  O2  — ►  Emissions  from  Air  Breathing  Missile  Plume  After  Burning 

•  AIO  +C02  — *■  Combustion  of  Aluminized  Composite  Propellants 

•  BO  +  02  — »  Combustion  of  Boron  Slurries 

•  AlCl  +  02  — ►  Microparticulates  Formation 

•  AIF  +  02  — Tendency  of  Fluorinated  Binders  to  Reduce  Agglomeration 

•  OBTAIN  THESE  Jt(T)  DATA  OVER  THE  300 -2.000  K  RANGE 
NON -ARRHENIUS  BEHAVIOR  PRECLUDES  EXTRAPOLATION  FROM 
NARROW  TEMPERATURE  INTERVALS 

TECHNIQUE 

I.  PRODUCE  REACTANT  IN  HIGH  TEMPERATURE  FLOW  TUBE 


O 
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ACCOMPLISHMENTS 

OBSERVATION  INTERPRETATION 
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WIDE  DISTRIBUTION  PROPELLANTS 

Robert  A.  Frederick,  Jr.,  John  C.  Matson  and  John  R.  Osborn 
School  of  Aeronautics  and  Astronautics 
Purdue  University 
W.  Lafayette,  IN  47907 


Composite  propellants  containing  a  wide  AP  particle  size  distribution  have 
exhibited  unpredictable  ballistic  properties.  Highly  loaded  propellants  of  this 
type  do  not  burn  continuously  but  show  local  stepwise  burning  patterns  at  the 
propellant  surface.  The  propellant  burns  quickly  for  a  short  period  of  time, 
then  "rests"  for  a  short  period  of  time.  The  increments  of  burning  occur  over 
distances  comparable  to  the  diameter  of  the  coarse  AP,  while  the  rest  period 
has  been  linked  to  the  energetic  properties  of  the  binder  for  IPDI  and  DDI 
formulations. 

These  observations  challenge  previous  conclusions  that  propellant  burning 
rate  is  controlled  largely  by  the  deflagration  of  individual  oxidizer  crystals 
and  their  surrounding  binder.  Scenarios  which  describe  propellant  combustion  as 
strictly  a  time  average  or  area  average  process  are  inadequate  to  describe  this 
local  intermittent  combustion.  It  now  appears  possible  that  AP  particle  ignition 
and  binder  pyrolysis  are  controlling  combustion  mechanisms,  at  least  for  wide 
distribution  formulations. 

The  goal  of  this  research  program  is  to  verify  that  local  stepwise  burning 
of  solid  propellants  occurs  and  to  identify  the  physical  processes  that  control 
the  duration  of  the  burn  periods  and  "rest"  periods.  Perfecting  a  technique  to 
continuously  measure  the  local  burning  rate  is  a  major  task  in  achieving  this 
goal.  Measurements  will  first  be  completed  using  high-speed  photography  as  a 
preliminary  screening  tool.  The  desired  technique,  however,  is  a  direct  measure¬ 
ment  using  a  laser  position  detector  (LPD).  This  LPD  is  being  designed  and 
developed  to  continuously  record  spatially  local  (40  micron  diameter  circle), 
high  resolution  (20  micron)  distance  measurements.  With  an  adequate  technique 
perfected,  systematic  composition  variations  will  be  examined  to  determine  the 
roll  of  AP  ignition  and  combustion  properties,  and  the  energetic  nature  of  the 
binder  on  the  local  and  average  burning  rate.  Specifically,  the  binder  curative, 
oxidizer  modal  size  ratio,  and  pressure  will  be  investigated  as  independent 
variables. 

The  significant  achievement  to  date  has  been  the  continuous  measurement  of 
average  burning  rate  of  solid  propellant  burning  at  atmospheric  pressure  by 
sensing  the  position  of  a  modulated  laser  beam  reflected  from  the  propellant 
surface.  The  instrument  tested  compensated  for  the  variable  optical  properties 
of  the  propellant  surface  and  the  combustion  gases. 
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FEATURES  OF  APPROACH 


o  HETEROGENEOUS  SURFACE  CHEMISTRY 

-  LARGE  SIZE  DIFFERENCE  BETWEEN  COARSE  AND 
FINE  AP 

o  FIRST  DIRECT  MEASUREMENTS 

-  MEASUREMENT  OF  LOCAL  MICROSCOPIC  NON-STEADY 
BURNING 

o  SYSTEMATIC  VARIATIONS 

-  AP  DECOMPOSITION  CHARACTERISTICS 

-  BINDER  ENERGETICS  ' 


DIRECT,  CONTINUOUS 
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PROPELLANT  HEIGHT  (IN) 
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PRIMARY  SCIENTIFIC 
ACCOMPLISHMENTS... 

o  ACHIEVED 

-  CONTINUOUS  MEASUREMENT  OF  SURFACE  HEIGHT  AT 
ATMOSPHERIC  PRESSURE 


o  ANTICIPATED 

-  EFFECT  OF  BINDER  ENERGETICS  ON  MICROSCOPIC 
NON-STEADY  COMBUSTION 

-  EFFECT  OF  SURFACE  HETEROGENEITY  ON  AP 
IGNITION  DELAYS 
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BEHAVIOR  OF  ALUMINUM 
IN  SOLID  PROPELLANT  COMBUSTION 

Edward  W.  Price  and  Robert  K.  Sigman 
Georgia  Institute  of  Technology 
Atlanta,  Georgia 

This  research  addresses  the  unique  behavior  and  role  of  low  volatility 
ingredients  in  combustion  of  solid  propellants.  Such  ingredients  include  metal 
powders  (fuel),  ballistic  modifiers,  and  combustion  stabilizers.  Because  of 
low  volatility,  such  ingredients  tend  to  concentrate  in  varying  degree  on  the 
burning  surface,  and  their  performance  is  strongly  affected  by  processes  that 
occur  at  the  surface  in  the  concentrate.  This  aspect  of  propellant  combustion 
is  not  embodied  in  any  combustion  theory,  and  its  importance  was  not  generally 
recognized  until  its  role  in  aluminum  behavior  forced  more  detailed  study. 

Because  of  the  complexity  of  surface  behavior,  a  variety  of  strategies  is 
required  to  clarify  underlying  mechanisms,  strategies  ranging  from  thermal 
testing  of  individual  propellant  ingredients  to  combustion  of  specially 
formulated  propellants  (Fig.  1).  In  the  present  paper,  results  are  presented 
from  studies  of  edge-burning  of  "sandwiches"  (laminates  of  two  sheets  of 
oxidizer  with  a  binder  lamina  between).  Combustion  is  observed  by  high  speed 
photography  and  by  microscopic  study  of  quenched  samples  (Fig.  2).  The  low 
volatility  (powdered)  ingredient  is  contained  in  the  PB AN  binder  lamina.  This 
combustion  system  is  used  because  of  the  extensive  correlation  of  observeable 
results  and  combustion  mechanisms  that  has  been  developed  on  a  companion 
project. 

Tests  have  been  run  on  sandwiches  with  several  additive  materials  that  are 
used  as  ballistic  modifiers  or  combustion  stabilizers  in  propellants.  The 

materials  tested  can  be  classified  (Fig.  2)  in  two  categories  that  give 

contrasting  results. 

Type  1^  A^O^,  B^C,  ZrC,  CuO:  These  additives  do  not  concentrate  on 

the  burning  surface  appreciably,  and  have  little  or  no  effect  on  burning 
rate  or  quenched  burning  surfaces. 

Type  2  Fe.Oj,  CUjCrjO^:  These  additives  concentrate  on  the  burning 
surface  (binder!.  They  cause  increased  burning  rate,  and  modify  the 
sample  surface  in  a  manner  that  has  been  found  to  result  from  increased 
proximity  of  the  oxidizer-binder  flame  lets  to  the  surface.  These  two 
ingredients  are  known  catalysts  of  burning  rate  or  propellants. 

The  results  indicate  that  effective  ballistic  modifiers  must  form 

"catalyst  beds"  on  the  burning  surface,  in  which  the  large  molecular  weight 
vapors  from  the  polymer  are  broken  down  into  more  easily  oxidized  fuel 

fragments.  The  quenched  samples  indicate  that  the  oxidizer-binder  flame  stands 
correspondingly  closer  to  the  surface,  explaining  the  higher  burning  rate. 

This  interpretation  was  tested  by  addition  of  Fe^ 0.  to  a  specially 
formulated  control  propellant,  and  an  aluminized  version  or  that  propellant, 
with  results  consistent  with  the  catalyst  bed,  fuel-cracking  interpretation. 
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AERODYNAMIC  BREAKUP  -  METAL  DROPLETS 
James  E.  Craig 

SPECTRON  DEVELOPMENT  LABORATORIES,  INC. 

Costa  Mesa,  CA  92626 

The  breakup  of  AL/AL2O3  agglomerates  in  solid  propeLlant  rocket  nozzles 
is  a  critical  process  effecting  combustion  efficiency  and  two-phase  flow 
losses.  While  direct  observation  of  agglomerate  breakup  has  been  made  in 
subscale  nozzles,  the  dynamics  and  mechanisms  of  breakup  remain  to  be  charac¬ 
terized.  These  phenomena  for  high  surface  tension  metal  droplets  may  differ 
significantly  from  those  associated  with  conventional  liquid  droplets.  Proper 
scientific  scaling  techniques  must  be  developed  to  extrapolate  subscale  labor¬ 
atory  data  to  full  scale  rocket  motor  pressures,  temperatures,  and  sizes.  The 
principle  goal  of  this  investigation  is  to  provide  increased  understanding  of 
the  fundamental  breakup  mechanisms  of  high  surface  tension  liquids  in  aerody¬ 
namic  nozzles. 

Droplet  breakup  experiments  in  aerodynamic  nozzle  contractions  have 
been  conducted  in  which  conventional  liquids  and  higher  surface  tension 
liquids  (Mercury)  were  examined.  A  key  element  of  the  experiment  is  the  use 
of  laser  diagnostics.  Pulsed  laser  holography  has  provided  droplet  and  frag¬ 
ment  observations  with  resolutions  not  previously  possible.  Laser  velocimetry 
has  provided  droplet  dynamics  data  revealing  dramatic  accelerations  prior  to 
breakup. 

A  series  of  experiments  were  designed  to  investigate  the  breakup 
mechanisms  of  high  surface  tension  liquid  metals.  The  droplet  aerodynamics 
were  determined  from  trajectory  estimates  of  mercury  droplets.  For  breakup 
(We  >  30),  the  slip  velocity  requirement  ranged  from  270  m/s  to  109  m/s  for 
droplet  diameters  of  100  pm  to  600  pm,  respectively.  For  the  experiments,  the 
nozzle  was  operated  at  the  sonic  velocity,  and  the  droplet  diameter  was  varied 
from  190  pm  to  560  pm.  The  droplet  breakup  mechanism  is  depicted  in  Figure  1. 
The  droplets  are  observed  to  expand  radially  while  contracting  in  the  flow 
direction.  As  large  expansions  are  reached  (3-5X),  a  sphere  forms  in  the 
center  which  is  connected  to  the  thick  outer  ring  by  a  thin  sheet.  Much  of 
the  original  droplet  mass  accumulates  in  the  inner  sphere  and  the  outer 
ring.  Since  the  sheet  is  flat  and  thin  the  aerodynamic  pressure  forces 
quickly  expand  the  sheet  downstream  eventually  bursting  the  sheet.  These 
initial  fragments  are  very  small  and  move  rapidly  downstream  from  the  larger 
primary  fragments.  The  primary  fragment  takes  the  form  of  a  cylinder  aligned 
with  the  flow.  The  frontal  diameter  is  about  30  to  40  percent  of  the  original 
droplet  size  and  deforms  as  higher  gas  velocities  are  reached  in  the  nozzle. 
Finally,  secondary  breakup  is  observed  as  the  primary  fragment  goes  through  a 
fairly  chaotic  failure  mechanism.  Hence,  for  metal  droplets  the  breakup 
mechanism  is  cyclic  in  nature,  and  in  these  experiments  two  clearly  defined 
cycles  have  been  observed. 

The  droplet  velocity  and  Weber  number  profiles  along  the  nozzle  con¬ 
traction  were  obtained  using  double-pulse  holographic  imaging.  The  position 
(upstream  of  the  throat)  for  the  primary  breakup  event  moved  from  8  to  14  mm 
for  the  initial  droplet  sizes,  190  and  560  pm,  respectively.  The  initial 
droplets  all  failed  in  the  Weber  number  range  of  15-20.  The  Weber  number 
histories  (Figure  2)  for  all  three  droplet  sizes  overlap  when  the  time  is 
normalized  by  the  droplets  first  natural  frequency. 
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Figure  1.  High  Resolution  Pulsed  Laser  Holographic  Images  of  Mercury 
Droplet  Breakup  Mechanisms 
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DROPLET  WEBER  NUMBER  HISTORIES 
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AGGLOMERATE  BREAKUP  CONTROL 
KEVIN  K.  NACK,  ILt,  U5AF 
AFRPL/DYC,  Stop  24 
EDWARDS  AFB,  CA  93523 


Two  phase  flow  losses  are  a  significant  loss  mechanism  with  regard  to  delivered 
specific  impulse  for  solid  rocket  motors  burning  metalized  fuels.  The  primary  losses 
occur  in  the  rocket  nozzle  due  to  thermal  and  velocity  lags  of  the  metal  agglomer¬ 
ated  relative  to  the  accelerating  gas  flow  field.  If  the  size  of  the  agglomerates  can 
be  reduced,  less  loss  would  occur. 

It  is  well  known  that  some  of  the  larger  agglomerates  break  up  as  they  enter  the 
nozzle  region.  Many  of  them  do  not,  however,  and  still  cause  the  previously 
mentioned  losses.  Metal  fuel  particles  accumulate  on  the  burning  surface  of 
propellant  into  agglomerates  10  to  100+  times  the  size  of  the  fuel  particles.  Surface 
tension  of  the  molten  metal  causes  accumulation  at  the  propellant  surface  and 
inhibits  agglomerate  breakup  in  the  nozzle.  If  surface  tension  could  be  reduced, 
accumulation  could  be  inhibited  and  breakup  enhanced  both  contributing  to  reduced 
two  phase  flow  losses. 

This  project  will  use  a  unique  approach  to  test  the  critical  breakup  characteris¬ 
tics  of  molten  metal  fuels,  in  addition,  methods  of  enhancing  breakup  will  be  sought 
by  testing  breakup  of  metal/alloy  and  metal/odditive  combinations  with  more 
favorable  breakup  characteristics.  The  experiment  will  consist  of  an  acoustic 
levitator,  laser  melter,  and  shock  tube.  The  levitator  will  suspend  an  aluminum 
particle  is  a  gas,  the  laser  will  melt  the  particle,  and  the  shock  tube  will  provide  the 
flow  conditions  to  cause  breakup  of  the  molten  metal  particle. 

The  project  is  currently  in  the  build-up  phase.  Development  work  is  underway  on 
the  acoustic  levitator.  Aluminum  has  been  levitated,  but,  more  work  is  needed  to 
increase  system  reliability.  The  laser  melting  device  consists  of  a  neodymium-glass, 
pulsed  laser  and  various  optics.  The  laser  has  been  received,  but,  all  the  optics  have 
not  and  are  still  out  on  order.  The  shock  tube  is  currently  awaiting  the  construction 
of  a  support.  Experiments  on  metallic  breakup  should  be  underway  in  late  summer 
1984. 

If  the  program  demonstrates  expected  results.  Scale  rocket  motors  will  be  fired 
in  order  to  demonstrate  increased  delivered  specific  impulse. 
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Motor  #2  Exhibits  Greater  total  impulse 
t!ue  ic  Decreased  2-phase  Flow  Losses. 
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PULSED-LASER  CINEMATOGRAPHY  OF 
PROPELLANT  DEFLAGRATION 


ROGER  J.  BECKER 

UNIVERSITY  OF  DAYTON  RESEARCH  INSTITUTE 
DAYTON,  OHIO  45469 


An  improved  understanding  of  the  microscopic  behavior  of 
grains  in  deflagration  is  needed  to  model  propellant  perform¬ 
ance.  The  temporal  character  of  deflagration  on  a  local  scale 
and  its  correlation  with  the  local  environment  is  desired. 
High-speed  photography  is  a  local  method  of  obtaining  the 
desired  data,  but  conventional  cinematography  is  limited  in 
resolution  by  motion  blur,  flame  brightness,  and  illuminating 
intensity. 

Due  to  the  intensity  and  monochromaticity  of  the  laser 
pulses  used  flame  incandescence  does  not  mask  the  target. 

Motion  blur  is  negligible  due  to  the  30  ns  pulse  width.  The 
resolution  of  these  experiments  should  be  diffraction  limited. 
Consequently,  if  a  resolution  exceeding  20  ym  is  to  be  obtained, 
there  will  be  a  severe  depth  of  field  problem.  To  meet  this 
difficulty  we  have  constructed  a  rigid  frame  that  holds  the 
camera,  pressure  chamber,  and  servo  system  together.  A  seven- 
port  window  bomb  designed  for  operation  at  1,000  psi  has  been 
constructed  for  the  experiments. 

The  strand  surface  will  be  maintained  within  the  depth-of- 
field  of  an  electro-optic  servo  system.  The  servo  circuit  fea¬ 
tures  a  linear  diode  array  and  a  proportional  drive  to  protect 
against  fouling  of  the  servo  windows  and  to  maintain  a  smooth 
drive . 

All  tests  will  be  made  on  1/4  inch  strands.  The  first  set 
of  experiments  will  be  on  a  matrix  of  AP  grain  mixes  at  a  series 
of  pressures.  These  runs  will  be  followed  by  attempts  to  make 
stereo  movies,  and  studies  of  induced  quenching. 


OBJECTIVES 
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SOLID  PROPELLANT  COMBUSTION 

DIAGNOSTICS 

David  W.  Netzer 
John  P.  Powers 

Naval  Postgraduate  School 
Monterey,  California  93943 


At  the  Naval  Postgraduate  School  (NPS),  a  continuing  investigation  is 
being  conducted  to  obtain  quantitative  data  that  can  be  used  to  relate 
propellant  composition  and  operating  environment  to  the  behavior  of  solid 
particulate  (A1 ,  A1 2^3 )  within  the  grain  port  and  exhaust  nozzle. 

These  data  are  needed  in  order  to  (1)  improve  solid  propellant 
performance  predictive  capabilities,  (2)  provide  needed  input  related  to 
AP-aluminum  interactions  for  current  steady-state  combustion  models,  and  (3) 
provide  in-motor  particle  size  distributions  which  will  allow  more  accurate 
predictions  of  damping  in  stability  analyses.  The  techniques  employed  are 
high  speed  motion  pictures  of  strand  burners  and  slab  burners  in  a  cross-flow 
environment,  SEM  analysis  of  post-fire  residue  (strand,  slab,  and  motor), 
determination  of  U32  across  the  exhaust  nozzle  using  measurements  of  scattered 
laser  light,  and  holograms  of  burning  propellant  strands  and  slabs  in  a 
cross-flow  environment.  In  addition,  work  has  been  initiated  to  develop 
automatic  data  retrieval  methods  for  obtaining  particle  size  distributions 
from  holograms  taken  of  the  combustion  of  solid  propellants.  Once  developed, 
these  techniques/diagnostic  methods  could  be  readily  employed  for  obtaining 
the  needed  data  discussed  above  from  a  series  of  tests  in  which  propellant 
composition,  motor  geometry,  and  operating  environment  are  systematically 
varied. 

Actual  particle  sizes  of  burning  aluminum  particles  were  obtained  in  high 
speed  motion  pictures  by  using  high  intensity  rear  illumination  of  the  burning 
propellants  to  eliminate  the  flame  envelopes  surrounding  the  burning 
particles. 

Measurements  of  di ffracti vely  scattered  light  were  successfully  made  for 
determination  of  changes  in  U32  across  a  solid  propellant  rocket  motor  exhaust 
nozzle.  However,  to  date,  propellant  composition  and  grain  geometry  have 
significantly  limited  the  range  of  obtainable  data. 

Two-dimensional  motors  have  been  successfully  employed  to  obtain  good 
quality  holograms  of  propellant  burned  in  a  cross-flow  environment.  Other 
efforts  have  included  reduction  of  speckle  in  the  recorded  hologram  and 
optimization  of  techniques  for  minimizing  excessive  smoke. 

A  Quantimet  720  has  been  used  in  an  initial  effort  to  obtain  particle 
size  distributions  from  a  photograph  taken  of  a  reconstructed  hologram. 

Current  efforts  are  directed  at  using  a  computer-controlled  Quantimet  720  and 
a  computer-controlled  x-y-z  translation  stage  to  obtain  particle  size  data 
directly  from  the  reconstructed  hologram. 


ally  formulated  propellants:  HTPB/AP/{A1,  ZrC,  or  G),  GAP/AP/A1 


Figure  2.  PRIMARY  ACCOMPLISHMENTS 
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FUEL -RICH  SOLID  PROPELLAOT  BORON  COMBUSTION  (ANALYSIS) 

Merrill  K.  King 
Atlantic  Research  Corporation 
Alexandria,  VA 

The  analytical  modeling  tasks  are  directed  at  evaluating  hypotheses 
regarding  mechanisms  of  Ignition  and  combustion  of  boron  particles,  agglomerates, 
and  clouds,  and  delineating  how  they  affect  combustor  characteristics.  The 
uniqueness  of  the  approach  lies  mainly  In  the  development  of  detailed  analytical 
descriptions  of  unit  processes  with  subsequent  "building-block"  combination  for 
definition  of  means  of  maximizing  use  of  boron's  high  heating  potential.  Emphasis 
is  placed  on  the  fact  that  boron  presents  unusual  Ignition/combustion  problems 
demanding  careful  attention  to  environmental  staging. 

A  model  of  boron  Ignition  In  wet  or  dry  atmospheres  has  been  developed: 
this  model  Is  unique  In  its  treatment  of  oxide-coated  boron  oxidation  as  Involving 
series  resistances  associated  with  oxygen  gas  diffusion,  dissolution  of  oxygen  Into 
the  oxide  film,  subsequent  solute  diffusion  to  the  metal/oxide  Interface,  and 
finite-kinetics  reaction  at  that  interface.  Evaporation  of  BpO,  and  reaction  of 
this  oxide  with  HpO  to  form  HBOp(g)  are  both  modeled  as  involving  both  kinetic  and 
diffusional  resistances.  Ignition-time  predictions  with  this  model  are  found  to 
agree  well  with  the  data  base  of  Macek.  In  addition,  a  quasi-steady  model  of 
clean-particle  combustion  with  COp,  Op,  and/or  HpO  with  a  unique  finite- 
rate-kinetics  treatment  of  the  surface  reactions  of  boron  with  these  oxidizers  (and 
other  Intermediates  with  0/B  ratios  >  1)  has  been  developed  and  calibrated  against 
Macek  data.  This  model  has  been  extended  to  treat  unsteady-state  effects  and 
coupled  with  the  Ignition  model  to  yield  a  unified  ignltlon/combustlon  model.  A 
unique  aspect  of  this  model  Is  definition  of  surface  gas  composition  and  temperature 
versus  time,  with  consequent  Identification  of  extinguishment  conditions  (occurring 
when  PB203>S  >  ^203,3 5  * 

Recently  the  Ignition  model  has  been  extended  to  treat,  possible  pfffWs  of 
halogens,  with  a  first-ever  examination  of  reactions  of  e  C lp,  HE,  and  HC1  wl^h 
the  oxide  layer.  Also,  a  mechanism  by  which  an  outer  layer  or  titanium  or  zirconium 
might  aid  Ignition  by  providing  particle  heating  without  production  of  additional 
boron  oxide  to  a  point  where  the  original  oxide  flashes  off  or  evaporates  quickly 
following  exhaustion  of  the  second  metal,  has  been  modeled.  Finally,  a  'inlque  ^odel 
of  the  effects  of  LIP  coatings  Involving  the  reaction  of  L1F  with  t0  corr’ 

and  BOP  has  been  developed. 

Two  models  of  a  boron  slurry  ramjet  combustor  uniquely  considering  th*1 
special  problems  associated  with  boron  Ignltlon/combustlon  through  Incorporation  of 
the  above  analyses  have  been  developed.  Both  Include  a  well-stirred  reactor  P\ 
Ignition  zone  In  which  the  fuel  and  part  of  the  air  are  mixed  and  partially  reacted. 
In  the  first  model,  downstream  addlti jn  of  the  remaining  air  Is  treated 
one-d linens lonally,  with  user-input  air  addition  schedule,  while  In  the  second, 
two-dimensional  coaxial  mixing  Is  treated. 

Major  accomplishments  of  this  effort  Include  Identification  or  several 
approaches  to  enhancing  boron  particle  Ignition  (Pig.  2  -  Item  I),  definition  of 
approaches  to  proper  boron  slurry  combustor  design  (’Fig.  ?  -  Item  II),  and 
delineation  of  causes  for  observed  deleterious  effects  of  low  combustor  operating 
pressure  on  boron-powered  airbreathing  systems  (Fig.  ?  -  Item  III). 
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I.  REACTIONS  TREATED  IN  MODELING  OF  SINCLE  BORON  PARTICLE  IGNITION  AND  COMBUSTION 


II. 


A.  IGNITION  IN  B-O-H-N-CI-F  ATMOSPHERES 
S02  +  2B(c)  -*  B203(1)  +  Q 

H20  +  B,03(l)-»  2HB02/2/3H3B306  +  Q 
B203(1)  -»  B203(g)  -  Q 

H20  +  4/3S(c)  +  1/JB203(1)~  2/3H3B303  +  Q 
b2°3(1)  +  Ci2—  SP2  +  2BOCl/2/3B3Cl303  +  Q 
B203(1)  +3F2  -•  3/202  +  2BFj  +  Q 

B203(1)  +  HC1  -*  B0C1/1/3B3C1303  +  HB02/l/3H.jB305  +  Q 
BjOjCl)  +  HF  — •  BOF  +  HB02/1/3H3B306  +  Q 

B.  IGNITION  WITH  COATINGS 


C. 


Tl(s)  +  02  —  T102(s)  +  Q 
Zr(s)  +  02  —  Zr02(s)  +  Q 
LiF(l)  +  B203(1)  —  LiB02(l)  +  BOF  -  Q 


CLEAN  PARTICLE  COMBUSTION 

B(c)  +  02,B02,H20,B203(g),HB02,C02  -* 


EQUILIBRIUM  GAS  MIX  OF  B2°3  +  0 

SUBOX  IDES  AND  HYDROXIDES  -  1  „„„ 

H,r  HBO-  +  Q 


SLURRY  COMBUSTOR  ANALYSES 


ZERO-DIMENSIONAL  ANALYSIS 

(1)  UTILIZE  MASS, ENERGY .SPECIES  BALANCES  OVER  IGNITION  ZONE 

(2)  CALCULATE  RESIDENCE  TIME  DISTRIBUTION  FROM  P  (TR  »  0  )  -  EXP<-*/TR) 

(3)  ANALYZE  IGNITION/COMBUSTION  OF  I  x  J  BINS  OF  PART ICLE-SIZE/RESiDENCE-TIME  COMBINATIONS 

(4)  USE  FL..L-UP  BORON  IGNITION  MODEL  AND  CORRELARIONS  BASED  ON  DETAILED  FINITE-KINETICS 
COMBUSTION  MODEL  TO  CALCULATE  FRACTION  VO RON  IN  EACH  BIN  IGNITED  AND  BURNED 

(5)  ITERATIVE  PROCEDURE  TO  CLOSE  LOOP  ON  TOTAL  FRACTION  BORON  BURNED  AND  IZ  TEMPERATURE,' 'CoMI'osl  I  ION 

ONE-DIMENSIONAL  ANALYSIS 

(1)  SPECIFIED  AIR  ADDITION  TO  EACH  AXIAL  INCREMENT 

(2)  OVERALLMASS, ENERGY, SPECIES  BALANCES  TO  CALCULATE  TEMPERATURE,  COMPOS  I TIoN,  AND  RESIDENCE 
TIME  IN  THE  INCREMENT 

(3)  IGNITION/COMBUSTION  MODELS  USED  TO  CALCULATE  INCREMENTAL  BORON  IGNITION  ANU  COMBUSTION  IN 

EACH  BIN 

TWO-DIMENSIONAL  ANALYSIS 

(1)  CORE-LENCTH,  CENTERLINE  DECAY  CORRELATIONS  USED  TO  DEFINE  VELOCITY .  ENTHALPY,  ATOMIC 
COMPOSITION,  AND  TRANSPORT  PROPERTY  DISTRIBUTIONS 

(2)  MASS  TRANSPORT  EOL  ATION  AND  CONSERVATION  EQUATION  FOR  BORON  OXIDE  INCLUDING  SOURCE  TERM  I  OR 
REACTION  USED  TO  DETERMINE  MOLECULAR  COMPOSITION  OK  I  ACH  SUCCESSIVE  RADI AL/AXI AL  ANNULAR  MU! 

(3)  ENTHALPY  EQUATION  USED  TO  CALCULATE  TEMPERATURE,  RESIDENCE  I  IMF  I  OR  EACH  SLICE 

(4)  IGNtTION/COMBUST  ION  MODELS  USED  TO  CALCULATE  BORON  CONSUMPTION  UnR  LACK  SLICE 
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IDENTIFICATION  OF  KEANS  OF  DECREASING  BORON  PARTICLE  IGNITION  DELAYS 


10  Micron  Particle  Radius,  0.2  Micron 
Initial  Oxide  Thickness,  P  *  1  Atm, 
Initial  Temperature  *  300*K 

1  Base  Case  P^g-O. 10,  *  0.  10 

2  PH20-  0.00,  P02  -  0.20 

3  PH20  -  0.20,  P02  *  0.00 

1  Base  +  Pf2  "  0.10 

3  Base  +  PHf.  *0.10 

6  Base  +  Pci2  ^ 

7  Base  +  ?HC1  ”0-10 

8  Titanium  Coating  (Ti/B-*0.15) 

9  Zirconium  Coating  (Zr/B— 0.25) 

10  LiF  Coating  (LiF/B** 0. 08) 


IDENTIFICATION  OF  IMPORTANCE  OK  TAILORING  AIR  DISTRIBUTION  IN  BORON  SLURRY  COMBUSTORS 


PRESSURE  (ATM) 


ftt truaa  i*rm) 


PRESSURE  1  *  run 


*» 
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FUEL-RICH  SOLID  PROPELLANT 
BORON  COMBUSTION  (Experimental) 

James  J.  Komar  and  Ronald  S.  Fry 
Atlantic  Research  Corporation 
Alexandria,  Virginia  22314 

This  experimental  research  is  aimed  3t  understanding  the  basic 
physics  of  boron  particle  ignition  and  combustion.  Experimental  data  are 
acquired  on  the  ignition  and  combustion  of  boron  particles,  agglomerates 
and  clouds,  as  well  as  combustion  and  ignition  phenomena  in  consolidated 
fuel  grains  in  high  temperature  air  crossflow.  This  research  is  vital  to 
identifying  means  of  using  the  full  potential  of  boron  fuels. 

Current  efforts  in  the  single  particle  combustion  investigations  are 
being  directed  towards  the  measurement  of  boron-oxide  reaction  rates, 
utilizing  small  single  beads  of  boron  oxide  mounted  on  micro-thermocouple 
substrates.  Determination  of  size  changes  as  functions  of  time,  operating 
conditions  and  initial  particle  sizes  are  being  made  with  microoptical 
pyrometry,  diode  array  photography,  particle  size  Interferometry,  and 
photographic  recording.  These  reaction  rates  are  not  well  known  and  the 
simultaneous  diagnostics  approach  is  unique  in  that  self-consistent 
measurements  are  possible  during  the  course  of  the  experiments. 
Representative  scanning  electron  micrographs  are  shown  of  oxide  particles 
of  varying  stages  of  reaction  completion. 

Combustion  of  consolidated  boron  fuel  grains  in  a  high  temperature 
air  crossflow  involves  stabilizing  a  flame  within  the  turbulent  boundary 
layer  above  the  ablating  grain  surface.  Mechanisms  fundamental  to  the 
problem  Include  grain  surface  ablation  characteristics,  ablated  material 
physical  properties  and  dynamics,  radiation  feed-back  to  the  fuel  surface 
from  burning  particles  in  the  main  stream  and  mean  and  fluctuating 
properties  of  the  main  stream  and  boundary  layer  regions. 

Investigation  of  these  mechanisms  are  conducted  in  a  unique 
two-dimensional  windowed  combustor  whose  design  comprises  two  facing  slab 
grains  to  allow  examination  of  radiative  coupling.  Flame  structure 
dynamics  are  studied  with  high  framing  rate  color  and  Rchlieren 
photography.  Particle  size  interferometry  is  being  employed  in  an  effort 
to  characterize  particle  size  distributions  of  the  ablated  material 
released  from  the  fuel  surface.  The  radiation  thermal  distribution  and  its 
interaction  with  the  convective  flow  above  the  fuel  surface  is  studied 
using  a  Retlcon  diode  array  camera.  Details  of  the  fuel  grain  conductive 
thermal  distribution  and  ablation  behavior  Is  studied  with  embedded  high 
response  thermocouples  and  fiber  optic  probes.  One  component  laser 
veloclmetry  is  employed  to  characterize  the  mean  flow  field. 

During  this  period,  the  experimental  apparatus  and  diagnostics  were 
checked  out  through  air-only  flow  and  hydrocarbon  fuel  combustion  testing. 
Study  of  the  combustion  details  of  boron-containing  fuels  is  ongoing. 
Preliminary  criteria  for  flame  stability  in  2-D  combustor  were  established. 
Definition  of  overall  flow  structure  is  possible  from  available  data. 
Detailed  features  of  the  combustion  dynamics,  identified  as  program  needs, 
are  being  sought. 
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BORON  COMBUSTION  DYNAMICS  APPROACH 

Boron  Oxide  Kinetics 


+  8,0,(l)  -  2  H80,(|) 


Boron  Grain  Flame  Structure 
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Seep 


Turbulent  Boundary  Layer 
Immersed  Flame  Structure 

Ablating  Grain 


APPROACH 


Phyalot  of  boron  coabuaelon  In  high 
taaperetora  In  croaa-flow  ara  not  undaratood. 

Hodallng  tha  coaplax  Intaractlon  batwaan 
blnatlea  and  alxlng  llaltad  phanoaana  la 
Intractabla  without  axpariaantal  data 

Prior  raaaareh  aeployad  dlagnoatlca  lnadaquata 
for  obtaining  aaaantlal  data. 


Unlqua  eoablnatloo  of  atperlaapt  daalgn  and 
dlagnoatlca  to  aaaaloa  boron  coabuatlon  in 
high  taaparatura  air  croaa  flow. 

Laboratory  acala  coabuatlon  tatting  In 
unlqual?  aultad  2 -Dlntna lonal  wlndowad 
coabuator. 

Application  of  advancad  non-lnt rual *a  optical 
dlagnoatlca  and  data  procaaalng. 
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OVERVIEW:  AFOSR  INTERESTS  IN  ROCKET  PROPULSION 
Leonard  H  Caveny 

Air  Force  Office  of  Scientific  Research 
Bolling  A  KB 
Wash  i  rig  tori ,  DC  ,'Oi'L1 

Present  and  future  rocket  propulsion  systems  present  critical  challenges  that 
are  not  addressed  adequately  by  present  technologies.  Propulsion  systems  must 
undergo  continuous  advancements  to  assure  that  the  winning  edge  is  maintained. 

As  the  conventional  propulsion  technology  matures,  the  need  persists  to  attain 
yet  another  round  of  performance  gains.  Thus,  the  technology  challenges  become 
even  greater  since  the  obvious  advancements  have  been  implemented.  Also,  broad 
classes  of  problems  (e.g.,  combustion  instability,  service  life,  signature, 
low- temperature  structural  integrity,  light  weight  thermal  protection)  are  never 
solved  entirely  and  recur  often  as  important  considerations  (and  compromises) 
prior  to  propulsion  system  qualification. 

Powerful  motivations  continue  for  basic  research  on  topics  such  as  reacting 
flows,  energetic  materials,  stability,  refractory  materials,  propulsive 
processes.  Many  phenomena,  not  explained  adequately,  are  accommodated  by 
specific  designs.  Periodically,  this  lack  of  understanding  is  the  precursor  to 
major  setbacks.  Clearly,  rocket  motor  processes  present  scientific  challenges 
suggesting  that  properly  conceived  research  will  lead  to  higher  performance, 
reduced  risk,  lower  development  costs,  longer  service  life,  etc. 

This  meeting  will  be  the  first  review  of  the  research  relating  to 
nonconventional  primary  propulsion  for  space.  Accordingly,  some  introductory 
discussion  will  help  put  portions  of  this  research  into  perspective. 
Nonconventional  propulsion  concepts  periodically  receive  attention  but 
invariably  the  propulsion  systems  chosen  have  been  adaptations  or  extensions  of 
conventional  liquid  and  solid  rocket  technology.  The  dominant  consideration  i  ri 
previous  years  was  that  the  missions  could  be  performed  using  conventional 
chemical  propulsion.  Consequently,  major  initiatives  to  provide  technology  and 
to  overcome  specific  barriers  could  not  be  justified.  The  recent  emphasis  on 
space  and  the  reusable  launch  vehicle  capability  to  low  Earth  orbit  present  new 
considerations  for  advanced  propulsion  for  orbit  transfer.  Transferring  the 
projected  large  payloads  using  conventional  propulsion  would  impose  severe 
limitations  on  the  missions.  For  example,  75$  of  the  mass  delivered  to  low 
Earth  orbit  may  be  the  chemical  propulsion  system  required  to  raise  the  other 
25$  (i.e.  the  active  payload)  to  geosynchronous  Earth  orbit;  nonconventional 
propulsion  offers  the  promise  of  reversing  this  ratio  of  propulsion  and  payload 
masses.  The  technology  issues  are  necessarily  complex;  controversial  approaches 
are  to  be  expected.  Much  of  the  research  is  coupled  at  optimistic  projections 
for  technological  advances  by  the  1995  to  2000  time  period. 

The  Wednesday  morning  session  treats  fundamental  processes  related  to  beamed 
energy  and  solar  propulsion,  i.e.,  radiation  transmission  and  absorption  in 
flowing  media,  plasma  initiation,  and  chamber  configurations.  The  immediate 
advantage  of  beamed  and  solar  energy  is  high  specific  impulses  obtained  by  using 
low  molecular  weight  working  fluids  heated  by  external  power  sources  to 
temperatures  greatly  above  combustion  gas  temperatures.  The  premise  of  the 
beamed  energy  approach  is  that  suitable  megawatt  laser  sources  will  be 
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available  and  justified  for  applications  other  than  propulsion.  Recent  studies 
support  the  thesis  that  ground  based  free-electron  lasers  operating  at  suitable 
short  wave  lengths  will  enable  attractive  transmission  efficiencies  and  curm-vt 
collection  optics.  Conversely,  a  scenario  for  using  orbiting  loser  sy.  t-m.-. 
beamed  energy  propulsion  is  considered  to  be  premature. 

The  Wednesday  afternoon  session  on  electric  propulsion  emphasizes  concepts  w'n.b 
lend  themselves  to  sustained  operation  at  megawatt  power  levels,  possibly  using 
clusters  of  thrusters.  Thi3  represents  a  major  departure  from  millinewtm 
thrust,  pulse-mode  electric  propulsion  considered  for  station  keeping  or 
planetary  probes,  e.g. ,  millisecond  pulses  driven  by  a  capacitor  bank  charged  by 
a  dedicated  kilowatt  power  source.  Present  mission  analyses  are  exploring  the 
dual  mode  premise  of  the  megawatt  nuclear  power  supply  onboard  for  the  main 
mission  being  available  for  propulsion  power.  Thus,  propulsion  does  not  take 
all  the  weight  penalty  for  the  power  source.  When  large  total  impulses  are 
required,  thruster  mass  is  small  compared  to  the  fuel  mass;  thus,  a  premium  is 
placed  on  increasing  fuel  efficiencies  and  the  low  thrust  densities  of  the 
thrusters  become  unimportant.  Electrode  and  insulator  lifetimes  have  beer, 
identified  as  primary  barriers  to  sustained,  high  power  density  operation  of 
magnetoplasraadynamic  thrusters.  Accordingly  a  major  portion  of  the  research  is 
directed  at  mass  loss  mechanisms  and  conditions  leading  to  abusive  environments 
and  inefficient  operation. 

The  Thursday  morning  session  includes  presentations  on  several  of  the  more 
challenging  and  higher  payoff  approaches  to  space  propulsion.  Such  a  forum 
stimulates  thinking  on  the  basic  research  approaches  which  must  accompany  major 
advancements  in  space  propulsion. 

Ample  opportunity  exists  for  new  approaches.  During  the  last  two  years,  several 
investigations  into  the  synthesis  of  energetic  ingredients  were  initiated; 
complementary  research  is  needed  to  improve  the  characterization  techniques  for 
these  higher  energy  (and  more  temperamental)  ingredients  and  propellants,  dood 
progress  is  being  made  on  understanding  the  origins  of  combustion  instability; 
more  attention  must  be  given  to  the  mechanisms  of  deliberate  suppression.  The 
advent  of  quantitative  flow  visualization  for  turbulent  reacting  flows  ms  not 
been  accompanied  by  corresponding  theoretical  treatments  to  fully  exploit  array 
data,  capable  of  revealing  rapid  evolution  of  flow  structure  and  flame  fronts. 
The  remarkable  advances  in  optical  diagnostic  techniques  present  new  approaches 
to  investigate  plasma  flows  which  must  be  understood  and  controlled  under 
magnetoplasmadynamic  and  beamed  energy  thruster  conditions.  The  research  on 
life  limiting  processes  which  occur  at  the  surfaces  of  electrodes,  electrical 
insulators,  and  refractories  will  benefit  from  advances  in  remote  sensing  of 
surface  temperatures,  composition,  and  structure.  Space  -issions  will  present 
entirely  new  autonomous  operation  challenges;  a  theoretical  "asis  must  be 
established  to  anticipate  and  guide  the  advances  in  sensors,  adaptive  control, 
thruster  configuratin  options,  etc.  Sustained  megawatt  operation  of  space 
thrusters  will  be  limited  by  the  inability  to  reject  heat;  research  is  needed  to 
enable  advances  in  light  weight  radiators  and  thermal  management  systems.  We 
will  welcome  discussion  on  these  and  other  topics  pertinent  to  the  Air  Force 
basic  research  program  on  energy  conversion. 


** 
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EVALUATION  AND  COMPILATION  OF  THE  THERMODYNAMIC 
PROPERTIES  OF  HIGH  TEMPERATURE  SPECIES 


Malcolm  W.  Chase 
The  Dow  Chemical  Company 
Thermal  Lab 

Midland,  Michigan  48640 


The  JANAF  Thermochemical  Tables  are  a  set  of  self-consistent  thermodynami c 
data.  These  tables  result  from  critically  reviewing  all  literature 
sources,  evaluating  the  accuracy  and  precision  of  the  experimental 
data/theoretical  calculations,  calculating  temperature-dependent 
thermodynami c  functions,  and  publishing  thermochemical  tables  for  use  in 
the  Air  Force  community. 

The  successful  application  of  these  tables  depends  on  the  user  being  aware 
of  the  uncertainties  in  the  numerical  values  and  the  possible  physical  and 
chemical  phenomena  peculiar  to  the  species  of  interest.  The  quality  of  the 
tabulation  is  often  temperature-dependent  and  is  a  compromise  between  the 
quality  and  extent  of  the  available  data  and  the  calculational  procedures 
used  to  generate  the  thermochemical  tables.  The  effort  directed  towards 
preparing  annotated  bibliographies  and  data  graphs  is  intended  to  aid  the 
users  in  judging  the  effects  of  the  uncertainties  on  the  end  results. 

The  current  contract  involves  the  re-examination  of  the  alkaline  earth 
metals,  their  gaseous  dimers,  oxides,  hydroxides,  halides,  carbonates, 
sulfides,  and  sulfates.  Bibliographies  and  data  graphs  are  being  generated 
for  each  of  the  75  species;  this  will  result  in  approximately  175  single- 
and  multi-phase  thermochemical  tables.  This  analysis  is  revealinq  a 
significant  lack  of  information  and  some  potential  problem  areas.  Three  of 
the  five  alkaline  earth  metals  are  not  well  cha racteri zed 
thermodynamically.  The  attached  graph  for  the  low  temperature  heat 
capacity  of  Mg,  Ca  and  Sr  illustrates  the  suspect  heat  capacity  of  calcium. 
A  similarity  in  the  temperature  dependence  of  these  three  metals  is 
expected.  A  similar  comparison  of  the  heat  capacities  of  the  alkaline 
earth  dihalides  suggests  that  the  CaCl2  data  may  be  in  error.  The  alkaline 
earth  metal  dimers  (gaseous)  pose  a  different  problem  in  that  the 
calculational  methods  must  be  carefully  chosen  to  properly  calculate  the 
thermal  functions  (see  attached  graph).  The  temperature  dependence  of  the 
heat  capacity  of  Ca2(q)  and  Be2(g)  would  be  expected  to  be  similar  to 
Mqx(q). 

The  examination  of  experimental/theoretical  data  by  itself  may  not  be 
sufficient  to  indicate  potential  problems  whereas  the  examination  of 
families  of  species  often  reveals  interesting  problems  and  discrepancies. 
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CRITICAL  EVALUATION  OF  HIGH  TEMPERATURE  CHEMICAL  KINETIC  DATA 


Norman  Cohen  and  Karl  Westberg 
Aerospace  Corporation 
P.O.  Box  92957 

Los  Angeles,  California  90009 


The  rapid  growth  in  chemical  kinetic  data,  and  simultaneously  in  the 
number  of  nonspecialists  who  require  such  data  in  DOD-related  work,  has  cre¬ 
ated  a  need  for  a  reliable,  easy-to-use  compilation  of  evaluated  rate  data  and 
recommended  rate  coefficients.  To  answer  this  need,  a  program  for  the  evalua¬ 
tion  of  kinetic  data  and  for  the  preparation  of  consistently  formatted  data 
sheets  patterned  after  the  JANAF  Thermochemical  Tables  was  undertaken  with  the 
joint  support  of  the  Air  Force  Office  of  Scientific  Research  and  the  National 
Bureau  of  Standards.  One  such  data  sheet  is  shown  in  Figs.  1  and  2. 

High-temperature  data  are  of  special  interest  to  the  Air  Force,  and  this 
is  the  first  data-evaluation  program  that  routinely  extrapolates  rate  coeffi¬ 
cients  for  bimolecular  reactions  using  methods  more  sophisticated  than  the 
familiar  Arrhenius  expression.  Theory  is  also  used  to  predict  rate  co¬ 
efficients  for  important  reactions  for  which  there  are  no  experimental  data. 

To  date  35  sheets  have  been  prepared  for  reactions  of  interest  of  various 
Air  Force  programs  and  published  in  a  technical  report.  These  include  six 
reactions  important  in  hydrogen-oxygen  combustion;  seven  in  hydrogen  halide 
chemical  lasers,  fifteen  in  hydrocarbon  oxidation,  five  in  aluminum  or  boron 
propellant  systems;  and  two  in  the  oxygen-iodine  chemical  laser.  Twenty-seven 
of  these  have  recently  been  published  in  J.  Phys.  Chem.  Ref,  Data. 

A  series  of  data  sheets  for  reactions  of  0  atoms  with  alkanes  (methane, 
ethane,  propane,  butane.  Isobutane,  cyclopentane,  pentane,  neopentane,  hexane, 
cyclohexane,  heptane,  cyclohepfane,  octane,  2,3-diraethylbutane,  neooctane)  has 
just  been  completed  and  this  work  will  be  discussed.  The  experimental  data 
have  been  reviewed  and  transition  state  theory  calculations  carried  out  for 
extrapolating  the  data  to  higher  temperatures.  Recurrent  problems  with  the 
experimental  data  for  CH^,  C2Hg,  neo-C^H^o*  and  neo-CgHjg  include  the  effects 
of  trace  impurities  in  the  hydrocarbon  and  the  effects  on  the  calculated  rate 
coefficient  of  secondary  reactions.  The  importance  of  these  effects  is  a 
result  of  the  primary  reaction  being  so  slow.  The  theoretical  approach  for 
reactions  with  such  large  molecules  is  quite  different  from  what  can  be  done 
with  simple  reactions,  such  as  F  +  Hj ,  for  which  fairly  reliable  a  priori 
potential  energy  surfaces  are  available. 

During  the  remainder  of  the  fiscal  year,  we  will  complete  data  sheets  for 
reactions  of  0  and  OH  with  NHg  and  N2H4,  and  for  reactions  of  OH  with  a  series 
of  halomethanes. 


** 


or  ♦  (c« j)4e 


** 

•CHjC(CHj)j  +  HjO 
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IhEWiGCHEMSlRY 

Tharmochamlcal  date  for  H20  ere  taken  from  the  second  edition  of  JANAP  Tharmochamlcal  labia*  (1971);  data  (or  ON  are  fio* 
unpublished  supplement  to  those  Cables  dated  Juno  1977.  Data  for  (CH^)^C  { 2, 2-dlmechyl propane,  or  noopontsno)  art  taken  '.</■  Stvii 
•e  al.J  oJ^neopantyl-H)  1*  reported2  aa  419.25  *  4.2  kJ  mol"1,  from  which  can  b«  calculacad  co  b*  35.7  kj  -o  "l  tor  the 

naopentyl  radical.  for  cha  radical  waa  calculacad2  to  be  329.4  J  mol ~ lK~ 1 . 

An  analytic  axpraaalon  for  R(  T)  haa  not  boon  calculatad  bacauaa  cha  ravaraa  reaction  will  navat  ba  important;  mi  1 
temperatures  it  la  too  endothermic  co  occur,  and  at  high  temperatures  cha  naopantyl  radical*  ara  unstable. 


MEASUREMENTS 

Pour  experimental  measurements  of  kj  hawa  baan  raportad.  Greiner*  made  obaarvatlona  ov*r  tha  temperature  rang*  of  292-493  K 
using  flash  photolysis  of  HjO  aa  an  OH  aourca,  monitoring  OH  disappearance  b y  kinetic  abaorptlon  apactroacopy .  JVo  room  tamparsrura 
measurements3, *  vara  aada  both  ualng  flaah  photolyala  of  H20  and  raaonanca  fluorescence  to  aonltor  OH  dlaappaarancc.  In  tha  cats  of 
ona  of  these,®  tha  measurement  waa  aada  ralatlva  to  raactloo  batwaan  OH  and  a-CgbjQ.  Ualng  our  own  recoamended  rata  coefficient  for 
tha  lattar  raactloo  (which  la  about  10%  aaallar  than  chat  uaad  In  Raf.  6),  wa  obtain  kj  •  5.7  a  2  0®  i  mol “ 1  a~ 1 .  Tha  thra*  room 

temperature  measurements  than  agraa  within  battar  than  10%;  tha  avaraga  valua  la  5.5  t  0.3  a  1U®  L  mol"1#"1.  Bakar  at  al.7 8  measured 

kj/k2,  vhara  k2  la  tha  rata  coafflclant  for  OH  ♦  H2  ♦  H  ♦  H20,  by  adding  traca  amount*  of  naopantana  to  Hj-Oj  mixtures  at  7  5)  K. 
Ualng  our  own  recommended  valua  for  k2  (aaa  Data  Shaat  for  that  raactlon)  wa  obtain  kj  -  3.2  a  10*  L  mol-1i“ * . 

CALCULATIONS 

Transition  atata  thaory  calculatlona  wara  carrlad  out  for  kj  In  ordar  to  extrapolate  tha  rata  coafflclant  ba’.rrn*  tha 
temperature  ranga  of  tha  axparlaantal  data.®  S*(298)  waa  calculatad  to  ba  388.3  J  mol"1  K"1  ualng  naopantana  aa  a  aodaj  compound, 
giving  AS®  “  -102  J  mol"1*"1.  Vibrational  frequencies  wara  takan  to  ba  tha  aaaa  aa  for  naopantana  axcapt  for  tha  removal  of  a  Oh 
•tretch  (3100  co'1).  ao  HCH  band  (1400)  and  a  CCH  band  (1100);  and  tha  addition  of  tha  following  frequencies:  3700  (OH  stretch). 
2200  (C*H)  atratch;  1000,  1000,  350  and  300  (deformations),  Tha  two  naw  hlndarad  lntarnal  rotora  wara  aaaiaud  to  hav«  fare' ctnn 
function*  OgOOO)  -  5.6  and  3.9,  and  barrlara  to  rotation  V  •  8.4  and  4.2  kJ  aol'1,  raapactlvaly .  Tha  alactronlc  degeneracy  of  the 
activated  complex,  g®,  waa  assumed  to  ba  2.  Tha  raaultlng  valuaa  of  kt(T)  ara  wall-daacrlbad  by  tha  expression  kj  •  6.3  *  102  T2*0 

axp(-10/T)  L  mol"1#"!.  Thia  axpraaalon,  darlvad  ualng  tha  axparlaantal  valua  of  kjUOO)  -  5.5  x  10®,  la  In  goed  agreement  with  the 

axparlaantal  data  at  higher  fmpmrmturea. 


DISCUSSION 

Thla  raactlon  la  of  thaoratlcal  lntaraat  bacauaa  naopantana  contains  only  ona  kind  of  h  a  ton  (prlaary),  so  chat  #h<~<  \,  -inly 
ona  reaction  Caking  place,  In  contrast  with  straight-chain  alkane*  with  mere  than  2  carbon  itoat. 

Tha  experimental  data  are  in  good  accord;  tha  discrepancy  among  tha  thraa  rooai  temperature  meaeutar  «nta  la  t.oc  a  ms  1  tc 
any  real  problems.  Since  our  calculatad  axpraaalon  agraaa  with  all  tha  axparlaantal  data  within  experiments*  unct.-r.tlni>,  vs*  c 
for  the  recommended  axpraaalon.  lhua  wa  recommend  kj  -  6.3  x  103  T2,0  «xp(-l0/T)  l  mol"1#'1,  with  an  uncertainty  m  log  kj  of  t  0. I 
at  300  X,  Increasing  co  ±0.5  at  2000  X. 


1.  D»  R.  Stull,  E.  F.  Uestna,  Jr.,  and  C.  C.  Sinks,  Tha  Qi  ami cal  The  rmod>n  sales  of  organic  Compound#  (Wiley,  New  York,  19b9), 
p.  247. 

2.  C.  W.  Larson,  E.  A.  Kardwldge,  and  8.  S.  Rablnovltch,  J.  Cham.  Phy*.  50 ,  2769  (1969). 

3.  H.  E.  O'Neal  and  S.  W.  Banaon,  " Tharmochamlat ry  of  Free  Radicals,'*  In  Free  Radicals,  ad.  J  K .  kochl  : Wiley,  hew  Y^-h,  1473), 

p.  275. 

4.  N.  R.  Crelnar,  J.  Cham.  Phya.  ^3,  1070  (1970). 

5.  C.  Paraakavopouloa  and  W.  S.  Nip,  Can.  J.  Chem.  58,  2146  (1980). 

6.  K.  R.  Darnell,  R.  Atklnaon,  and  J.  N.  Pitta,  Jr.,  J.  Phya.  Cham.  82,  1581  (1978). 

7.  X.  R.  Biker,  X.  X.  Baldwin,  and  X.  V.  Walker,  Coabuae.  Flame  Vj  Xkl  <^76). 

8.  N.  Cohan,  Int.  J.  Cham.  Klnet.  _14,  1339  (1982). 
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LASER  THERMAL  PROPULSION 


Dennis  Keefer 

The  University  of  Tennessee  Space  Institute 
Tullahoma,  Tennessee  37388 


The  principal  objective  of  this  research  investigation  is  to 
determine  experimentally  the  effects  of  a  forced  convection  environ¬ 
ment  and  optical  geometry  on  the  stability,  fractional  power 
absorption,  plasma  structure,  and  fluid  mixing  in  a  laser  sustained 
plasma  (LSP).  A  continuous,  1.5  kW,  axial  flow,  carbon  dioxide 
laser  is  used  to  create  the  LSP  in  a  cylindrical  quartz  flow  channel. 
The  convection  flowfield  surrounding  the  plasma  is  controlled  by  the 
volume  flow  through  the  test  chamber,  and  the  optical  geometry  is 
determined  by  the  focal  length  of  the  lens.  Digital  images  of  the 
plasma  in  selected  narrow  wavelength  intervals  are  obtained  using  a 
calibrated,  CID  digital  camera  and  a  VICOM  digital  image  processing 
computer.  These  images  are  then  Abel  inverted  to  give  a  spatial 
plasma  emission  coefficient  which  is  used  to  determine  the  spatial 
distribution  of  the.  plasma  temperature.  The  measured  temperature 
field  of  the  plasma  is  used  to  determine  the  spatial  distribution 
of  the  plasma  absorption  coefficient  for  the  plasma,  and  permits  a 
spatially  resolved  calculation  of  the  balance  between  energy  absorbed 
from  the  laser  beam  and  the  energy  lost  from  the  plasma  through 
radiation,  convection,  and  thermal  conduction.  Data  have  been 
obtained  for  argon  plasmas  at  a  nominal  pressure  of  two  atmospheres 
and  four  different  bulk  flow  velocities  from  0.4  to  2.9  m/s.  The 
temperature  isotherms  for  these  four  cases  are  shown  in  Figure  1, 

Note  that  the  plasma  forms  upstream  of  the  focal  point  of  the  laser 
beam.  The  distance  between  the  temperature  maximum  and  the  focal 
point  changes  as  the  flow  velocity  changes,  and  is  extremely  important 
in  determining  the  fractional  power  absorbed  from  the  laser  beam. 

As  the  flow  velocity  increases,  this  distance  decreases,  but  the 
temperature  isotherms  become  more  elongated  and  move  into  a  region 
of  greater  laser  intensity,  thus  producing  an  overall  increase  in 
the  optical  depth  of  the  plasma.  Figure  2  shows  the  volumetric  power 
absorption  of  the  plasma  as  determined  from  the  measured  temperature 
field.  Note  that  the  maximum  volumetric  power  absorption  occurs  at 
a  position  downbeam  from  the  plasma  temperature  maximum,  an  effect 
which  is  caused  by  the  focusing  properties  of  the  converging  laser 
beam.  Detailed  examination  of  the  complex  interactions  of  the  various 
energy  absorption  and  loss  mechanisms  will  lead  to  a  more  complete 
understanding  of  the  processes  which  control  plasma  stability, 
fractional  power  absorption  and  mixing  in  the  laser  sustained  plasma. 


1984  ROCKET  RESEARCH  MEKTiHJ 
Abstract  26  P/_*  2 


*# 
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Figure  1.  Anticipated  isotherm  profiles  for  an 
argon  plasma  at  a  nominal  pressure  of  two  Atmos 
pheres  and  the  indicated  bull,  flow  velocities. 
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RESEARCH  STUDIES  IN  SUPPORT  OR  A  PULSED  LASER-HEATED 
THRUSTER  -  LASER-INDUCED  GAS  BREAKDOWN  IN  THE  VISIBLE/UV 

DAVID  I.  ROSEN  AND  GUY  M.  WEYL 

Physical  Sciences  Inc. 

Research  Park 
Andover,  MA  01810 

In  the  RP  thruster  concept  shown  in  Fig.  1,  parabolic  nozzle  walls  focus 
the  incoming  beam  to  yield  propellant  breakdown  at  the  focal  point  of  the  parabola. 
Depicted  schematically  in  Figs.  1a-1d  are  the  various  stages  in  the  thrusting 
sequence:  breakdown,  laser  absorption,  blast  wave  expansion,  and  late  time 
expansion  (when  the  flow  becomes  nearly  1-D).  Past  experimental  studies  using  CO2 
lasers  have  demonstrated  the  validity  of  the  laser  heated  thruster  concept. 

Shorter  wavelength  lasers,  however,  are  advocated  because  of  the  smaller  optics 
required  for  propagation  over  large  distances. 

We  have  carried  out  experiments  to  determine  threshold  laser  intensities 
required  for  stages  a  and  b  -  ignition  and  absorption  -  to  occur.  We  used  the 
laser  facilities  of  the  MIT  NSF  regional  laser  center:  a  Nd-Yag  laser  beam  of 
pulse  energy  E  <  0.75  J  and  Tp  =  15  nsec  in  conjunction  with  a  frequency  doubler 
and  tripler  to  generate  0.53  and  0.35  micron  radiation.  Thresholds,  Ith»  were 
measured  in  Ne,  Ar,  Xe,  and  N2  gases,  either  pure  or  with  easily  ionizable  contam¬ 
inants,  as  a  function  of  pressure  p.  Ith  was  defined  as  the  lowest  I  for  which  a 
visible  spark  was  observed  in  the  gas  together  with  significant  absorption,  a>0.05, 
of  laser  energy. 

Breakdown  is  associated  with  electron  density  growth  through  two  main 
mechanisms:  multiphoton  ionization  (MPI)  of  the  atom/molecules  in  the  intense 

laser  field  and  cascade  breakdown  due  to  electron  impact  ionization  of  the  gas  or 
electron  impact  excitation  followed  by  photoionization  of  the  excited  states.  We 
have  modeled  both  processes  using  known  MPI  cross  sections  and  using  a  Boltzman 
code  to  solve  for  the  electron  distribution  and  calculated  excitation  and  ioniz¬ 
ation  rates.  Comparison  of  theoretical  and  experimental  results  indicate  that 
breakdown  as  defined  above  is  associated  with  cascade  ionization  and  not  MPI. 

The  measured  thresholds  are  plotted  in  Fig.  2  as  a  function  of  ionization 
potential,  €j,  of  the  gas.  They  show  that  lt^  increases  with  Ej ,  with  N2  showing 
an  anomalous  behavior.  These  thresholds  are  to  be  compared  with  those  at  10.6  ym 
also  plotted  in  Fig.  2.  The  data  for  Ne  were  obtained  by  scaling  to  10.6  pm  the 
microwave  breakdown  measurements  of  Kuston  and  Fuja1  the  data  for  Xenon  was  taken 
from  Generalov  et  al.2,  N2  from  Weyl3,  Ar  from  our  previous  work,4 5  and  the  data  for 
He  from  Yablonovick . 3  Pt  scaling  was  used  to  reduce  the  results  to  p  =  3atm, 
t  =  15  nsec.  For  breakdown  dominated  by  cascade  ionization,  one  would  expect  the 
thresholds  to  scale  as  the  inverse  bremstrahlung  absorption  coefficient,  i.e.  as 
X-2.  Though  such  a  scaling  from  10.6  urn  to  0.53  ym  is  nearly  obeyed  by  Neon, 
thresholds  are  found  to  be  much  lower  for  most  gases  with  smaller  ionization  poten¬ 
tial.  This  is  because  excited  states  of  these  gases  become  rapidly  photoionized  at 
short  wavelength. 

1.  R.  Kuston  and  R.  Fuja,  J.  Appl.  Phys.  47,  498  (1976). 

2.  N.  Generalov,  V.  Zumakov,  G.  Kozlov,  V.  Masyukov  and  Y.  Raizer, 

JETP  Lett,  J_1_,  228(  1971  ). 

3.  G.  Weyl,  J.  Phys.  D,  Appl.  Phys.  Y2_,  33(  1979). 

4.  DiRosen  et  al.  PSI  TR-184  (undated). 

5.  E.  Yablonovick,  Phys.  Rev.  Lett,  32,  1101  (1974). 
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Fig.  1  Different  stages  in  pulsed  laser  thruster.  (Sequence  a-d  is  repeated 
after  initial  gas  has  expanded  out  of  nozzle  and  frest  gas  enters.) 
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HIGH  TEMPERATURE  MOLECULAR  ABSORBERS  FOR  CW  LASER  PROPULSION 

David  T.  Rosen 

David  O.  Ham  and  Lauren  H,  Cowles 
Physical  Sciences  Inc, 

Research  Park 
Andover,  MA  01810 

The  physical  processes  involved  in  heating  the  working  fluid  of  a  rocket  en¬ 
gine  with  a  high  power  CW  laser  beam  are  diagramed  in  Fig.  1.  Gas  is  injected  into 
the  absorption  zone  at  a  temperature  most,  probably  determined  by  regenerative  cooling 
requirements.  As  the  gas  flows  toward  the  throat,  it  is  heated  by  absorption  of  laser 
radiation.  With  hydrogen  as  the  primary  propellant  constituent,  the  equivalent  of 
nine  10.6  um  photons  per  molecule  must  be  absorbed  to  reach  a  stagnation  condition 
that  yields  a  specific  impulse  of  ~  1000  seconds. 

The  absorption  scheme  originally  considered  required  the  laser-induced  break¬ 
down  of  the  H2  "fuel"  followed  by  the  formation  of  a  stable  laser-supported  combustion 
(L SC)  wave.  The  principal  absorption  mechanism  in  this  case  is  inverse  electron 
bremsstrahlung,  which  requires  significant  ionization  levels  in  the  carrier  gas.  It 
has  been  suggested  that  the  introduction  of  alkali  seeds,  which  will  begin  to  thermal¬ 
ly  ionize  at  temperatures  of  ~  3000-3500  K,  would  allow  operation  at  temperatures  of 
~  6000  K,  well  below  the  temperature  required  for  ionization  in  pure  hydrogen.  It  re¬ 
mains  to  be  shown,  however,  that  stable  LSC  waves  can  be  formed  at  such  low  tempera¬ 
tures  . 

Although  the  use  of  alkali  seeds  appears  promising,  a  technique,  such  as 
laser-induced  breakdown,  is  still  required  to  heat  the  gas  to  T  -  3000-3500  K  (to  ini¬ 
tiate  alkali  ionization).  One  alternative  is  to  use  "seed"  molecules  to  absorb  the 
laser  radiation  via  vibration-rotation  band  transitions.  Such  absorbing  molecules  can 
provide  for  gas  heating  to  temperatures  of  ~  3000-3500  K,  so  that  heating  from  the  in¬ 
itially  "cold"  gas  to  stagnation  conditions  can  be  continuous  rather  than  via  laser- 
induced  breakdown.  Furthermore,  if  such  species  can  absorb  to  T  *  4500-5000  K,  then 
specific  impulses  of  1000-2000  s  can  be  achieved  without  the  need  for  ionization  (and 
thus  alkali  seeds). 

The  high  resolution  absorption  properties  of  potential  absorbers  for  this  ap¬ 
proach  are  not  well  known  at  temperatures  exceeding  1000  Important  absorption 

paths  at  elevated  temperatures  involve  hot  band  transitions  which  have  not  been  stud¬ 
ied.  Furthermore,  most  molecules  will  dissociate  in  the  temperature  rang*.'  of  inter¬ 
est,  so  that  they  no  longer  provide  effective  absorption. 

We  will  report  experimental  results  for  the  absorption  of  C02  laser  radiation 
by  the  three  most  promising  candidate  absorbers  in  the  10  un  region,  WHj.  HjO  and  cio . 
Measurements  were  made  in  the  low  intensity,  unsaturated  absorption  regime  using  sev¬ 
eral  lines  of  a  line  tunable,  CW  0 O-  laser.  Absorption  profiles  ;'••••  sun  time  we---  >b- 
tained  in  both  the  incident  and  reflected  shocks  for  ve-y  -Uluke  mixtures  of  tne  ab¬ 
sorbers  in  Ar,  shock  heated  to  temperatures  of  600-4000  K  and  pressures  of  0.5-120  at¬ 
mospheres.  Curves  drawn  through  our  experimental  data  joints  are  shown  in  Fig.  2  for 
the  P(  20)  CO  2  laser  line  at  10.59  jim.  The  absorption  coefficients  plotted  in  Fig.  .1 
are  calculated  per  arnagat  of  the  absorber  as  it  none  > i  t.n  ■  absorber  dissociates ,  Any 
dissociation  appears  in  the  measurements  -is  reduced  absorption  unless  dissociation 
fragments  also  absorb.  The  nigh  temperature  lai  j  u ,  tne  absorption  coefficient 

for  NH3  is  extremely  sensitive  to  hen tj  nq  1 *1 Rapid  sen  ting  allows  continued  ab¬ 
sorption  up  to  the  tempe  rat  e  -  .  win-/--  trie  1  r.soe  i  .1 1 ;  t -u  tat*-  e  o  -npa  rabi  o  to  th»*  -.eat¬ 
ing  rate.  These  small  siqn.t1  resu'ts  radii  -r-  *>„-•  .. •. , j ,1  laser  heating,  adequate 

absorption  to  heat  a  thruster  gas  amove  4(1D0  k  -'at :  be  attained  by  seeding  w,  ch  a/i 
H20  -  NH3  mixuture.  Sattir.it-  -  on  a"  rw-e,;  i-  t  1  ,,  ■  spe<--,  es  at  hi  oh  laser 

intensities  will  be  investigated  j  subsequent  vmr'  . 
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a)  Continuous  Absorption  Mode 


Fig.  i  Two  possible  modes  of  operation  of  a  C\-i  later-heated  thruster. 


Absorption  Coefficient  (cm-amagat) 
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Propagation  of  Laser  Supported  Plasmas 
in  Flowing  Media 

Charles  L.  Merkle  and  Thomas  M.  York 
The  Pennsylvania  State  University 
University  Park,  PA 


The  absorption  of  electromagnetic  energy  in  a  laser  thermal  thruster 
is  shown  schematically  in  Fig.  1.  An  incoming  collimated  beam  with  a 
gaussian  intensity  distribution  is  focussed  inside  the  nozzle.  The 
energy  absorption  is  determined  by  a  strongly  coupled  interaction  between 
the  nozzle  flowfield  and  the  radiation.  This  interaction  is  being  modeled 
in  a  detailed  manner  using  both  one-dimensional  and  two-dimensional 
analyses.  The  one  dimensional  analysis  is  useful  for  understanding  such 
fundamental  phenomena  as  the  propagation  mechanisms  cf  the  laser  absorp¬ 
tion  volume  and  their  sensitivity  to  the  sh.  pe  of  the  absorptivity-tem¬ 
perature  curve.  Propagation  mechanisms  being  considered  include  con¬ 
duction,  radiation,  and  non-equilibrium  effects,  individually  and  in  com¬ 
bination.  Sensitivity  to  the  absorptivity-temperature  curve  is  to  be 
found  both  by  parameterizing  in  terms  of  a  power-law  dependence  on  tem¬ 
perature  and  density,  and  by  using  specific  absorptivity-temperature 
variations  corresponding  to  selected  H2  and  IH-seedant  mixtures. 

One-dimensional  results  to  date  have  indicated  that  the  behavior 
of  the  absorption  process  inside  a  choked  nozzle  is  considerably  diffeient 
from  that  in  a  constant  area  duct.  Some  representative  results  of  c'nis 
nature  are  given  on  Fig.  2. 

Although  such  one-dimensional  analyses  are  invaluable  for  obtaining 
initial  understanding  and  establishing  trends,  they  cannot  be  expected 
to  yield  better  than  order  of  magnitude  predictions.  More  realistic 
predictions  must  take  into  account  the  many  strong  two-dimensional  effects 
noted  in  Fig.  1.  The  effects  of  parameters  such  as  beam  convergence, 
radial  intensity  gradients,  radiant  heat  loading  on  the  walls,  mixing 
between  the  hot  central  core  and  cool  outer  gases,  and  variations  in  the 
size  and  shape  of  the  absorption  volume  cannot  even  be  estimated  from  one¬ 
dimensional  results.  In  particular,  the  two-dimensional  formulation 
allows  tradeoffs  between  the  size  and  the  peak  temperatures  of  the  absorption 
volume  to  be  made.  Realistic  modeling  of  the  propagation  mechanisms  also 
becomes  possible  in  two-dimensions.  Some  preliminary  two  dimensional 
results  which  have  been  obtained  to  date  are  also  shown  on  Fig.  ?.  These 
results  are  for  beam  absorption  in  an  unchoked,  subsonic  nozzle.  The 
beam  convergence,  Mach  number  and  temperature  con lour 5;  are  presented. 


Thruster 
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A  HIGH  POWER  Nd -GLASS  LASER  INSTRUMENT 
FOR  ADVANCED  PROPULSION  AND  DIAGNOSTICS 

Charles  L.  Merkle  and  Thomas  M.  York 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


A  high  power  short  wavelength  (1.06pm)  source  is  being  developed 
by  adding  two  Nd-glass  amplifiers  to  an  Nd-Yag  oscillator  as  shown  in 
Fig.  1.  The  source  is  to  be  used  in  radiation  absorption  studies  for 
laser  propulsion  applications,  and  for  advanced  optical  diagnostics. 

A  short  wavelength  source  is  attractive  for  laser  propulsion  because 
it  is  near  a  transmission  window  in  the  atmosphere  and  because  it  allows 
smaller  collector  optics.  The  Nd-glass  instrument  also  promises  to  be 
a  more  flexible  laboratory  tool  than  typical  gas  lasers. 

The  integration  of  the  oscillator  and  the  two  high  gain  amplifiers 
must  be  accomplished  with  extreme  care  because  of  the  near-ultimate 
loading  in  the  rods  and  amplifiers.  The  combined  two-stage  amplifier 
has  a  gain  of  nearly  10  ,  and  spurious  triggering  caused  by  uncontrolled 
feedback  will  destroy  the  amplifiers.  The  maximum  power  delivered  will  be 
dependent  on  achieving  a  high  degree  of  effective  isolation  between  the 
amplifiers.  This  isolation  will  be  accomplished  by  means  of  a  Faraday 
rotation  device.  Precise  control  of  beam  nonuniformities  generated  by  the 
pump  lamps  and  by  rod  doping  must  also  be  maintained  to  prevent  amplifier 
damage.  The  concept  described  here  is  also  being  tested  at  Los  Alamos 
National  Laboratory,  but  run  time  requirements  for  the  present  application 
require  an  order  of  magnitude  extension.  A  maximum  run  time  of  10ms  is 
being  used  as  a  goal. 

The  oscillator-amplifier  system  has  two  distinct  duty  cycles.  In 
the  first,  the  instrument  operates  in  a  quasi-CW  mode  as  outlined  in 
Fig.  1.  In  this  duty  cycle,  the  oscillator  is  operated  mode-locked,  pro¬ 
ducing  a  100  nJ  pulse  in  100  ps  every  10  ns.  The  output  goal  after  am¬ 
plification  is  to  deliver  800  J  in  10  ms  in  this  mode,  an  average  power  of 
80  KW.  This  energy  delivery  represents  a  quantum  jump  over  presently 
available  systems.  As  shown  in  Fig.  2,  the  individual  pulses  in  this  wave 
train  will  be  so  closely  spaced  in  time  that  gas  molecules  have  no  time  to 
move  between  pulses  and  hence  see  an  effectively  continuous  radiation  input. 

The  second  duty  cycle  is  also  shown  on  Fig.  1.  Here,  the  oscillator 
is  Q-switched  to  produce  a  pulsed  output  at  frequencies  up  to  50  KHz. 

These  signals  are  again  amplified  to  produce  a  high  power  repetitively 
pulsed  output,  although  net  gains  in  the  pulsed  mode  are  expected  to  be 
lower  because  of  the  higher  instantaneous  intensities.  This  repetitively 
pulsed  output  may  be  of  interest  in  pulsed  laser  propulsion,  but  also 
offers  promise  for  advanced  instrumentation  concepts.  In  such  applications, 
frequency  doubling  to  0.53pm  will  probably  be  required.  The  frequency 
doubling  can  also  be  accomplished  in  the  first  duty  cycle,  and  in  addition 
wavelength  doubling  to  2.12pm  can  be  performed. 


Characteristic  Time 
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Applications 

-  Radiation  Absorp 
tion  Mechanisms 

-  Advanced  Optical 
Diagnostics 


Pulse  Duration,  s 


Time,  Sec. 

Fig.  2  Power  Output  Goals  in  Quasi-CW  Mode 
and  Use  for  Absorptivity  Measurement 
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PLASMA  INITIATION  MECHANISMS 
FOR  CW  LASER  PROPULSION  (AFOSR  #83  -  0041) 

Herman  Krier  and  Jyoti  Mazumder 
University  of  Illinois  at  Urbana-Champaign 
Department  of  Mechanical  and  Industrial  Engineering,  Urbana,  IL  61801 

Successful  utilization  of  laser  propulsion  requires  a  thorough  under¬ 
standing  of  laser-sustained  plasmas;  most  of  the  fundamental  processes  found 
in  such  plasmas  are  summarized  in  Figure  1.  A  high  energy  (10  kW  CO2)  laser 
is  used  to  provide  high  flux  densities  (>lCr  W/cnr)  to  a  gas  inside  a  pres¬ 
surized  flow-through  chamber.  Many  important  questions  remain,  particularly 
in  regard  to  plasma  initiation  and  stability,  plasma  properties,  convective 
mixing,  and  energy  losses.  It  is  our  intention  to  resolve  these  questions 
by  1)  investigating  plasma  initiation  using  targets  and  seedants,  2)  pro¬ 
ducing  detailed  temperature  and  density  measurements  in  the  plasma  core, 

3)  extending  the  temperature  data  downstream  to  define  mixing  behavior,  using 
infrared  thermography,  thermocouples,  and  laser-induced  fluorescence  (LIF), 
and  4)  estimating  global  absorption,  using  a  high-flux  calorimeter.  Figure  2 
summarizes  some  of  the  quantitative  data  we  plan  to  collect  in  the  near  future; 
each  graph  represents  a  current  analytic  estimate. 

During  the  past  year,  our  efforts  have  focused  on  the  design  and  con¬ 
struction  of  a  pressurized  flow  chamber  to  enclose  the  plasma,  and  on  instal¬ 
ling  diagnostic  equipment.  Our  first  experiments  are  planned  for  early 
April,  '84.  The  chamber  has  been  designed  for  maximum  flexibility  in  data 
acquisition,  permitting  studies  of  plasmas  under  high-pressure  (up  to  10  atm) 
steady-state  flow  conditions.  Considerable  attention  was  given  to  the  optical 
system  that  will  focus  the  laser  beam  into  desired  locations  inside  the  chamber 
in  order  to  initiate  the  plasma.  A  wide  range  of  flow  velocities  are  possible, 
and  the  five-inch  diameter  allows  complete  study  of  downstream  mixing.  The 
laser  optics  permit  precise  focusing  in  the  range  2.2  <  f  <  3.4  to  test  for 
stability  effects. 

To  measure  plasma  temperatures,  an  EG  &  G  0MA-2  (spectroscope)  system 
is  being  purchased  to  analyze  spectral  line  broadening  and  relative  line 
intensities.  Work  is  now  underway  to  correlate  spectroscopic  output  with 
temperature.  At  lower  temperatures,  we  will  be  using  an  infrared  thermo¬ 
graphy  system  to  examine  mixing  behavior.  This  system  has  recently  been 
calibrated  to  3500  K,  but  an  even  higher  upper  limit  appears  likely.  We 
also  plan  to  purchase  a  tunable  dye  laser  to  induce  LIF  in  the  gas,  allowing 
accurate  temperature  resolution  in  the  crucial  2,000  -  10,000  K  range. 

With  the  data  we  obtain  in  the  PIFC  (Plasma  Initiation  Flow  Chamber) 
it  will  be  possible  to  accurately  determine  the  processes  of  plasma  heating 
including  the  loss  mechanisms,  information  required  as  input  into  our  two- 
dimensional  steady-flow  computer  model.  Clearly  the  data  will  show  whether  the 
theoretical  absorption  coefficients  used  in  the  model  are  valid  in  the  flowing 
system  required  for  propulsion  application. 


Gas  Velocity  (cm/s) 
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AN  EXPERIMENTAL  STUDY  OF  LASER  SUPPORTED  HYDROGEN  PLASMAS 


T.  D.  McCay,  D.  M.  VanZandt  and  R.  H.  Eskridge 

Turbomachinery  and  Combustion  Devices  Branch 
Structures  and  Propulsion  Laboratory 
NASA,  George  C.  Marshall  Space  Flight  Center 
Huntsville,  Alabama 


ABSTRACT 


A  research  program  has  been  underway  since  1981  at  the  Marshall  Space  Flight 
Center  to  develop  an  understanding  of  the  physical  processes  pertinent  to  a 
laser  supported  hydrogen  plasma.  A  comprehensive  experimental  program  has  been 
undertaken  which  uses  a  30  kilowatt  continuous  wave(cw)  C02  electric  discharge 
laser  to  produce  laser  supported  plasmas  (LSP)  In  air,  argon  and  hydrogen.  A 
pulsed  laser  spark  Ignition  scheme  has  been  developed  to  permit  an  initial 
highly  ionized  plasma  to  be  produced  at  the  focus  of  the  cw  laser  and  thus  initiate 
an  LSP.  The  spark  breakdown  scheme  which  uses  a  6  Joule  per  pulse  Transverse 
Excitation  Atmospheric  (TEA)  CO2  laser  has  been  thoroughly  documented.  The 
temporal  behavior  of  the  plasma  spark  emission,  its  absorptivity,  and  the  behavior 
of  the  breakdown  threshold  with  pressure  have  all  been  examined.  Very  limited 
data  on  hydrogen  laser  breakdown  has  previously  existed  and  this  small  amount  of 
available  data  compares  favorably  with  the  current  measurements.  A  breakdown 
intensity  threshold  of  10  Gigawatts  was  measured  for  one  atmosphere  hydrogen  for 
the  100  nanosecond  pulse  length  typical  of  the  TEA  laser  used.  In  addition,  the 
mode  locking  behavior  of  the  TEA  laser  was  found  to  significantly  affect  the 
character  of  the  spark,  a  finding  which  disagrees  with  some  previous  work. 

The  experiments  which  document  the  transition  of  the  laser  spark  into  a  sustained 
plasma  supported  by  the  30  kw  continuous  wave  laser  have  only  recently  begun,  but 
data  which  documents  the  transition  and  also  the  steady  state  properties  of  the 
plasma  should  be  available  for  inclusion  in  the  paper.  This  data  would  include 
fast  frame  and  streak  photography  of  the  transition  process,  spectrographic 
surveys,  and  temperature  maps  obtained  using  a  viedo  camera  system  with  automated 
digitization  and  Abel  Inversion  data  reduction  schemes.  Details  of  the  diagnostics 
and  the  results  of  the  study,  especially  the  parametric  studies  of  plasma  absorptivity 
will  be  highlighted  in  the  article.  Schematics  of  the  simple  spark  experiments  and 
the  overall  laser  supported  plasma  experiment  are  attached. 

The  ultimate  goal  of  this  line  of  research  is  to  develop  a  sufficient  data  base 
from  which  to  design  a  thruster  which  would  employ  this  concept  for  propellant 
heating.  Based  on  concurrent  chemical  kinetics  studies  of  such  a  system  it 
appears  specific  impulse  values  as  high  as  2000  seconds  at  thrust  levels  of 
hundreds  of  pounds  could  be  obtainable  using  such  a  concept. 
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AFRPL  SOLAR-THERMAL  ROCKET  ACTIVITIES 
C.  C'.  Selph 

AF  Rocket  Propulsion  Laboratory 
Edwards  AFB  CA  93523 

ABSTRACT 

Earth  orbital  maneuvers,  particularly  orbit  raising,  provide  the  most  important 
and  extensive  propulsion  tasks  in  space  and  are  not  ideally  served  by  either  chemical 
rockets  with  their  comparatively  low  Isp  or  electric  rockets  with  their  low  thrust. 
For  several  years  the  Air  Force  Rocket  Propulsion  Laboratory  (AFRPL)  has 
investigated  solar  thermal  propulsion  as  a  potential  way  of  filling  this  gap  in 
propulsion  capability.  The  concept  makes  use  of  concentrated  sunlight  to  heat 
propellant  to  temperatures  approaching  5000  Deg  R  for  expansion  through  a 
converging-diverging  nozzle.  Typically  the  propellant  is  hydrogen,  although 
ammonia  might  be  used  where  long  term  storage  is  more  important  than  ultimate 
Isp.  No  oxidizer  is  needed,  so  molecular  weight  is  low,  and  Isp  values  of  870  seconds 
for  H2  and  500  seconds  for  ammonia  are  believed  reasonable  goals.  More  advanced 
thrusters  might  add  two  hundred  seconds  more.  The  most  characteristic  feature  of 
the  propulsion  system  is  the  pair  of  large,  paraboloidal  mirrors,  which  collect  and 
concentrate  the  sunlight  several  thousand  fold  into  the  absorber  (Fig.  I).  The 
collectors  are  45  degrees  off-axis  parabolas,  capable  of  tracking  the  sun  by  rotating 
around  an  axis  perpendicular  to  the  thrust  axis,  and  also  by  rotating  the  vehicle 
around  the  thrust  axis  as  necessary.  They  are  envisioned  as  inflatable  structures  to 
save  weight. 

A  combined  contractual  and  in-house  study  was  launched  to  provide  a 
demonstration  and  evaluation  of  the  concept  within  the  limits  of  a  ground  test 
facililty.  The  design  and  fabrication  of  a  nominal  one  pound  thrust  solar  rocket  was 
contracted  out,  and  an  in-house  facility  was  initiated  to  design  and  build  a  facility 
for  testing  the  solar  engine.  In  addition  the  analysis  of  alternative  thruster  concepts 
has  been  accomplished. 

The  ground  test  thruster  design  and  thrust  stand  facility  is  shown  in  figure  2. 
The  thruster  is  composed  of  a  cavity  and  heat  exchanger  formed  by  windings  of 
rhenium  tubing,  and  a  chamber  and  nozzle,  also  of  rhenium.  Other  than  facility 
windows  to  enable  vacuum  operation  the  thruster  may  be  classified  as  "windowless". 
The  main  research  issues  are  the  durability  of  rhenium  parts  at  these  operating 
temperatures,  and  heat  loss  control. 

Alternative  thruster  designs  for  possible  future  improvements  include  the 
rotating  bed  thruster  and  solar  plasma  thruster.  These  are  windowed  thrusters  ,and 
may  require  improvements  in  window  materials  over  what  is  available  today.  They 
share  concerns  with  the  rhenium  thruster  in  terms  of  heat  loss  control,  particularly 
the  rotating  bed,  since  there  are  bearings  ond  seals  to  protect.  Greater 
understanding  of  the  coupling  of  sunlight  via  plasmas  is  being  sought. 
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BASIC  PROCESSES  OF  PLASMA  PROPULSION 


Herbert  0.  Schrade  and  Michael  Kirschner 
Institut  fur  Raumf ahr tantr i ebe  der  Universitat  Stuttgart 


The  fundamental  understanding  of  magneto  plasmadynamic  arc  thrusters  is  still 
incomplete.  In  order  to  assess  the  applicability  of  these  thrusters  for  space 
propulsion  purposes,  and  in  order  to  improve  their  performance,  a  continuous, 
rigorous  analytical  and  experimental  investigation,  together  with  a  profound 
knowledge  of  related  and  competent  literature  is  required. 

This  research  addresses  itself  to  the  two  mechanisms  involved  which  are  of 
great  concern  in  these  thruster  devices:  1)  the  arc  electrodes  with  their 
erosion  mechanisms  and  2)  the  flow  arc  discharge  mechanism  with  its  stability 
behavior.  Recently  a  unique  analytical  approach  has  been  developed  and  so  far 
successfully  applied  to  both  problem  areas.  This  work  will  be  continued  and 
extended,  first,  in  order  to  gain  deeper  insight  into  the  very  complex  physics 
of  the  electrodes  and  the  midstream  flow  discharge  regimes,  and  second,  in 
order  to  come  up  with  a  realistic  analytical  model  which  predicts  not  only  the 
thruster  performance  limits,  but  also  allows  development  of  design  criteria 
for  improved  thruster  types . 

The  scientific  approach  towards  that  goal  is  shown  in  Fig.  1  and  is  based  on 
the  following  conception.  According  to  todays  understanding,  the  arc  cathode 
attachment  consists  of  one  or  more  plasma-channels  of  high  current  densities 
(>1017A/m2)  which  electrically  connect  one  or  more  microspots  with  the  main 
interelectrode  plasma  body.  These  channels  or  microspots  are  highly  nonsta¬ 
tionary;  they  propagate  and  disappear  and,  depending  on  the  electrode  material 
(oxidized  surface  layer,  surface  chemistry),  surface  roughness,  overall  tem¬ 
perature,  and  type  and  pressure  of  the  gas,  these  micro  discharges  cluster 
together  to  macrospots  of  different  sizes  with  different  erosion  rates.  It  is 
very  likely  that  even  a  so  called  thermionically  emitting  cathode  of  fairly 
large  overall  emitting  site  consists  of  many  highly  nonstationary  micro  dis¬ 
charges  . 

The  main  plasma  body  which  extends  from  the  cathode  to  the  anode  attachment 
area  is  confined  on  the  upstream  side  by  the  higher  pressure  cold  gas  and 
ends  down  stream  in  the  jet  plume.  The  conditions  within  this  midstream  flow 
discharge  region  can  be  considered  likewise  unstable  or  turbulent  despite  the 
overall  discharge  area  assuming  a  time  averaged  quasi  steady  configuration. 

At  first  the  overall  discharge  area  will  be  discussed  by  means  of  our  new 
stability  theory.  Secondly  a  simplified  turbulence  model  based  on  local  sta¬ 
bility  considerations  shall  be  developed  and  implemented  into  a  two-dimension¬ 
al  axial  flow  discharge  calculation.  / 

In  Fig.  2  certain  initial  results  of  an  analysis  about  a  single  micro  spot 
discharge  occuring  on  the  cathode  surface  of  a  copper  electrode  are  shown. 
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Scientific  Approach 
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a)  Crater  Radius 

b)  Crater  Pressure 

c)  Erosion  Rates 


as  function  of  time  for  3  different  initial  cathode  temperatures 


at  a  current  rate  of  2. 667*10y  A/s 


Hemispherical  Crater  Model  which 
Accounts  for  Losses  of  Evaporated 
and  Fluid  Material 


Fig.  2 :  Primary  Accomplishments. 

Initial  Arc  Cathode  Spot  Calculation  for  Copper 
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MPD  THRUSTER  PERFORMANCE 

R.  Joseph  Cassady,  R.  J.  Vondra,  and  W.  M.  Schmidt 

School  of  Aeronautics  and  Astronautics 
Purdue  University 
West  Lafayette,  Indiana  47907 

Air  Force  Rocket  Propulsion  Laboratory 
Edwards  AFB,  CA  93523 


The  thrust  delivered  by  the  Magnetoplasmadynamic  (MPD)  thruster  is  comprised 
of  two  components,  the  electromagnetic  and  the  electrothermal.  While  research  has 
been  conducted  on  the  electromagnetic  component,  the  electrothermal  contribution  to 
thrust  has  remained  largely  neglected.  Depending  on  the  operating  parameters  of  the 
thruster,  the  electrothermal  component  may  account  for  between  10  percent  and  50 
percent  of  the  total  thrust.  Therefore,  it  is  desirable  to  determine  the  pressure 
distribution  upon  the  backplate  of  the  thruster  arc  chamber  thereby  determining  the 
electrothermal  contribution  to  the  thrust. 

In  order  to  measure  the  pressure  exerted  by  the  expanding  plasma  during  the  short 
(I  millisecond)  pulse,  a  piezoelectric  pressure  transducer  was  selected  as  the  primary 
diagnostic  device.  To  overcome  the  dificulties  of  acoustic  and  electromagnetic  noise 
associated  with  such  a  high  energy  discharge,  a  nylon  cop  and  sleeve  arrangement  was 
used  to  protect  the  transducers  from  the  plasma.  The  cap  was  designed  with  a  cavity 
separating  it  from  the  transducer  diaphragm.  This  cavity  was  filled  with  glycerin  to 
act  as  a  transmission  medium  for  the  pressure  pulse.  The  front  face  of  the  cap  was 
mounted  flush  with  the  backplate  surface  and  capped  with  backplate  material.  This 
will  be  the  first  time  such  a  pressure  sensing  device  has  been  utilized  at  this  location 
within  the  thruster  arc  chamber.  These  measurements  will  provide  investigators 
modeling  the  arc  chamber  physics  with  fundamental  data  on  the  arc  chamber  pressure. 

The  effort  was  comprised  of  two  phases.  The  experimental  portion  was  performed 
at  the  AFRPL  electric  propulsion  laboratory.  Analysis  of  the  AFRPL  data  was 
performed  at  Purdue  University.  A  Benchmark  geometry  MPD  thruster  was  modified 
to  accept  the  pressure  transducers  in  July  1983.  Transducer  modifications  were 
completed  in  August,  1983.  Initial  pressure  data  were  gathered  in  September,  1983. 
Analysis  was  performed  at  Purdue  from  September  1983,  to  December  1983.  The 
conclusions  drawn  from  the  intial  data  prompted  further  measurements.  Difficulties 
with  magnetic  interference  in  the  transducer  electronics  cast  some  doubt  on  the 
precision  of  the  pressure  measurements. 

Future  plans  call  for  shielding  the  modified  transducers  with  Mu  -  metal  magnetic 
shielding  and  repeating  the  pressure  measurements.  Additionally,  B-field 
measurements  and  Kerr  cell  photographs  will  be  taken  to  supplement  the  pressure 
data. 
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PRESSURE  DECREASE 
SEEN  AT  OUTER  RADIUS 


PRESSURE  INCREASE 
SEEN  AT  INNER  RADIUS 
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MAGNETOPLASMA  DYNAMIC  THRUSTER  EROSION  RESEARCH 

David  Q.  King 
Jet  Propulsion  Laboratory 
California  Institute  of  Technology 
Pasadena,  California 

Magnetoplasmadynamic  thruster  (MPDT)  lifetime  at  sustained  multi-megawatt 
power  levels  is  unknown  but  will  be  limited  by  plasma  erosion  of  thruster 
surfaces.  Before  the  thruster  can  be  developed  for  an  orbital  propulsion 
application  the  physics  of  the  erosion  mechanisms  must  be  studied.  The 
research  proposed  herein  addresses  the  following  key  questions  that  must  be 
resolved  in  order  to  understand  erosion: 

1)  What  are  the  erosion  mechanisms  involved  and  what  are  the  quantitative 

rates  for  each  mechanism? 

2)  What  are  the  cathode  emission  mechanisms  and  heat  transfer  rates? 

3)  What  are  the  surface  temperatures? 

1)  How  does  the  work  function  of  the  cathode  change  with  time  and  what 

effect  does  this  have  on  erosion? 

The  approach  is  unique  in  that  it  aims  at  developing  an  understanding  of 
the  erosion  of  the  electrode  and  insulator  surfaces  by  conducting  experiments 
on  a  scaled-down  electrode  configuration  or  model,  which  is  under  construc¬ 
tion.  The  most  significant  physical  parameters  are  electrode  current  densi¬ 
ties  and  attachment  mode  (spot  or  diffuse),  surface  temperatures,  electron 
densities  and  temperatures  close  to  the  surface  and  total  erosion  rates  for 
each  surface  (Fig.  1). 

The  design  and  fabrication  of  a  pyrometer  has  been  completed.  Testing  of 
this  instrument  using  a  known  temperature  source  shows  that  temperature  meas¬ 
urements  in  the  range  1600-2800°k  have  a  maximum  error  of  +9*  -21  with  the 
biggest  uncertainty  being  that  of  the  spectral  emissivity. 

A  simple,  first-order  model  for  heat  transfer/electron  emission  at  the 
cathode  has  been  developed.  Despite  its  simplifying  assumptions  from  it  one 
can  make  the  following  important  conclusions:  1)  The  surface  temperature  and 
work  function  are  key  parameters, and  2)  electron  cooling  due  to  thermionic 
emission  is  significant  and  at  densities  greater  than  10^5/CB3  radiation 
cooling  can  be  neglected. 

Most  importantly,  the  model  combined  with  measurements  of  the  surface 
temperature  and  the  electron  density  close  to  the  cathode  can  be  used  to 
calculate  the  cathode-fall  voltage  and  ion  current  density.  These  parameters 
are  vital  in  determination  of  the  electron  emission  mechanisms.  Measurement 
of  the  surface  temperature  is  critical  in  that  the  evaporation  rate,  therm¬ 
ionic-emission/electron-cooling  and  sputtering  rate  all  show  an  exponential 
dependence  on  temperature;  in  addition  surface  temperature  plays  a  major  role 
in  the  overall  heat  balance  of  a  radiation-cooled  cathode. 

Similarly,  since  electron  cooling  and  thermionic  emission  show  an  expo¬ 
nential  dependence  on  work  function  experiments  will  be  conducted  using  both 
pure  tungsten  and  thoriated  tungsten  (W-Th)  which  have  work  functions  of  4.5eV 
and  2.63eV  respectively.  Under  the  same  operating  conditions  (mass  flow  and 
total  current),  measurements  of  the  temperature  and  erosion  rate  for  these 
materials  will  yield  important  results  on  cathode  erosion  (Fig.  2).  Thorium 
depletion  in  the  W-Th  cathode  will  be  systematically  studied  using  surface 
chemical  analysis  and  correlated  with  changes  in  work  function,  temperature 
and  erosion  rates. 


APPROACH:  OPERATE  SCALED-DOWN  ELECTRODE  MODEL  TO  DETERMINE 
CONDITIONS  CAUSING  MOO-SPOT-MODE  OR  DIFFUSE  MODE  CURRENT 
EMISSION,  AND  CATHODE  TEMPERATURE 
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Erosion  and  Onset  Mechanisms 

in 

Magnetoplasmadynamic  Thrusters 

John  L.  Lawless 
Carnegie-Mellon  University 
Pittsburgh,  Pennsylvania  13213 

Although  erosion  is  a  major  concern  in  MPD  thruster  design,  very  little  is 
presently  known  about  it.  To  minimize  erosion,  it  is  desirable  to  operate  in  the  diffuse 
as  opposed  to  spot  mode  of  current  conduction.  The  goal  of  this  work  is  ( 1 )  to 
determine  the  limits  of  diffuse  mode  operation  and  (2)  to  quantify  its  erosion  rates. 
Since,  at  present,  there  are  no  theories  to  answer  either  question,  this  work  is  unique. 
Also,  it  is  hoped  to  clarify  the  relationship  between  erosion  and  the  experimental 
phenomenon  of  ’onset.' 

Erosion  is  believed  to  be  due  to  evaporation  and  sputtering.  Evaporation  is  a 
serious  concern  with  tungsten  electrodes  if  the  surface  temperature  reaches  3000K.  as  it 
might  in  steady-slate  operation.  The  surface  temperature  is  determined  by  a  heat  transfer 
problem  involving  the  plasma,  the  sheath  structure,  and  the  electrode  cooling  mechanisms. 
The  importance  of  sputtering  is  determined  by  the  rate  with  which  energetic  ions  strike 
electrode  or  insulator  surfaces.  The  major  unknown  here  is  how  such  high  energy  ions 
are  created.  One  hypothesis,  due  to  Kuriki  and  Onishi,  is  that  they  are  created  m  a  high 
voltage  region  of  the  anode  sheath.  Another  possibility  is  that  they  are  created  by 
thermal  means  within  the  plasma.  The  relative  importance  of  these  two  mechanisms,  is 
determined  by  the  mean  free  path  for  energetic  ions  and  the  plasma  temperature.  The 
above  considerations  are  illustrated  in  figure  1. 

Before  erosion  can  be  quantified,  more  information  on  the  bulk  plasma  flow  is 
needed.  To  this  end,  a  flow  model  is  being  studied.  This  model  includes  the  effects  of 
pressure  gradients,  magnetic  forces,  ohmic  heating,  and  convection  in  a  quasi-one- 
dimensiona!  approximation  similar  to  that  of  King  et.  al.(  1 98 1  >.  Various  assumptions 
about  the  thermodynamic  functions  have  been  considered.  To  understand  erosion,  this 
model  is  being  extended  to  yield  more  information  on  plasma  conditions  near  the 
electrodes. 

The  primary  scientific  accomplishment  to  date  has  been  the  prediction  of  a  new 

onset  mechanism.  This  mechanism  predicts  that  the  back  emf  rises  rapidly  enough  to 

limit  the  current  to  some  critical  value.  This  is  illustrated  in  figure  2.  This 
mechanism  is  distinct  from  the  anode  sheath  failure  mechanism  analyzed  by  Baksht 
et.  al.(  i 974 >  among  others.  It  should  be  possible  for  experiments  to  distinguish  between 

these  two  mechanisms.  The  Baksht  mechanism  would  be  observed  as  a  rapid  increase  in 

the  voltage  drop  of  the  anode  sheath  near  the  exit  as  onset  is  approached.  The  proposed 
new  mechanism  would  be  observed  as  a  rapid  drop  in  current  density  in  the  middle  of 
the  thruster.  The  preliminary  evidence  exists  to  support  the  existence  of  both 

mechanisms  but  is  inconclusive  about  which  occurs  first  under  different  flow  conditions. 


Plasma  conditions 
near  the  electrodes 
may  be  affected  by 
thermal,  viscous, 
and  compositional 
boundary  layers 
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In  the  Kurikl-Onishi  mecnanism, 
high  energy  ions  which  are 
created  in  the  anode  sheath  may 
travel  to  the  cathode  and  cause 
sputtering . 


Hot  Plasma 
~20,000k 


The  Kuriki-Onishi 
mechanism  can  only 
work  if  the  mean 
free  path  for  past 
ions  is  long. 


If  the  plasma  is  hot 
enough,  ions  belonging 
to  the  high  energy  tail 
of  the  Maxwellian  distri¬ 
bution  may -cause 
sputtering. 


In  steady-state  ore.-., 
the  electrodes  may  b*. 
enough  that  evircrat  i 
an  important  cr.siou 
mechanism.  . 


FIGURE  1 

The  lifetimes  of  MPD  thrusters  are  limited  by 
electrode  erosion.  Scientific  questions  related  to 
quantifying  erosion  are  illustrated. 
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Electric  field  The  back  EMF  opposes 

drives  the  current  current  flow  in  a  thryste 
in  gii  HPD  Tiiruster 


FIGURE  2 


For  a  smooth  solution  it  is  necessary  that  uB<E. 
•Flow  modeling  predicts  that  uB  rises  faster  than  E.  This 
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Performance -Limiting  Factors 
for  MPD  Thrusters 


Research  Abstract 
for  AFOSR  Rocket  Propulsion 

Research  Meeting,  Lancaster,  CA 
March  12-15,  1984 

by 

M.  Martinez-Sanchez 
Massachusetts  Institute  of  Technology 

Both  self-field  and  applied-field  MPD  thrusters  are  subject  to  destruc¬ 
tive  instabilities  at  high  currents  and  low  flow  rates.  Several  plausible 
correlations  exist  for  the  onset  condition  for  these  instabilities,  but  they 
tend  to  be  too  sweeping  and  lack  enough  mechanistic  detail  to  allow  design 
modifications  aimed  at  avoiding  or  delaying  this  onset.  Our  research  has  as 
one  primary  goal  the  clarification  of  the  onset  mechanism(s)  by  means  of  rela¬ 
tively  simple  models  of  the  plasma  flow.  A  second  goal  is  the  extension  of 
these  models  to  a  computer  program  that  can  be  used  as  a  design  and  data  re¬ 
duction  tool,  and  as  a  means  to  study  the  tradeoffs  between  various  thruster 
conf igurations . 

In  the  course  of  this  year  we  have  considerably  extended  the  formulation 
of  a  quasi-two-dimensional  model  (Bakhst  et  al.)  which  had  had  some  success 
in  describing  onset.  This  model  was  restricted  to  fully  ionized  gas,  and 
relaxing  this  restriction  has  entailed  addition  of  equations  to  describe 
Ionization  rate,  ambipolar  diffusion  and  electron  energy  balance.  The  com¬ 
puter  program  embodying  this  extended  model  is  still  being  refined,  but  the 
major  convergence  difficulties  have  been  overcome  and  results  are  expected 
within  this  month.  These  results  will  consist  of  predicted  onset  currents 
for  given  flow  rates,  type  of  gas  and  major  geometrical  parameters,  and,  in 
addition,  predicted  distributions  of  current,  velocities,  densities,  etc., 
within  the  channel. 

We  have  also  performed  plasma  local  stability  analyses  at  various  levels 
of  detail.  Results  for  partially  ionized  regions  show  the  high  Hall  para¬ 
meter  instability,  with  little  attenuation  due  to  the  strong  magnetic  inter¬ 
action.  Extension  to  regions  near  full  ionization  is  under  way. 

We  also  plan  to  present  at  the  meeting  a  detailed  proposal  and  plan  of 
action  for  the  establishment  of  an  MPD  data  base.  We  envision  this  work  as 
a  series  of  cycles,  the  first  of  which  will  be  the  systematization  of  data 
in  the  literature  into  something  like  a  loose-leaf  data  book.  Next  will 
probably  come  a  computerized  version  of  this;  finally,  a  mechanism  for  up¬ 
dating,  maintenance,  and  access  will  be  devised. 


GDI  O 
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SCIENTIFIC  APPROACH: 

MODELING  PLASMA  DYNAMICS  AND  PHYSICS 


CONVECTED 

STRIATIONS 


APPROACH: 

(A)  EXTEND  CLASSICAL  MODEL  OF  BAKSHT  TO  INCLUDE 

-  PARTIAL  IONIZATION 

-  SELF-CONSISTENT  T 

6 

-  PRE-ONSET  REGIME 

(B)  STUDY  LOCAL  STABILITY  IN  HIGH  HALL  PARAMETER, 
PARTIALLY  IONIZED  REGIONS 

(C)  INCORPORATE  SAME  PHYSICS  INTO  NEW  3-D  CODE,  INCLUDING 

-  GEOMETRY  EFFECTS 


EXTERNAL  B-FIELD 
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EXPECTED  ACCOMPLISHMENTS 

ANALYTICAL/NUMERICAL  PREDICTIONS  OF  CONFIGURATION  EFFECTS: 


NOZZLE 

VARIATIONS 


INJECTION  PATTERN 
FUEL  CHOICE 

t 


BASIC 

CONFIGURATION 


SCALING 

I 


EXTERNAL 

B-FIELO 

FIRST  YEAR'S  EFFORT  DEVOTED  TO 
(A)  CLARIFYING  PHYSICS 


(B)  DEVELOPING  NUMERICAL  CODE 
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PULSED  INDUCTIVE  ENERGY  TRANSFER 

Peter  Mongeau  and  Henry  Kolm 
Electromagnetic  Launch  Research,  Inc. 
625  Putnam  Avenue 
Cambridge,  MA  02139 


Electromagnetic  induction  is  useful  for  converting  stored 
electrical  energy  to  the  kinetic  energy  of  a  conductive 
propellant.  Energy  stored  in  a  pulsed  power  device,  such  as  a 
capacitor,  is  discharged  into  a  high  performance  pulse  coil,  the 
rapidly  varying  magnetic  field  induces  a  current  in  a  tightly 
coupled  conductive  reaction  mass.  The  two  currents  are 
magnetically  repelled  thus  producing  thrust  (see  figure  1). 

Magnetic  induction  is  particularly  efficient  when  used  with 
solid  metals  for  reaction  mass  such  as  aluminum.  The  combination 
of  high  electrical  conductivity  and  density  of  such  materials 
enables  efficient  operation  in  the  500  -  1500  second  specific 
impulse  range  while  demonstrating  a  force  density  usually 
associated  with  chemical  propellants.  The  reaction  forces  are 
generated  with  no  direct  electrical  or  mechanical  connection 
between  the  driving  pulse  coil  and  the  exhaust  mass. 

Overall  performance  is  especially  sensitive  to  the  geometric 
matching  of  the  drive  coil  with  the  reaction  mass.  Induction 
analysis  has  shown  that  both  the  energy  transfer  efficiency  and 
thrust  are  maximized  for  thin  annular  coils.  In  this  case  the 
reaction  mass  is  a  foil  ring  placed  directly  on  top  of  the  drive 
coil.  During  the  acceleration  impulse  induced  currents  in  the 
reaction  mass  generate  substantial  ohmic  heating  that  can  lead  to 
melting  and  ultimately  vaporization  during  the  discharge.  The 
enthalpy  in  many  metals,  notably  aluminum,  is  sufficient  to 
permit  exhaust  velocities  to  beyond  10  km /sec  before  vaporization 
is  reached  (see  figure  3). 

The  operating  voltage  of  metallic  induction  reaction  engines 
limits  the  maximum  thrust,  and  consequently  exhaust  velocity, 
that  can  be  generated  while  the  reaction  mass  remains  in  the 
optimum  magnetic  coupling  region  of  the  drive  coil. 
Unfortunately,  as  higher  velocities  are  reached  the  temperature 
rise  in  the  exhaust  mass  resisitvely  damps  the  induced  current, 
thus  degrading  performance  (see  figure  4). 

Initial  experiments  have  been  conducted  in  a  glass  vacuum 
chamber  (see  figure  2).  Six  inch  diameter  drive  coils  of  various 
configurations  have  been  used  in  conjuction  with  a  90  uF,  20  kV 
capacitor  bank.  Induced  melting  has  been  demonstrated  in  0.0002 
inch  thick  aluminum  foil  at  600  m/s. 
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UNIFIED  STUDY  OF  PLASMA-SURFACE  INTERACTIONS 
FOR  SPACE  POWER  AND  PROPULSION 

P.J.  Turchi 

Washington  Research  Laboratory 
R  &  D  Associates 
30 1A  S.  West  St. 

Alexandria,  VA  22314 


Control  of  working  fluids  in  high  specific  power  and  high  specific 
impulse  devices  often  involves  plasmas  in  contact  with  solid  surfaces.  A 
three-year  basic  research  effort  is  being  pursued  in  which  plasma-surface 
interactions  important  in  space  power  and  propulsion  systems  will  be 
examined.  A  matrix  of  experimental  tests  has  been  suggested  that  will  allow 
a  broad  range  of  operating  parameters  to  be  considered  in  a  systematic 
fashion  as  indicated  in  Fig.  1.  Measurements  will  include  both  optical 
techniques  (e.g.,  spectroscopy,  combined  laser-schlieren,  -interferometry, 
-absorption)  and  probe  methods  (e.g.,  electrostatic,  magnetic,  piezoelectric 
probes),  to  determine  densities,  temperatures,  potentials,  and  flow  conditions 
near  electrically  conducting  and  insulating  surfaces.  The  use  of  a  consistent 
set  of  idealized  experimental  arrangements  will  provide  detailed  understanding 
of  processes  found  in  a  variety  of  high  specific  power  devices,  such  as  MHD 
generators  or  MPD  arcjets.  As  shown  in  Fig.  2,  low  temperature  (-0.5-I0eV) 
plasma  is  established  adjacent  to  metal  and  dielectric  surfaces,  in  the  absence 
of  flow  and  conduction,  with  flow,  and  with  both  flow  and  conduction.  The 
experiments  use  geometries  that  are  more  accessible  for  diagnostic  techniques 
and  more  tractable  for  modelling  than  arrangements  typically  developed  based 
on  performance  optimization  for  mission  needs. 

During  the  first  six  months  of  the  program,  the  basic  apparatus  was 
assembled,  including  a  0.6  x  6m.  vacuum  system  and  a  20kJ  pulse  forming 
network  (PFN) .  Preliminary  tests,  performed  in  conjunction  with  related 
work  for  AFRPL,  created  a  low  temperature  plasma  flow  (T  -  leV,  u  -  lcm/ysec) 
by  discharging  the  PFN  across  a  Teflon  surface  separating  two  brass  electrodes. 
The  PFN  current  pulse  time  was  150  psec  with  a  level  that  was  varied  from 
4-10kA.  Magnetic  probes  indicate  that  about  half  the  total  current  remains 
along  the  Teflon  surface  and  the  rest  propagates  to  the  downstream  end  of 
the  channel.  Time-integrated,  spatially-resolved  spectra  of  the  plasma  near 
the  cathode  show  several  lines  of  copper  and  zinc  that  diminish  in  relative 
intensity  by  different  degrees  with  increasing  distance  from  the  electrode. 

A  rather  uniform  carbon  and  fluorine  (intensity)  distribution  is  indicated 
across  the  channel.  Time-resolved  spectroscopy  is  planned  using  a  micro¬ 
channel-plate  image-intensif ier  system  (vs .  present  rotating-mirror  framing- 
camera  photography  that  suffers  from  inadequate  exposure  levels)  in  order  to 
monitor  the  motion  of  surface  material  into  the  freestream.  Understanding 
of  fundamental  requirements  for  charge  and  mass  transport  at  the  interfaces 
between  solid  surfaces  and  low  temperature  plasmas  should  lead  to  improved 
efficiencies  (e.g.,  lower  electrode  falls  in  MHD  and  MPD  devices),  and 
greater  performance  life  (i.e.,  less  electrode  erosion  and  insulator 
degradation) . 


SCHEMATIC  OF  PLASMA-SURFACE  INTERACTION  EXPERIMENTS 
PLACEMENT  AND  ORIENTATION  OF  THE  SURFACE  SAMPLE  IN  A  QUASI-STEADY 
PLASMAJET  EXHAUST  PROVIDES  A  RANGE  OF  PLASMA  AND  FLOW  VALUES 
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N0NC0NVI.NT1ONAL  ELECTROMAGNETIC  I’KOITLS  1  ON  CONCEPTS 


Michael  M.  Micci 

The  Pennsylvania  State  University 
University  Park,  PA 


The  electromagnetic  acceleration  of  a  nonionized  dipolar  gas  by  means  ot 
the  force  produced  by  the  cross  product  of  a  magnetic  field  and  an  alternating 
polarization  was  analyzed  to  assess  its  feasibility  as  an  advanced  propulsion 
concept.  Experimental  evidence  was  found  in  solid  dielectrics  which  verified 
both  the  existence  and  magnitude  of  the  force.  Using  benchmark  values  of  1000 
seconds  of  desired  specific  impulse  with  a  0.5  meter  long  acceleration  length, 
the  required  product  of  rotational  frequency  and  magnetic  field  strength,  c)B, 
was  calculated  for  several  types  of  propellants.  Water  has  the  highest  naturallv 
occurring  permanent  dipole  moment  to  mass  ratio  of  any  atom  or  molecule  and  gives 
a  value  of  wB  of  5-10iJ  Tesla/s  to  obtain  the  specified  performance.  However 
such  an  alternating  field  would  induce  an  electric  field  three  orders  of  magnitude 
greater  than  the  field  ionization  limit  for  water,  thus,  destroying  its  dipolar 
nature.  Use  of  light  gases  with  high  ionization  limits  such  as  helium  was 
examined.  The  acceleration  of  helium  to  exhaust  velocities  of  10^  m/s  was  found 
to  be  achievable  before  ionization  breakdown  would  occur,  however,  the  required 
product  of  wB  would  still  be  quite  high.  Electronically  excited  or  Rydberg  atoms 
or  molecules  were  studied  as  potential  propellants  due  to  their  extremely  large 
dipole  moments,  thereby  lowering  the  oB  required  for  a  given  exhaust  velocity.  It 
was  found,  however,  that  the  increased  dipole  moments  are  accompanied  by  a  greater 
decrease  in  the  ionization  field  strengths.  Since  the  maximum  acceleration 
obtainable  with  a  dipole  gas  is  limited  by  the  ionization  of  the  gas,  the  use  of 
Rydberg  enhanced  dipole  moments  actually  decreases  the  performance  of  such  a  device 
When  compared  with  other  electromagnetic  propulsion  schemes,  the  alternating  dipole 
thruster  was  found  to  give  gas  accelerations  many  orders  of  magnitude  less. 


An  investigation  was  also  inititated  to  study  the  magnetic  induction  acceler¬ 
ation  of  electrically  conducting  solid  phase  materials.  Electrical  energy  stored 
in  a  capacitor  is  discharged  through  an  inductor  which  is  in  close  proximity  to 
the  conducting  exhaust  material,  inducing  a  current  in  the  conducting  material 
which  generated  an  opposing  magnetic  field  and,  therefore,  thrust.  The  coupled 
differential  equations  describing  the  current  flows  in  the  thruster  circuit  and 
exhaust  material,  the  acceleration  of  the  exhaust  material  and  the  heating  due  to 
resistive  losses  in  the  thruster  circuit  and  exhaust  material  were  numerically 
integrated  to  determine  thruster  performance.  Previous  analyses  had  not  included 
variable  temperature  effects  and  it  was  found  that  the  heating  of  the  exhaust 
material  degraded  the  potential  performance.  A  theoretical  analysis  of  the  a>  o  1  - 
oration  process  was  undertaken  to  examine  the  coupling  of  electrical  to  kinet . c. 
energy.  A  simple  analysis  shows  that  exhaust  velocity  is  proportional  to  initial 
capacitor  voltage  and  experimental  evidence  supports  this  with  a  constant  of 
proportionality  of  0.1.  A  more  detailed  analysis  indicates  that  the  constant 
cannot  be  increased  more  than  a  factor  of  2.  Another  concern  is  the  state  of 
the  material  after  it  is  exhausted  from  the  thruster.  If  the  exhaust  material  is 
heated  to  its  melting  temperature  by  the  acceleration  process  a  vapor  pressure  i  an 
be  obtained  such  that  the  material  will  evaporate.  The  vacuum  evaporation  process 
was  examined  and  coupled  to  the  numerical  model  to  determine  the  conditions  require 
to  resistively  heat  the  exhaust  material  until  complete  evaporation  occurs. 


Finally,  the  heating  of  a  high  pressure  (>  1  atm.)  gas  by  moans  of  microwave 
radiation  was  examined.  A  literature  search  showed  that  the  vast  n;a i or i tv  of 
tile  lesenreh  in  this  area  has  been  londucted  in  the  Soviet  Union.  Botii  analytical 
and  experimental  studios  indicate  that  this  may  be  a  very  efficient  means  of 
adding  energy  to  a  working  fluid. 
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ACCELERATION  OF  DIPOLAR  GASES 


TO  ACHIEVE  10^  IN  0,5  M,  ACCELERATION  MUST  BE  10^ 


Optimal  natural  substance  in  1^0: 

-i-  =  2 •  10-i*  ^  THUS  *B  =  5-1011  KfA 

BUT  Eind|JCED  >  103  X  H2O  IONIZATION  POTENTIAL. 

Thus  H2O  ionizes  before  sufficient  acceleration, 


Optimal  nonionizing  substance  is  He: 

1^  =■  2.7'10"16'E  thus  wB  =  2.7-1012  I£fLA 

AND  ^INDUCED  I0NIZATI0N  POTENTIAL. 

Acceleration  might  be  possible  before  acceleration. 
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CLOSE-SPACED  TEMPERATURE 
KNUDSEN  FLOW 

John  B.  McVey 
Rasor  Associates,  Inc. 
253  Humboldt  Court 
Sunnyvale,  CA 


The  objective  of  this  research  is  to  investigate  and  improve  methods  for 
reducing  space-charge  in  collisionless  (Knudsen  Mode)  cesium-vapor  thermionic 
energy  converters.  The  reduction  of  the  space-charge  due  to  thermionical ly 
emitted  electrons  is  essential  for  efficient  generation  of  electric  power 
by  such  converters.  Three  approaches  will  be  investigated  both  experimentally 
and  by  computer  simulation.  These  are  illustrated  in  Fig.  1.  They  are: 

1.  Very  close  interelectrode  spacings  (10  microns  or  less). 

2.  Enhanced  thermal  ionization  of  cesium  on  the  emitter  surface.  This 
enables  the  interelectrode  gap  to  be  much  larger  (25-300  microns). 
Soviet  researchers  (Babinin,  Mustafayev,  Sitnov,  and  Ender;  1975) 
have  indicated  that  this  is  possible  by  using  structured  emitters. 

The  theoretical  basis  for  the  enhancement  is  not  well  understood, 
however.  Addition  of  barium  or  cesium  oxide  vapor  may  be  necessary. 

3.  Externally  generated  cesium  ions.  In  our  experiment,  we  plan  to 
use  radiant  emitter  heating.  This  heater  can  then  act  as  a  source 
of  positive  cesium  ions,  producing  them  by  thermal  ionization. 

Our  experimental  program  will  focus  on  the  use  of  SAVTEC  (Self-Adjusting, 
Versatile  Thermionic  Energy  Converter)  -  type  thermionic  diodes.  This  is  a 
type  of  converter  which  uses  small  diameter  electrodes,  radiant  heating  of 
the  emitter,  and  typically,  uses  thermal  expansion  of  the  emitter  electrical 
lead  to  obtain  the  interelectrode  gap.  Such  converters  have  been  shown  in 
previous  tests  to  be  capable  of  6-10  um  spacings. 

Calculational  analysis  will  be  used  to  correlate  the  experimental  data. 

An  analytical  model  will  be  developed  which  can  analyze  the  Knudsen  Mode 
discharge,  including  the  effects  of  charged  particle  trapping  and  spatially 
oscillating  interelectrode  potentials.  Existing  models  will  be  upgraded  to 
study  the  production  and  transport  of  cesium  ions  from  an  external  source 
through  a  collision-dominated  region  into  the  interelectrode  space.  The 
interaction  between  the  ion-producing  discharge  and  the  Knudsen  discharge 
in  the  interelectrode  space  will  be  handled  in  an  approximate  way,  guided 
by  experimental  results. 

The  expected  departure  of  the  experimental  volt-ampere  characteristics 
from  the  ideal,  non-space-charge  limited  diode  is  shown  in  Figs.  2a  and  2b. 

In  Fig.  2a,  the  effect  of  decreased  interelectrode  spacing  is  indicated. 

For  smaller  spacings  the  curve  is  much  closer  to  the  ideal  one.  In  Fig.  2b, 
the  expected  curve  due  to  the  unignited  mode  is  shown  for  a  spacing  of  2.0 
mils  (51  jim).  The  addition  of  ions  from  an  external  source  or  from  emitter 
structuring  could  probably  double  the  output  current  of  the  unignited  mode. 
Experimental  tests  will  determine  just  how  effective  these  ionization 
mechanisms  are. 
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heater  surface.  Cooling 


Fig.  1  Some  Mechanisms  for  Space-Charge  Reduction 
in  a  SAVTEC  -  Type  Thermionic  Converter 
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Electrode  Voltage 


Fig.  2a  Expected  Volt-Ampere  Characteristics 
for  Close-Spaced  Diodes 
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VARIABLE  FIELD-EFFECTS  RELATED  TO 
ROTATING  FLUIDIZED  BED  SPACE  POWER  SYSTEMS 

Owen  C.  Jones,  Jr. 

Rensselaer  Polytechnic  Institute 
Troy,  New  York  12181 


The  overriding  goal  of  the  program  is  to  analytically  model  and  experimental 1y 
verify  the  effects  of  variable  g-field  behavior  of  fluidized  beds  in  rotating 
systems.  The  initial  emphasis  of  the  program  is  to  evaluate  questions  which 
may  provide  limiting  or  obviating  conditions  on  the  use  of  such  systems  in 
space  power  generation. 

During  the  first  six  months  of  the  program,  questions  relating  to  the  coupled 
thermal  hydraul i c-neutroni c  stability,  thermal  stress,  and  bed  expansion 
limits  have  been  addressed.  It  has  been  determined  that  an  infinitely  long 
reactor  is  extremely  stable  with  neutronic  reactor  period  dominating  the  low 
frequency  response,  and  propellant  and  fuel  thermal  time  constants  contribut¬ 
ing  at  the  higher  frequencies.  This  evaluation  provides  a  basis  for  research 
into  the  stability  of  a  finite  reactor  when  neutronics  will  be  affected  by 
variable  neutron  leakage  due  to  the  effects  of  propellant  flow-perturbed, 
fluidized  fuel  bed  thickness.  Analysis  of  thermal  stress  in  the  blind  end- 
plate  has  resolved  the  question  regarding  the  importance  of  heat  transfer 
characteristics  from  the  core-zone  to  the  plate,  and  focused  attention  on 
mechanisms  for  active  cooling  of  this  rotating  cylindrical  plate  subjected  to 
sharp  extreme  temperature  difference  profiles. 

Finally,  a  simple  first-principals  kinematic  model  has  been  devised,  based  on 
mass  conservation  principals,  to  predict  the  initiation  and  extent  of  fluidi¬ 
zation  in  rotating  systems  as  pictured  in  Fig.  1.  The  model  developed 
confirms  in  form  a  well-known  empirical  observation  used  for  liquid-solid 
fluidization.  The  utility  of  this  model  for  gas-solid  fluidization  has  been 
demonstrated  and  the  important  scaling  criteria  have  been  identified.  The 
model  indicates  that  the  local  bed  voidage  a  (fraction  of  differential  volume 
perpendicular  to  the  local  g-field)  may  be  related  to  the  local  volume  flux  of 
gas  jg  and  the  local  terminal  velocity  v„  of  the  particles  by 

« -  (y*j1/n  m 

where  n  is  an  empirical  exponent  relating  local  relative  velocity  to  terminal 
velocity  through  the  voidage  (vr  =  anVa,).  Contrary  to  current  wisdom,  exami¬ 
nation  of  existing  data  has  shown  that  the  coefficient  in  (2)  is,  in  fact,  vm 
(Fig.  2a).  Current  controversy  regarding  the  behavior  of  n  has  been  resolved 
showing  it  to  be  independent  of  a.  Data  long  referred  to  in  the  literature 
has  been  judged  spurious  in  comparison  with  all  other  data  (the  indicated  ter¬ 
minal  velocities  of  hundreds  or  thousands  of  meters  per  second  for  50-1000-un- 
sized  particles  being  patently  false).  A  new  correlation  for  n  given  by 

n  =  5.05  [1  -  0.5  tanh  (Re„,/90ln*"  1  (?) 

is  shown  to  accurately  represent  the  existing  data  for  the  range  of  terminal 
Reynolds  numbers  of  0.06  <  Re„  <  1000  (Fig.  2b).  The  data  base  must  be 
extended  to  high  pressures  and  larger  Re*,. 


EXPANSION  EXPONENT- 
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l  itjure  2,  Comparison  of  fluidization  coefficient  (a)  and  expansion  exponent 
(b)  with  data  in  the  literature  for  gas-fluidi zed  particulate 
beds . 
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TIIEOKET I CAL  AND  KXl'K K  [MENTAL  STUDIES  OF  STABILIZED  METASTABLE  HELIUM 

Jonas  S.  Zmuldzinas 

Jet  Propulsion  Laboratory,  California  Institute  of  Technology 
Pasadena,  California  91109 

The  general  question  how  to  stabilize  metastable  excited  states  of  atoms 
and  molecules  individually  and  collectively  in  bulk  is  being  addressed  in  a 
multiphase  research  program  in  the  context  of  metastable  triplet  helium.  The 
known  long  radiative  lifetime,  weak  spin-spin  and  spin-orbit  forces,  and 
simple  atomic  structure  of  metastable  triplet  helium  make  it  the  prime  cand¬ 
idate  for  studying  this  question.  The  principal  objectives  of  the  program 
are  to  acquire  an  understanding  of  the  physics  of  bulk  metastable  triplet 
helium,  to  devise  methods  to  stabilize  it  for  long  periods  of  time,  and  to 
create  and  stabilize  a  novel  state  of  energetic  matter.  He  IV-A/M,  which  is 
a  bulk  aggregate  of  spin-polarized  metastable  triplet  helium  atoms/molecules. 
The  basic  scientific  question  and  a  three-step  approach  to  it  are  illustrated 
in  Figure  1. 

In  the  first  step,  the  surface  of  superfluid  helium  at  1-2  K  is  bombarded 
by  50  keV  electrons.  Most  of  the  excited  ^pecies  produced  in  the  liquid  decay 
fairly  rapidly,  the  survivors  being  He2  |  He2^a*5'u  or<^er  to  *SUP“ 
press  the  dominant  destruction  mechanism  He2  +  He2  ♦  4He  +  VUV,  the  He2  are 
spin-polarized  by  circularly  polarized  4650  A  CW  laser  light.  In  this  man¬ 
ner,  which  is  the  key  and  unique  feature  of  the  experiment,  concentrations 
of  He2  may  be  built  up  much  in  excess  of  thos^  previously  achieved.  Because 
of  attractive  forces  between  individual  microcrystals  of  He  1V-M 
are  expected  to  form  at  sufficiently  high  He2  t  concentrations. 

To  suggest,  guide,  and  interpret  experiments,  a  theory  of  He  IV  is  being 
developed,  first  for  He  IV-A  and  later  for  He  IV-M.  The  starting  point  is  a 
Hamiltonian  for  a  system  of  a-particles  and  electrons  interacting  via  Cou¬ 
lomb  forces.  The  weak  spin-orbit  and  spin-spin  interactions  are  initially  ne¬ 
glected,  to  be  included  in  later  stability  studies  of  He  IV.  The  formulation 
of  the  theory  employs  modern  methods  of  quantum  field  theory  at  finite  tempera¬ 
tures.  An  effective  action  for  a  system  of  electrons  and  ^He+  Ions  has  been 
obtained  to  serve  as  the  basis  for  future  numerical  calculations  on  He  IV. 
The  theory  is  expected  to  be  valid  for  a  wide  range  of  temperatures  and  pres¬ 
sures.  Preliminary  calculations  suggest  an  fee  lattice  structure  for  He  IV-A. 

Two  key  mechanisms  have  been  identified  to  deexcite  He  IV.  The  second 
step  of  the  approach  (Figure  1)  calls  for  devising  methods  to  suppress  or 
inhibit  these  mechanisms  in  order  to  realize  long-term  stability  of  He  IV. 
The  final  step  is  experimental  test  of  both  theory  and  stabilization  methods. 


The  primary  scientific  accomplishment  expected  to  be  realized  at  the  end 
of  the  program  is  illustrated  in  Figure  2.  One  way^  of  inhibiting  relativistic 
magnetic-dipole  radiative  decay  of  an  isolated  He  =  He(2JSj)  atom  is  by  the 
quantum  interference  effect.  Successful  suppression  of  this  decay  mode  would 
result  in  an  8-year  lifetime  of  He*.  Suppression  of  the  second  deexcitation 
mechanism  in  Figure  2  is  required  to  achieve  long-term  stability  of  He  IV.  To 
this  end,  a  number  of  possible  schemes  are  being  explored  at  present  by  means 
of  simple  physical  models  which  neglect  complications  such  as  lattice  vibra¬ 
tions  and  imperfections. 
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QUESTION 


•  HOW  TO  STABILIZE  He  IV -A/M  (BULK  METASTABLE  TRIPLET  He*  ATOMS/He  *  MOLECULES) 
FOR  LONG  PERIODS  OF  TIME?  * 


APPROACH 

STEP  1  -  PRODUCE  HIGH  CONCENTRATIONS  OF  SPIN -POLARIZED  HE,*  IN  SUPERFLUID  HELIUM 
(He  II)  AND  CRYSTALLIZE  He  IV -M 


50  keV  ELECTRON 
BEAM  TO 
EXCITE  He  * 


4650  A  CIRCULARLY 
POLARIZED  LIGHT  TO 
SPIN -POLARIZE  He2* 


HIGH  He2*  f  CONCENTRATIONS 


:  i  !  i  T  r 


REGULAR 
LATTICE  OF 


IV,  I 

— ) 

\  / 

—  He  II  AT  Y 

1-2  K 

o  °  °  . 

MICROCRYSTALS  OF  He  IV-M 
FORMED  BECAUSE  OF 
ATTRACTION  AMONG  He  *  f 


•  STEP  2  -  DEVELOP  THEORY  OF  He  IV  AND  DEVISE  METHODS  FOR  LONG-TERM  STABILIZATION 
OF  He  IV  AGAINST 

•  SPIN-ORBIT  AND  SPIN-SPIN  DEPOLARIZATION 

•  RELATIVISTIC  MAGNETIC -DIPOLE  DECAY  OF  He*Ate2* 

•  STEP  3  -  TEST  THEORY  AND  STABILIZATION  METHODS  IN  THE  LABORATORY 


Figure  1.  The  scientific  approach  used  in  the  theoretical  and  experimental 
studies  of  stabilized  metastable  helium. 
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SCIENTIFIC  ACCOMPLISHMENT 


•  UNDERSTANDING  ACHIEVED  HOW  TO  STABILIZE  METASTABLE  HELIUM  ATOMS/MOLECULES 
IN  ISOLATION  AND  IN  BULK  AGAINST  DOMINANT  DECAY  MODES: 


•  MAGNETIC-DIPOLE  RADIATIVE  DECAY 


ELECTRON  FLIPS  SPIN 
AND  EMITS  PHOTON 
DURING  TRANSITION  TO 
Is  LEVEL 


H- 

He  +  PHOTON 


•  SPIN  DEPOLARIZATION  DUE  TO  SPIN-ORBIT  AND  SPIN-SPIN  FORCES  WITH 

SUBSEQUENT  DOUBLE -ELECTRON  EXCHANGE  ELECTRIC-DIPOLE  RADIATIVE  DECAY 


2i_rf 

\4 —  “h- 


He*  (f)  +  He*  (4) 


He*  (f)  +  He*  (|) 


He  +He  +  2  PHOTONS 


=>  => 

RIGHT-HAND  He*  DOUBLE -ELECTRON 

ATOM  FLIPS  SPIN  EXCHANGE  +  RADIATIVE 

DECAY 


Figure  2.  The  scientific  accomplishment  expected  by  the  end  of  the 
multiphase  research  program  on  stabilized  metastable  helium. 
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ANTIPROTON  ANNIHILATION  PROPULSION 

Dr.  Robert  L.  Forward 
Forward  Unlimited 
34  Carriage  Square 
Oxnard,  California  93030  USA 

Antimatter  represents  the  most  highly  concentrated  form  of 
propulsion  energy  possible  since  the  antimatter  converts  all 
of  its  mass  to  energy  upon  annihilation  with  normal  matter. 
Antimatter  fuel  is  a  million  times  more  energetic  than 
chemical  fuel  and  a  thousand  times  more  energetic  than  nuclear 
fuel.  A  spacecraft  which  uses  antimatter  as  its  source  of 
propulsion  energy  could  "drive"  at  high  speed  anywhere  in  the 
Solar  System  and  even  travel  to  the  nearest  stars  in  a  small 
fraction  of  a  human  lifetime.  A  cruise  missile  could  circle 
the  world  on  the  "fumes"  from  its  super power ful  antimatter 
warhead . 

The  antimatter  should  be  in  the  form  of  antihydrogen  or 
antiprotons  since  antiprotons  do  not  convert  directly  into 
gamma  rays  upon  annihilation  as  do  antielectrons  (positrons.) 
As  is  shown  in  Figure  1,  when  an  antiproton  (P-)  annihilates 
with  a  regular  proton  (P  +  ),  the  annihilation  energy  is 
converted  into  elementary  particles  called  pions.  One  third 
of  the  energy  is  in  neutral  pions,  which  do  convert  into  gamma 
rays.  The  other  two-thirds  of  the  energy,  however,  is  emitted 
as  high  energy  charged  pions.  Although  shortlived,  these 
particles  travel  many  tens  of  meters,  allowing  plenty  of  time 
to  convert  their  kinetic  energy  into  thrust  by  interaction 
with  a  magnetic  field  nozzle  or  a  working  fluid.  The  pions 
decay  into  muons  and  then  electrons,  all  of  which  are  charged 
and  all  of  which  can  supply  additional  propulsion  energy.  A 
major  portion  of  our  initial  research  effort  will  be  to 
determine  the  exact  number  of  pions  and  spallation  products 
produced  as  antihydrogen  interacts  with  different  nuclei  such 
as  hydrogen,  helium,  and  heavy  elements  like  uranium  and 
mercury.  The  goal  is  to  find  a  reaction  that  distributes  the 
annihilation  energy  into  a  large  number  of  low  mass  charged 
particles  such  as  alpha  particles  (helium  nuclei). 

Antiprotons  are  already  being  generated,  captured,  cooled,  and 
stored  at  a  number  of  particle  physics  laboratories  around  the 
world,  albeit  in  small  numbers  and  with  low  efficiency. 
Figure  2  shows  one  concept  for  an  antihydrogen  factory  based 
on  extrapolations  of  the  present  research.  Each  one  of  the 
elements  has  already  been  demonstrated  on  a  small  scale.  A 
second  phase  of  our  research  is  to  study  the  physics  of  the 
present  processes  of  generation,  capture,  cooling,  and 
deceleration  of  antiprotons,  and  the  formation,  cooling,  and 
trapping  of  an t i hydrogen .  This  basic  research  will  lead  to 
techniques  for  the  efficient  generation,  long-term  storage, 
and  effective  utilization  of  milligram  quantities  of 
antiprotons  for  aerospace  propulsion. 
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NON-C ONVENT ION AL  PROPULSION 

FRANKLIN  B.  MEAD,  Jr. 

AFRPL/LKCS 
EDWARDS  AFB  CA  93523 

The  identification  and  development  of  new  propulsion  concepts  is  fundamental  to  tiie 
continued  effectiveness  of  the  U.S.  Air  Force  as  a  military  entity.  As  new  concepts 
lead  to  the  evolution  of  improved  propulsion  capability,  new  dimensions  are  added  to 
strategies  and  operations.  A  properly  conducted  advanced  concepts  program  can 
stimulate  and  vitalize  propulsion  thinking  to  the  extent  that  transitions  will  occur  at  a 
more  rapid  rate  than  otherwise  expected. 

The  non-conventional  propulsion  program  at  the  Rocket  Propulsion  Laboratory  is 
designed  to  seek-out  and  technically  assess  the  most  advanced  energy  and 
revolutionary  propulsion  concepts  which  show  promise  of  leading  to  major  advances  for 
space  propulsion;  to  identify  limitations  that  prevent  these  concepts  from  being 
feasible;  and  to  propose  research  programs  which  will  scientifically  attack  these 
limitations. 

Some  of  the  most  recent  advancements  in  non-conventional  propulsion  occurred  as  a 
result  of  a  study  of  alternate  propulsion  energy  sources  conducted  by 
Dr  Robert  Forward  during  the  calendar  year  of  1983.  The  objective  of  this  program 
was  to  look  for  breakthroughs  in  propulsion  energy  sources  that  might  lead  to  major 
advances  in  space  propulsion  in  the  next  century.  The  program  was  conducted  in  two 
phases.  Phase  I,  consisted  of  a  literature  search  and  preliminary  evaluation  that 
resulted  in  the  selection  of  a  limited  number  of  concepts  for  indepth  study  in  the 
second  phase.  Phase  II,  identified  critical  technical  research  issues  and  identified 
people  or  groups  best  qualified  to  attack  the  critical  technology.  As  a  final  task,  a 
research  program  plan  was  provided  for  each  of  the  selected  Phase  II  concepts. 

The  work  on  this  program  was  greatly  aided  by  a  number  of  previous  surveys  and  the 
cooperation  of  a  large  number  of  individuals  willing  to  discuss  their  ideas  with 
Dr  Forward.  Although  the  program  itself  is  complete,  the  Air  Force  Rocket 
Propulsion  Laboratory  continues  to  search  for  new  revolutionary  concepts  worthy  of 
its  support. 
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COUPLING  BETWEEN  VELOCITY  OSCILLATION’S 
AND  SOLID  PROPELLANT  COMBUSTION 

Robert  S.  Brown,  Anthony  M.  Blackner,  Roger  Dunlap 
United  Techonolgie s/Chemical  Systems  Division 
San  Jose,  CA  95150-0015 

Studies  are  being  conducted  to  define  and  characterize  the  basic  fluid 
mechanics  and  heat  transfer  mechanisms  controlling  the  coupling  between 
acoustic  velocities  and  solid  propellant  combustion.  The  coupling  produces 
substantial  effect  on  the  overall  stability  of  the  combustion  pressure,  yet 
substantial  deficiencies  exist  in  the  current  basic  knowledge.  Oscillatory 
heat  fluxes  and  radial  profile  of  the  mean  and  oscillatory  velocity  are 
being  measured  at  several  axial  stations  in  a  cold  flow  rocket  simulator. 

At  low  acoustic  pressure  (ie,  0.05%),  oscillatory  heat  fluxes  are 
generated  near  the  head-end  of  the  motor,  but  are  not  discernible  towards 
the  aft-end.  At  the  higher  acoustic  pressure  (0.4%)  nonlinearities  first 
appear  at  the  head-end  and  at  the  driving  frequency  further  downstream. 
Velocity  and  turbulence  profile  measurements  show  that,  at  the  lower  acoustic 
pressure,  the  acoustic  waves  exist  across  the  entire  cross-section  at 
the  head-end,  but  downstream,  are  destroyed  in  the  near  wall  region  by  the 
turbulence.  At  higher  acoustic  amplitudes,  nonlinear  flow  behavior  is 
observed  near  the  wall  towards  the  head -end  and  the  linear  oscillations  reach 
the  wall  at  the  downstream  location. 

These  data  indicate  when  and  how  the  heuristic  models  for  this  coupling 
are  deficient.  The  principal  coupling  first  occurs  at  the  head-end,  not  the 
aft-end  as  predicted  by  the  models.  Further,  the  nonlinear  behavior  also 
first  appears  at  much  lower  acoustic  pressures  than  expected  from  the  models. 


Figure  1.  VELOCITY  COUPLING  WITH  PROPELLANT  COMBUSTION 
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ANALYSIS  OF  COMBUSTION  OSCILLATIONS  IN 
HETEROGENEOUS  SYSTEMS 

MOSHE  BEN-REUVEN  AND  MARTIN  SUMMERF I  ELD 
PRINCETON  COMBUSTION  RESEARCH  LABORATORIES,  INC., 

MONMOUTH  JUNCTION,  NJ 

This  study  is  aimed  at  elucidation  of  the  physical  mechanisms 
capable  of  driving  acoustic  instability  in  solid  propellant 
rocket  motors,  particularly  of  the  type  termed  velocity-coupled 
instability.  Previous  studies  on  the  coupling  between  acoustic 
oscillations  and  the  combustion  processes  in  solid  propellant 
motors  have  demonstrated  the  complexity  of  the  overall 
phenomenon,  but  have  not  vet  defined  the  basic  mechanisms,  nor 
how  such  combustion-acoustic  interactions  occur  under  the  flow 
conditions  prevailing  inside  rocket  chambers.  The  first  and 
second  tasks  of  this  research.  Critical  Literature  Review,  and 
Order  of  Magnitude  Anlayses  of  velocitv-coupling  mechanisms,  are 
now  complete.  The  third  part  of  the  study.  Analytical  Simulation 
of  the  Interior  Flowfield  Within  a  Solid  Pronellant  Grain,  is 
near  completion.  The  subject  of  the  oresent  analysis  is  the 
cold-flow  simulation  by  Dr.  Brown  at  UTC/CSD,  in  which  nitrogen 
is  injected  through  the  porous  sidewalls  of  a  cylindrical  pipe, 
creating  an  internal  axisymmetric  flow  field.  A  comprehensive 
analytical  model  of  the  nonsteady  flow  processes  entails  a  system 
of  four  partial  differential  equations  for  continuity,  radial  and 
axial  momentum,  and  thermal  enthalpy.  The  flowfield  is  considered 
compressible  and  viscous,  so  that  all  of  the  associated 
dissipative  terms  are  included.  A  particular  focal  point  of  the 
analysis  has  been  the  thin  viscous  sublayer  adjacent  to  the 
porous  surface,  where  visco-acoustic  interactions  prevail.  Order 
of  magnitude  analysis  (based  on  previous  theoretical  works  by 
Lighthill  and  Illingworth,  precluding  injection  at  the  wall),  as 
well  as  earlier  experimental  data  (regarding  great  enhancement  of 
wall  heat  transfer  in  pipes  at  the  presence  of  a  sufficiently 
intense  acoustic  field),  indicate  that  this  type  of  nonlinear 
coupling  can  influence  the  nonsteady  surface  heat  feedback 
appreciably,  and  hence  is  a  highly  plausible  velocity-coupled 
instability  mechanism.  The  analytical  approach  to  visco-acoustic 
coupling  is  summarized  in  Fiq.  1.  A  singular  perturbation 
analysis  was  carried  out  for  the  viscous  sublayer;  a  simple 
algebraic  expression  for  the  first  order  axial  pressure  drop 
within  the  port  was  obtained,  correlating  the  available  UTC/CSD 
experimental  (steady  state)  data  excellently,  demonstrating  the 
inherent  Reynolds  No.  effects,  as  shown  in  Fig.  2.  A  numerical 
finite-difference  aloorithm  has  been  developed,  using  a  modified 
(un-split)  MacCormack  integration  scheme;  the  code  utilizes 
nonuniform  radial  stepsize,  with  verv  fine  mesh  near  the  Dorous 
surface.  Results  of  converged  steadv-state  cases  are  currentlv 
generated.  A  second  numerical  algorithm,  using  a  split 
MacCormack  scheme  (to  accomodate  for  the  great  differences  in 
radial  steosize,  while  satisfying  the  necessary  CFL  stability 
condition),  is  near  completion,  and  will  be  used  to  simulate  the 
nonsteady  cases,  where  the  exit  flow  is  perturbed  per iod ical ly. 
Results  will  be  available  in  March  1984. 
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APPBQACIk 

ANALYSIS  OF  VELOCITY-COUPLED  INSTABILITY  MECHANISMS. 
ONE  DOMINANT  MECHANISM  IS  VISCO-ACOUSTIC  COUPLING: 


FREQUENCY-DEPENDENT  SURFACE 
HEAT  TRANSFER  (FIRST  ORDER) t 
WITH  PHASE  RELATIVE  TO 
ATTENDANT  PERTURBATION. 


iToKIfS 


DC-COMPONENT:  ACOUSTIC  STREAMING, 
WHEN  U,  V  NOT  OUT  OF  PHASE ,  NET 
X-MOMENTUM  TRANSFER  OCCURS. 
(SECOND  ORDER), 
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ORDER  OF  MAGNITUDE  ANALYSIS  INDICATES  EXPECTED  EFFECTS  OF 
STROUHAL  NO.  AND  MEAN  FLOW  MACH  MO.  ON  SURFACE  HEAT  TRANSFER. 
THESE  INCLUDE  EFFECTS  OF  FREQUENCY ,  X/D  AND  (IMPLICITLY) 
PRESSURE. 


FIG.  1 
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SOLID  PROPELLANT  GAS  PHASE  FLAME  DRIVING  OF 
AXIAL  COMrUISTION  INSTABILITIES 

Ben  T.  Zinn  and  Brady  R.  Daniel 
School  of  Aerospace  Engineering 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332 

This  recently  initiated  research  program  consists  of  theoretical  and  experimental 
investigations  of  the  interaction  between  solid  propellant  gas  phase  flames  (to  be  denoted 
by  SPGPF)  and  the  flow  oscillations  in  rocket  motors  experiencing  axial  instabilities.  The 
objectives  of  this  program  are  (L)  elucidate  the  processes/mechanisms  which  contribute  to 
the  damping  or  driving  of  the  combustor  oscillations  by  SPGPF  and  (2)  determine  whether 
state  of  the  art  models  are  capable  of  describing  the  interaction  of  SPGPF  with  unstable 
combustor  oscillations. 

Since  detailed  diagnostics  of  actual  SPGPF  are  presently  impossible  because  of 
their  extremely  small  dimensions,  this  research  program  investigates  the  characteristics 
of  flat,  premixed  flames  which  possess  some  of  the  important  features  of  SPGPF  and  are 
amenable  to  experimental  probing.  The  proposed  experimental  setup,  shown  in  Fig.  1,  is 
designed  to  simulate  the  conditions  in  an  unstable  rocket  combustor.  It  consists  of  a  tube 
with  a  rectangular  cross  section.  The  investigated  flame  is  positioned  along  one  of  the  side 
walls  to  simulate  the  burning  of  a  solid  propellant.  A  second  flat  flame,  located  at  one 
end,  provides  a  steady  flow  of  hot  combustion  products  past  the  investigated  flame.  An 
acoustic  driver  at  the  opposite  end  excites  a  standing  wave  of  desired  properties  in  the 
tube.  The  expected  structure  of  the  investigated  flame  is  shown  on  the  bottom  of  the 
figure.  Proposed  diagnostics  will  consist  of  measurements  of  acoustic  pressures,  the 
steady  and  unsteady  distributions  of  u  and  v,  the  steady  state  temperature  distribution 
T(y)  and  the  injection  plate  admittance  Y  .  The  measured  T(y),  Y  ,  p'  ,  V(y  =  0)  will  be 
input  into  the  developed  theoretical  model  which  will  be  used  to  predict  v'  (y).  The 
validity  of  the  model  will  be  checked  by  comparing  the  predicted  and  measured 
distributions  of  v'  (y).  These  studies  will  be  conducted  at  various  wave  frequencies  and 
amplitudes,  different  flame  mixture  compositions,  different  magnitudes  of  the  mean  flow 
u  and  different  locations  of  the  flame  along  the  standing  longitudinal  wave. 

To  date,  using  the  boundary  layer  conservation  equations,  an  approximate  model  of 
the  investigated  flame  and  an  exact  model  of  the  oscillatory  cold  tlow  near  the  injection 
plate  in  the  absence  of  combustion  have  been  developed.  Figure  2  presents  predictions  of 
the  cold  flow  model  for  the  distributions  of  the  real  parts  of  the  normal  velocities  at 
different  axial  locations  within  that  half  of  the  combustor  where  p'  is  negative.  The 
computations  assumed  that  the  real  part  of  the  injection  plate  admittance  Y  is  negative 
(i.e.,  damping)  and  they  show  that  due  to  viscous  effects  v^fy)  (and  the  damping)  increase 
with  distance  from  the  injection  plate.  Furthermore,  these  results  show  that  the 
magnitude  of  the  damping  is  strongly  dependent  upon  axial  location  within  the  combustor. 
Additional  studies  to  date  showed  that  the  magnitude  of  the  damping  is  a  strong  function 
of  Y  and  a  weak  function  of  the  mean  flow  Mach  number  M.  Future  studies  will 
invesrfgate  whether  the  presence  of  combustion  within  the  acoustic  boundary  layer  will 
result  in  net  driving  by  the  near-wall  region  which  includes  the  viscous  layer  and  the 
flame.  Also,  these  studies  will  concentrate  on  identifying  the  SPGPF  processes  which 
exert  the  greatest  influence  upon  the  flame  driving.  The  validity  of  these  predictions  will 
be  checked  by  comparing  them  with  experimental  data. 


COMBUSTIBLE 

MIXTURE 


Proposed  Experimental  Setup  (above)  and  the  ^.nricipat^o 


(Dimensionless) 
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Figure  2.  Dependence  of  Normalized  Real  Part  of  v'Cy)  Upon  Normalized 
Distance  from  the  Injector  Plate. 
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RELATION  BETWEEN  COMPOSITE  PROPELLANT 
HETEROGENEITY  AND  NONSTEADY  BURNING 


Leon  Strand  and  Norman  S.  Cohen 
Jet  Propulsion  Laboratory/California  Institute  of  Technology 
Pasadena,  California  91109 


Rocket  motor  development  experience  and  laboratory  combustion  data  have 
indicated  that  ammonium  perchlorate  particle  size  distribution  has  a  major  ef¬ 
fect  on  combustion  instability.  The  effect  cannot  be  explained  by  classical 
theories  which  relate  the  combustion  response  of  a  solid  propellant  to  its 
ballistics  properties.  This  research  addresses  the  question  of  why  this  is 
so,  and  has  hypothesized  that  the  answer  is  contained  within  the  heterogeneity 
of  the  composite  propellant  structure.  The  answer  is  of  long-range  signifi¬ 
cance  because  combustion  instability  is  a  pervasive  problem  and  propellant 
formulation  adjustment  is  the  most  convenient  remedy.  Further,  nonsteady  com¬ 
bustion  behavior  and  composite  propellant  fine  structure  are  of  interest  in 
many  other  contexts. 

The  approach  combines  analytical  modeling  and  experiments  to  characterize 
the  heterogeneity  of  composite  propellants  and  its  relationship  to  the  combus¬ 
tion  response.  Particularly  unique  features  include  an  analysis  of  the  com¬ 
bustion  response  to  fluctuations  in  propellant  composition,  measurement  of  the 
fluctuations  in  composition  as  inherently  contained  within  unburned  propel¬ 
lant,  and  measurement  of  dynamic  burning  at  constant  pressure  which  purported¬ 
ly  manifests  such  fluctuations.  Additional  measurements  to  complete  the  ap¬ 
proach,  showis  in  Fig.  1,  are  of  size  distribution  and  the  combustion  response 
under  oscillatory  pressure  conditions.  It  is  desired  to  see  if  characteri st i c 
frequencies  exist  in  propellants,  if  they  appear  in  the  course  of  dynamic 
burning  at  constant  pressure,  and  if  they  influence  the  combustion  response  to 
pressure  perturbations. 

Theoretical  accomplishments  have  appeared  in  AIAA  papers  83-0476  and  83- 
1Z70,  and  a  nonlinear  analysis  is  in  progress.  The  major  conclusion  is  that 
the  combustion  response  to  compositional  fluctuations  can  be  relatively  large. 
Thus  if  compositional  fluctuations  exist  at  the  same  frequency  as  pressure 
perturbations,  the  combustion  response  to  the  pressure  perturbations  will  be 
augmented  significantly  (Fig.  2).  The  planned  series  of  experiments  have  been 
completed  and  the  data  are  being  processed  and  studied  at  this  writing. 


\wu. 
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MECHANISM  DRIVES  COMBUSTION  INSTABILITY 

i 


MEASURE  BY  MICROWAVE 
DETECTION  OF  BURNING 
UNDER  OSCIUATORY  v 
PRESSURE  \  \ 


MEASURE  BY  MICROWAVE 
DETECTION  OF  BURNING 
AT  CONSTANT  PRESSURE 
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ANALYSIS  IN  AIAA  PAPER 
83-0476 


MEASURE  BY  INTEGRATED 
CHLORINE  SCANS  OF 
PROPELLANT  SURFACES 


BURNING  YIELDS  COMPOSITIONAL  (a  I 
FLUCTUATIONS  AT  CHARACTERISTIC  FREQUENCIES 
WHICH  FEED  BACK  TO  AFFECT  BURNING 


COMPOSITE  PROPELLANT 
FINE  STRUCTURE 
(HETEROGENEITY) 


SIZE  DISTRIBUTION  LOADED 
INTO  PROPELLANTS 


Sature  of  the  approach 


Figure  1  . 
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OSCILLATION  FREQUENCY 
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INTERACTION  OF  MULTI-DIMENSIONAL  MEAN  AND 
ACOUSTIC  FIELDS  IN  SOLID  ROCKET  COMBUSTION  CHAMBERS 


Joseph  D.  Baum  and  Jay  N.  Levine 
AFRPL/DYC 
Edwards  AFB,  CA  93523 


The  objective  of  this  research  program  is  to  develop  understanding  of  the  physical 
mechanisms  by  which  energy  is  exchanged  between  the  mean  and  acoustic  flow  fields  in 
a  combustion  chamber,  specifically  the  mechanisms  by  which  mean  flow  energy  is 
transformed  into  acoustic  energy  through  vortex  shedding  and  the  mechanisms  by  which 
energy  is  transferred  between  the  acoustic  field  and  the  mean  flow  of  combustion 
products  entering  the  chamber. 

A  new  numerical  code  that  couples  the  solutions  for  the  flow  and  combustion  in 
the  chamber  is  currently  under  development.  The  flow  in  the  rocket  combustion 
chamber  is  described  via  the  solution  of  the  two-dimensional,  time  dependent,  turbulent 
compressible  Navier  Stokes  equations.  The  transient  burn  rate  of  the  solid  propellant  is 
modeled  using  the  modified  transient  one-dimensional  heat  conduction  equation  utilized 
in  the  previous  one-dimensional  nonlinear  combustion  instability  analysis.  A  algebraic 
grid  generation  procedu-e  is  incorporated  in  the  program  to  provide  both  fine 
distribution  of  mesh  points  in  the  boundaries  and  in  the  free  shear  layer.  The  program 
currently  includes  the  explicit  unsplit  MacCormack  finite  difference  integration 
scheme.  Using  an  explicit  scheme  results  in  a  very  small  integration  time  step.  Several 
options  are  currently  under  examination,  including:  (a)  zonal  approach,  where  an 
explicit  scheme  is  utilized  near  the  center  line  and  implicit  schemes  at  the  nozzle  and 
near  the  boundaries;  and  (b)  utilizing  an  impficit  scheme  everywhere. 

The  developed  code  will  be  utilized  to  study  the  interactions  between  vortices 
shed  from  a  sharp  downstream  facing  edge  and  the  acoustic  field  inside  the  chamber. 
Specifically,  the  objective  is  to  understand  the  mechanisms  by  which  vortex  energy  is 
transformed  into  acoustic  energy  and  the  conditions  under  which  an  acoustic  field  can 
influence  the  existence,  frequency,  and  strength  of  vortex  shedding.  The  program  will 
examine  the  coupling  between  the  mean  and  acoustic  fields  in  a  chamber  in  which 
combustion  products  are  provided  by  propellant  burning  at  the  head-end,  while  the 
chamber  walls  are  non-interacting.  Several  chamber  geometries  will  be  examined.  The 
results  from  these  numerical  experiments  will  be  examined  based  upon  the  effect  of 
mean  flow,  cavity  length  to  diameter  ratio,  and  cavity  harmonic  content  upon  acoustic 
energy  transfer  between  vortices  and  acoustic  pressure  oscillations;  reflection  of  waves 
from  the  bottom  of  the  cavity;  reflection  and  amplification  or  attenuation  of  waves 
from  the  trailing  wall;  significance  of  trailing  wall  height;  spatial  and  temporal 
propagation  of  large  scale  vortices;  interaction  of  shear  layer  with  reflected  acoustic 
waves  and  large  scale  vortices;  and  the  effect  of  blowing  from  the  surface  of  hot 
combustion  products  upon  the  large-scale  eddies  and  shear  layer  stability. 

The  developed  code  will  be  utilized  to  study  the  interaction  between  the  mean 
flow  of  combustion  products  leoving  the  burninq  surface  ond  the  ocoustic  oscillations  in 
the  combustion  chamber.  Conservation  of  mass,  momentum  and  energy  will  be 
calculated  at  several  grid  points  close  to  the  burning  propellant  surface,  in  order  to 
determine  the  momentum  and  energy  flux  rates  in  and  out  of  the  cell.  Following 
several  cells  across  the  chamber  area  cross-section  (including  cells  in  the  acoustic 
boundary  layer)  will  enable  understanding  and  determination  of  the  acoustic  energy  loss 
due  to  flow  turning. 


» 


** 


•  How  is  Vortex  Energy  Transformed  into  Acoustic  Energy 

•  Influence  of  Acoustic  Field  or.  Sheer  Layer  Instability 

•  Energy  Exchange  Between  rustic  f  ield  and  Transverse  Mas 
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FINITE  ELEMENT  ANALYSIS  FOR  COMBUSTION  INSTABILITY 
OF  SOLID  PROPELLANT  ROCKET  MOTORS 

T.  J.  Chung 

Department  of  Mechanical  Engineering 
The  University  of  Alabama  in  Huntsville 
Huntsville,  AL  35899 

The  purpose  of  this  research  is  to  provide  computational  schemes  via  finite 
elements  to  analyze  combustion  instabilities  which  may  occur  in  solid  propellant 
rocket  motors.  Specifically,  rigorous  flow  field  calculations  and  determination 
of  the  effect  of  correct  flow  field  on  combustion  instability  constitute  the  most 
important  part  of  the  current  research. 

Achievements  of  the  past  year  consist  of  (1)  derivations  of  stability  inte¬ 
grals  which  allow  acoustic  oscillations  to  interact  with  vortical  oscillations 
and  (2)  development  of  computer  programs  to  implement  both  the  flow  field  equa¬ 
tions  and  stability  integrals.  The  numerical  approach  used  is  based  on  two-  and 
three-dimensional  isoparametric  finite  elements  (Fig.  1)  which  lend  themselves  to 
efficient  numerical  integrations  by  means  of  Gaussian  quadrature. 

For  irregular  geometries  containing  inhibitors  or  transitions  of  cross  sec¬ 
tions  (Fig.  1),  it  has  been  shown  that  there  are  significant  interactions  between 
the  acoustic  and  hydrodynamic  oscillations.  It  has  also  been  made  evident  that 
such  a  study  can  not  be  performed  unless  accurate  flow  field  is  determined,  inclu¬ 
ding  both  mean  flow  and  oscillatory  flow  fields  (Fig. 2, 3).  Depending  on  a  parti¬ 
cular  vortical  frequency  mode,  the  fluctuating  vortical  components  of  the  velocity 
are  expected  to  vary  widely.  The  most  significant  aspect  of  the  present  study  is 
the  determination  of  the  stability  boundaries.  Fig.  4  shows  variations  of  hydro- 
dynamic  instability  growth  constants  (a^)  vs  vortical  oscillation  frequencies 
interacting  with  selected  acoustic  frequencies.  Although  these  graphs  indicate 
frequency  ranges  of  instability,  it  is  meaningful  to  show  stability  boundaries 
(Fig.  5).  The  general  trend  is  that  with  an  increase  of  Reynolds  numbers,  the 
stability  boundaries  shift  to  the  lower  frequencies,  a  trend  similar  to  a  flat- 
plate  boundary  layer  flow.  Note  that  multiple  "islands"  develop  due  to  inhibitors. 
Critical  Reynolds  number  appears  to  be  around  Re=700.  In  order  to  compare  the 
acoustically  interacting  vortical  instability  with  "acoustics  free"  hydrodynamic 
instability  (Fig.  6),  the  corresponding  growth  constants  are  shown  in  Fig.  7. 

It  is  seen  that  the  lower  frequency  instabilities  are  absent  in  the  lower  ranges 
of  Reynolds  numbers  in  hydrodynamic  oscillations  alone  not  interacting  with  acous¬ 
tics,  a  possibility  that  less  stable  conditions  arise  when  vortical  oscillations 
do  interact  with  acoustic  fluctuations.  However,  this  trend  may  not  be  true  in 
the  case  of  other  combinations  of  acoustic  and  vortical  frequencies. 

To  determine  the  effect  of  flow  fiejd  due  to  tapered  cross  sections  with 
various  transition  angles  on  stability,  we  have  analyzed  an  axi symmetric,  motor 
with  transition  angles  of  30°,  45°,  60°,  and  90°.  Although  the  results  are 
not  shown  in  figures  it  was  concluded  that  the  larger  the  transition  angle,  tin- 
less  stable  is  the  motor  as  a  result  of  unfavorable  interactions  of  vortical 
fluctuations  with  acoustic  fluctuations.  The  unstable  hvdrodvnamic  instability 
growth  constants  are  positive  but  they  are  brought  to  negative  side  as  summed 
to  the  acoustic  instability  constants. 
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Fig.  1  Axisymmetric  motor  with  inhibitors 
and  Finite  Element  modeling 
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ELECTROLYTIC  PREPARATION  OF  NOVEL  AZIDODINITRO  COMPOUNDS 
Milton  B.  Frankel  and  James  F.  Weber 

Rockwell  International  Co  .  oration 
Rocketdyne  Division 
Canoga  Park,  CA  91304 

Future  higher  performance  propellants  will  require  new  classes  of  more  ener  - 
getic  ingredients.  These  will  include  new  oxidizers,  plasticizers,  and 
binders.  Polynitro-  and  azidoalkyl  compounds  have  either  been  or  are 
currently  being  evaluated  for  this  purpose.  Highly  substituted  geminal  dinitro 
compounds  [e.g.,  R-C(N02)2X,  where  X=CN,  CH-j,f]  are  of  particular  interest, 
and  substitution  of  an  azido  group  (X=Nj)  will  further  enhance  the  energetics 
of  such  compounds.  This  should  result  from  the  fact  that  this  substituent  can 
contribute  ^80  Kcal/group  of  energy  to  the  compound  without  a  decrease  in  the 
0/C  ratio. 

Only  two  compounds  have  been  reported  which  contain  the  -CfNC^^^  group.  The 
approach  under  investigation  is  the  synthesis  of  1, 1 , 1-azidodini tromethyl  com¬ 
pounds,  containing  various  functional  groups,  by  electrochemical  techniques. 

The  program  is  divided  into  four  interrelated  phases:  1)  synthesis,  by  classi¬ 
cal  means,  of  the  requisite  reactants;  2)  electroanalytical  studies  on  the 
reactants;  3)  exploratory  electrochemical  syntheses  of  the  model  compounds; 
and  4)  synthesis  of  the  target  -C(N02)2^3  substituted  propellant  ingredients. 
The  second  and  third  phases  tend  to  direct  the  other  two.  The  electroanalyti¬ 
cal  studies  delineate  suitable  electrode  materials,  and  optimum  reaction 
(temperature,  solvent,  and  pH)  conditions  and  also  provide  mechanistic  infor¬ 
mation.  The  model  reactions  direct  the  selection  of  other  possible  reactants 
and  means  by  which  the  desired  products  can  be  obtained.  Model  reactions  also 
yield  information  about  optimum  concentrations  of  substrates,  order  of 
addition,  and  current  density, 

Electroanalytical  studies  on  the  anodic  oxidation  of  dinitrocarbanions  and 
the  azide  ion  in  water  and  acetonitrile  have  shown  both  to  be  electrochemi- 
cally  irreversible  and  relatively  insensitive  to  temperature  (10  to  50C)  and 
pll  (6  to  11).  Voltammetry  also  indicates  a  simultaneous  oxidation  of  the 
dinitro  carbanions  and  the  azide  ions.  This  means  that  both  radicals  formed 
are  present  in  the  reacti.e  layer  near  the  electrode  surface  and  that  the 
product  is  formed  by  a  radical-radical  coupling  mechanism  rather  than  a  step¬ 
wise  oxidation  -  addition-oxidation  (ECE)  mechanism.  Although  graphite 
appears  to  be  the  best  electrode  material  in  an  electrochemical  sense,  its 
undesirable  physical  and  mechanical  properties  have  led  to  a  preferential  use 
of  platinum. 

Electrochemical  azidization  of  the  model  compound  1, 1-dini troethane  has  given 
1 , 1 , 1-azidodinitroethane  in  40%  yield.  This  yield  was  obtained  by  conducting 
the  synthesis  in  a  divided  cell  with  substrate  concentration  greater  than  one 
molar,  and  a  current  density  maintained  at  100  mA/cm^.  The  results  obtained 
from  the  model  substrate  1, 1-dinitroethane  are  currently  being  applied  to 
other  dinitro  compounds  containing  ether,  ester,  and  acetal  functioning  groups. 
Since  the  target  compounds  contain  such  groups,  it  is  of  importance  to  under¬ 
stand  the  effects  of  these  groups  on  the  desired  azidization  reaction. 


ELECTROCHEMICAL  SYNTHESIS  OF  ENERGETIC  RC(N09)9N?  COMPOUNDS 
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SYNTHESIS  AND  CHEMISTRY  OF  ENERGETIC  METALLOTET RAAZADI ENE S 

William  C.  Trogler 
Department  of  Chemistry,  D-006 
University  of  California,  San  Diego 
La  Jolla,  California  92093 

The  goals  of  our  research  are:  1)  explore  the  products  obtained  from 
the  reaction  between  organic  azides  and  organotransition  metal  compounds; 

2)  develop  new  approaches  to  the  synthesis  of  unsaturated  (and  hence  ener¬ 
getic)  metal-nitrogen  ring  compounds;  3)  study  physical,  chemical  and  photo¬ 
chemical  properties  of  the  new  compounds.  To  date  there  are  few  known 
organometallic  compounds  that  contain  unsaturated  nitrogen  ligands,  and 
virtually  nothing  is  known  of  their  physical  and  chemical  properties.  Three 
research  problems  have  been  investigated  in  the  last  year.  Our  synthesis  of 
a  new  ligand  useful  for  the  preparation  of  unsaturated  metal-nitrogen  ring 
compounds  is  nearly  complete  (Figure  1).  Because  the  organic  synthesis  of 
unsaturated  nitrogen  compounds  has  been  largely  neglected  since  the  late 
1800's  we  have  had  to  develop  conditions  for  optimizing  each  step  in  Figure 
1.  The  preparation  of  a  series  of  metal lotetraazadienes  (Figure  2)  and 
their  chemical  or  electrochemical  reduction  to  the  radical  anions  has  been 
completed^-.  All  of  the  radical  species  exhibit  an  isotropic  epr  spectrum 
(g  =  2.16-2.21)  in  solution  characteristic  of  cobalt  centered  radicals 
(ACo  =  50-58  Gauss).  The  acceptor  orbital  is  ~50%  cobalt  3d  and  the  remain¬ 
der  consists  of  nitrogen-based  pir  orbitals  of  the  metalloycylic  ring.  The 
exploratory  synthesis  of  early  transition  metal  complexes  has  also  been  ini¬ 
tiated.  We  have  preliminary  results  that  show  Mo(CO)^,  V(C0)k,  (n-M^Hj) 
V(C0)  ^  and  (r^-C^^^V  all  react  with  xylyl  azide  or  methyl  azide  to  yield 
intensely  colored  products.  Characterization  of  these  reactions  will  form 
much  of  the  basis  for  next  year's  research. 
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Figure  1:  Synthesis  of  a  New  Energetic  Ligand 


1,4-d 11 ithio-1, 4-diphenyl  tetrazene 
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Reduction  potentials  vs.  NHE  =  -1. 53v(CH^) ,  -l.Olv(C^H^) , 


-0.71v(C,Fr) ,  -1. 31v(2 ,6-(CH_) ,,C,H_)  . 
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HMX  COMBUSTION  MODIFICATIONS 

Joseph  E.  Flanagan,  Milton  B.  Frankel,  Dean  0.  Woolery 

Rockwell  International  Corporation 
Rocketdyne  Division 
Canoga  Park,  CA  91304 

An  understanding  of  the  decomposition  mechanisms  at  heating  rates  and  temp¬ 
eratures  similar  to  those  experienced  in  actual  propellants  is  required  so 
that  the  flame  chemistry  can  be  properly  altered  via:  1)  modification  of 
the  ingredient's  molecular  structure;  or  2)  incorporation  of  additives  which 
will  yield  decomposition  products  that  interact  with  the  cyclic  nitramines 
decomposition  products. 

Experimental  evidence  from  pyrolysis  decomposition  studies  suggests  that 
structural  symmetry  of  HMX  derivatives  must  be  maintained  to  arrive  at  com¬ 
parable  rates  of  C-N  and  N-NO^  bond  cleavage  resulting  in  increased  energy 
release  and  generation  of  reactive  gas  phase  species.  In  addition,  it  has 
been  found  that  structural  symmetry  is  a  key  to  obtaining  increased  melting 
or  liquefaction  temperatures. 

The  effort  during  the  last  year  has  focused  upon  expanding  the  data  base  for 
the  systems  which  generate  amidogen  radicals  (N^ • )  via  NIL  and  the  inter¬ 
action  of  NH^’  with  nitric  oxide.  TAGZT  was  the  primary  NH^  source.  The 
burning  rate  acceleration  via  TAGZT  was  correlated  in  terms  of  the  equation 

r  =  A  exp  (-E/RT,.) 

Selected  experiments  were  carried  out  in  which  TVOPA  was  utilized  to  generate 
dif luoroamino  radicals  (NF  • )  to  compare  the  relative  augmentation  rates. 
Propellants  containing  TVOPA  and  TAGZT  had  higher  augmentation  rates  than 
-0N0„  plasticized  systems,  but  the  non-TAGZT  containing  systems  burned  at  an 
equal  rate  for  comparable  flame  temperature  (NG  vs.  TVOPA).  The  NF0’  radicals 
appear  to  have  a  "catalytic"  effect  on  the  overall  controlling  mechanism. 

The  structural  modification  of  HMX  to  replace  one  -N0?  group  with  a  -CH^N^ 
group  (AZMMTC)  was  scaled  up  to  the  20-gram  level  based  upon  a  unique  reaction 
mechanism  which  employs  acetyl  azide.  Attempts  to  convert  AZMTTC  to  the 
corresponding  amine  (N^  NH^)  via  a  variety  of  reactions  have  been  unsuccess¬ 
ful.  Efforts  continue  in  this  case  with  the  ultimate  objective  being  the  con¬ 
densation  of  the  amine  with  various  energetic  alcohols. 

The  bromo-precursor  (BrMTTC)  to  AZMTTC  has  not  been  found  to  be  especially 
reactive  to  most  standard  reagents.  Therefore,  the  iodo-precursor  (IMTCC)  is 
being  investigated  as  a  potential  intermediate. 

Attempts  to  synthesize  a  structurally  symmetrical  dibromo-HMX  derivative  from 
the  recently  reported  HMX  precursor  DADN  were  unsuccessful. 
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■PRIMARY  ACCOMPLISHMENTS 


CATALYTIC  EFFECT  OF  NF2*  VS.  NO 


1/T  X  104  CK'1) 


RDX  AUGMENTED  MORE  READILY  THAN  HMX  AT  LOT*7  PRESSURE  DUE  TO 
I<UCH  T7EAKER  LATTICE  FORCES  IF  RDX  VERSUS  HMX ,  THUS  ALLOWING 
FASTER  RATE  OF  RDX  GASIFICATION. 


pressure.  mpi 


FIGURE  2. 
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FUNDAMENTAL  STUDIES  OF  AZIDE  DECOMPOSITION  &  COMBUSTION 

Joseph  E.  Flanagan  and  Milton  B.  Frankel 

Rocketdyne  Division 
Rockwell  International 
Canoga  Park,  CA  91304 

The  utilization  of  azide  (-N^)  ingredients  is  rapidly  increasing  in  a  wide 
spectrum  of  Air  Force  applications.  Azide  materials  are  primarily  characterized 
by  their  high,  positive  heats  of  formation,  but  several  unique  combustion 
features  have  been  observed  such  as  higher  burning  rates  and  extinguishment 
at  pressures  substantially  above  ambient. 

This  study  will  elucidate  the  decomposition  and  combustion  characteristics  of 
azide  compounds  and  the  interaction  of  the  compounds  with  selected  solid  pro¬ 
pellant  ingredients.  The  final  output  from  the  study  will  be  twofold: 

1.  Definition  of  the  azido-chemistry ,  which  leads  to  the  desired  end 
product  (stop-start,  higher  burn  rate,  etc.) 

2.  Insight  into  suggested  new  chemical  compounds  (structures),  which  would 
yield  the  desired  results  and  serve  as  a  guide  for  future  synthesis  efforts 

The  experimental  investigation  will  include: 

1.  High  heating  rate  pyrolysis  investigations  at  pressures  ranging  from  atmospheric 
up  to  at  least  1.38  MPa.  Data  analysis  will  be  performed  via  mass  spectrometer ,  gas 
chromatographic  techniques,  and  isotopic  labeling,  to  elucidate  the  azide-chemistry 
and  key  ractive  species  that  control  the  combustion  process.  The  following  three 
categories  of  materials  will  be  studied: 

a.  Pure  azide  polymers  (mono  and  diazido  types) 

b.  Azide  polymers  plasticized  with  inert  or  energetic  materials 

c.  Solids  loaded  azide  polymer/plasticizer  systems  containing  oxidizers  (HMX, 

RDX,  AP,  and  AN)),  metallic  fuels,  and  alkali  seedants. 

2.  Combustion  investigations  will  be  carried  out  in  standard  strand  burners,  windowed 
combustion  vessels,  and  micromotors  on  each  of  the  three  categories  described 

above  to  delineate: 

a.  Regions  of  stable  combustion  in  both  strand  and  micromotor  conditions 

b.  Close  in  (less  than  2  mm)  flame  zone  characteristics  (both  temperature 
and  geometry) 

c.  Fuel  and  seedant  behavior  (consumption)  at  the  burning  surface 

d.  Degree  of  ionization  and  related  flow-field  conductivity 

3.  Mechanistic  descriptions  of  the  important  rate  controlling  and  intermediate 
product  producing  reaction  steps  will  be  developed. 

Sufficient  data  should  be  generated  to  aid  the  combustion  modelling  community 
in  model  formulation  and  validation. 


*# 
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SYNTHESIS  OF  DIFLUORAMINOXY- ,  DIFLUORAMINO-  or  FLUOROD I AZON I UM-CONTA I N E NG 
MATERIALS 

Jean'ne  M.  Shreeve 

University  of  Idaho 
Moscow,  Idaho  83843 

Future  rocket  propulsion  advances  require,  stable  high  energy  pro¬ 
pellant  oxidizers.  The  research  is  concerned  with  the  synthesis  of  stable, 
highly  oxidizing  compounds  which  contain  nitrogen- fluorine  bonds,  such  as 
d i f luor aminoxy ,  d i f luoramino ,  or  f luorod iazonium  groups  and  the  subsequent 
investigation  of  their  physical  and  chemical  properties. 

Several  routes  to  these  compound  types  are  currently  under  investiga¬ 
tion  or  will  be  attempted  during  the  next  several  months. 

(a)  The  Lewis  acid  catalyzed  introduction  of  -ONF2  into  F-olefins  and 
functionalized  F-olefins  via  reactions  of  NF3O. 

(b)  Photolysis  or  thermolysis  of  N2F4  with  hypofluor ites  and/or 
hypochlorites  to  form  mono  or  bis  -ONF2  compounds,  RfONF2  or  RfCF(ONF2)2- 

(c)  Mixed  -ONF2  and  -NF2  compounds  via  F-acyl  fluorides  and  HNF2 . 

(d)  Deoxygenation  of  N- f luoroazoxy  compounds  to  f luorodiazonium  com¬ 
pounds,  RfN=NF. 

(e)  Synthesis  of  RfN=NR  compounds  where  R  is  a  good  leaving  group 
which  can  be  displaced  by  fluorine. 

New  compounds  containing  the  N ,N-dif luorohydroxylamino  group  have  been 
prepared  by  the  addition  of  NF3O  to  olefins  in  the  presence  of  catalytic 
amounts  of  a  Lewis  acid.  The  new  compounds  RfONF2  (Rf  =  CF(CF3)C5Fi  jn, 
-CF(CFa)C(Q)F,  CF(CFt)SFs,  -CF( CF3 ) CF2OCF2CF( S02F)CF3 , 

-CF(CF3)CF26CFOCF(CF3)CF2CFlCF3) ,  -CF2CF2OC2F5 ,  CF2CF20CF2CF(CF3)OC3F7n) 
were  obtained  and  charac te> ized .  The  Markownikov  addition  products  may  be 
explained  by  a  simple  electrophilic  addition  of  NF2<1+  to  the  polar  double 
bond,  followed  by  fluorination  of  the  terminal  carbonium  ion.  The  anti- 
Markownikov  products  arise  from  NF2 0+  addition  to  the  tautomeric  form  of 
the  olefin  (F2C-CF=0Rf)  in  which  the  normal  polarity  of  the  double  bond  has 
been  reversed. 

Hexa f l uor ocyc lobutene  and  oc t a f I uorocyc 1 opent ene  were  reacted  with 
nitrosyl  chloride  or  nitrogen  dioxide  in  the  presence  of  KF  in  acetonitrile 
to  give  80%  yields  of  the  blue  hept af luoronit rosocyc lobutane( 1 )  and 
nonafluoronitrosocyclopentane  (II),  respectively.  Thermal  decomposition  of 

( I )  and  ( 1 1 )  in  Pyrex  glass  resulted  in  the  analogous  nitro  compounds, 
CK2CF2CF2SFNO2  and  CF2CF2CF2CF2CFNO2 .  Cycloaddition  react  ions  of  ( I )  and 

(II)  with  tetraf  luoroethene  gave  the  oxazet id ines ,  5F2CF2CIF2CFNOCF2CF2  and 
CF2CF2CF7CF2CFNOCF2CF2 ,  and  with  1 ,3-hex a f luorobutad iene ,  formed  the 
oxazines,  CF2CF2CF2tFNOCF2CF=CFCF2  and  fcF2'CF2CF2CF2£FNOCF2CF=CFCF2 •  With 
tetrafluorohydrazine  in  metal  reaction  vessels,  either  (I)  and  (II)  or  the 
cycloolefins  gave  di f luoroamines ,  CF2CF2CF2tFNF2  and  iF2CF2CF2CF2CFNF2 . 
However,  in  the  presence  of  Pyrex  glass,  N-heptafluorocyclobutyl-N'-fluoro- 
diimide  N-oxide  and  N-nonaf luorocyc lopentyl-N 1  - f 1 uorod i  imide  N-oxide 
(RfN(0)=NF)  resulted  from  heating  (I)  and  (II)  with  N2F4. 
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STABLE  HIGH  ENERGY  PROPELLANT  OXIDIZERS 
NEW  NITROGEN-FLUORINE-CONTAINING  COMPOUNDS 

PROPERTIES 

STABLE  TO  HEAT,  SHOCK,  MOISTURE 
OXIDIZING 

ROUTES 

LEWIS  ACID  CATALYZED  ADDITIONS 
OF  NF3O  TO  OLEFINS 

IN  SITU  GENERATION  OF  AN 

ACTIVE  FLUOR  I  NAT  I NG  REAGENT 

TETRAHYDRAZINE  REACTIONS  IN  THE 
PRESENCE  OF  SILICA 
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ACCOMPLISHMENTS 

ANTI MARKOWN I KO V  ADDITION 

+  F  F 

+  NFv0  AsFr  - >  FoNOC-CR 

2  6  2  F  F 

R  ~  -0C2F3/  Cl/  Br 
-0CF2CF(CF3)0C3F7 

MARKOWNIKOV  ADDITION 

F  F 

FC-CONFg 
F  F  2 

R=  CF3,  C5F1r  C(0)F ,  SF5 
CF20CF2CF(S02F)CF3 
CFoO(:F(CF,.)CFoOCF(CF7) 


**  1984  ROCKET  RESEARCH  MEETING  ** 
Abstract  60  Pe  1 


SYNTHESIS  AND  CHEMISTRY  OF  POLYNITROALKANES 
AND  POLYNITROOLEFINS 

Clifford  D.  Bedford  and  Robert  J.  Schmitt 


SRI  International 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 


Research  programs  in  energetic  materials  synthesis  and  mechanisms 
have  been  hampered  by  the  lack  of  general  methods  for  synthesizing 
aliphatic  polynitro  compounds  and  geminal  dinitroalkenes .  Efforts  to 
identify  and  overcome  these  research  barriers  have  thus  far  been  directed 
toward  special  situations  and  have  met  with  only  limited  success. 

The  objective  of  this  research  program  is  to  provide  a  better  under¬ 
standing  of  the  chemistry  of  polynitro  and  fluoronitro  compounds  and  to 
lead  to  the  design  of  new  classes  of  energetic  materials  having  the  desired 
properties  for  propellant  and  explosives  applications.  Basic  questions  con 
ceming  the  nature  of  specific  molecular  requirements  and  reaction  condi¬ 
tions  giving  to  mild  nitration  organic  materials  will  be  answered  through 
the  development  and  understanding  of  these  new  energetic  unsaturated  inter¬ 
mediates  and  nitration  techniques. 

During  the  first  six  months  of  this  contract  we  have  been  concerned 
with  developing  new  synthesis  methods  for  preparing  gem-dinitroolefins  and 
gem-dinitroalkanes .  In  line  with  this  objective,  a  new  nitrodealumination 
reaction  was  developed  for  the  preparation  of  nitroolef ins .  This  novel 
reaction  will  lead  into  the  more  extensive  area  of  nitrodemetallation 
reactions . 

Work  to  date  is  summarized  below: 

•  We  demonstrated  the  first  example  of  a  nitroolefin  preparation 
by  a  nitrodealumination  reaction. 

•  New  synthesis  of  nitroolefins  were  developed  using  the  olefin- 
aluminate  complex  in  the  nitrodealumination  reaction. 

•  Electrochemical  reduction  of  bromoalkanes  in  the  presence  of 
NO  leads  to  the  alkane  instead  of  the  nitroalkane. 
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The  vinylaluminum  intermediate  was  prepared  via  a  known  procedure. 
The  mixture  of  I  and  II  is  then  treated  with  tetranitromethane  (TNM) 
[equations  (2)  and  (3)]. 


Ph-CiC-H  +  <1-Bu)2A1H 


52  °C 


2-3  Hr*. 
Hexane/Toluene 


oc  + 

Ph-CSC-Al(i-Bu) 

\ 

Al(i-Bu) 

'I  * 

II 

672 

33% 

(1)  TNM 


-20°C  to  RT 

(2)  H20 


(2) 


(1)  N02BF4 

- H - -  No  Reaction 

(2)  H20 


NMR  and  mass  spectral  analyses  of  partially  purified  reaction  mixtures 
indicated  the  possible  presence  of  various  cross-coupled  (phenyl  and 
isobutyl  groups  in  the  same  molecule)  compounds  containing  nitro,  oxime, 
hydroxyl,  and/or  carbonyl  groups.  Efforts  are  underway  to  characterize 
these  interesting  and  unusual  product  mixtures.  No  6-nitrostyrene  was 
evident  in  the  product  mixture  resulting  from  the  reaction  with  N0oBF  . 
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Preparative  gas  chromatography  of  the  vinyl  aluminate  salt  reaction 
mixture  afforded  adequate  separation  of  several  reaction  products.  Com¬ 
bined  analyses  of  mass  spectral,  IR,  and  NMR  data  of  the  volatile  reaction 
products  have  enabled  us  to  characterize  the  following  compounds  in 
addition  to  8-nitrostyrene: 

FhC=CCH(OH)CH(CH3)2  PhCH2CH2N02  PhCH(i-Bu)CH2N02 

III  IV  V 

Hypothetical  mechanisms  for  the  formation  of  these  products  are  indi¬ 
cated  in  Scheme  1. 
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MODIFICATION  OF  PROPELLANT  BINDER  NETWORK 
FOR  IMPROVEMENT  OF  MECHANICAL  PROPERTIES 

AFRPL  Contract  F - 461 1  - 8 3- COOOl 

C.  Sue  Youn  Kim 

California  State  University,  Sacramento 


One  of  the  disadvantages  of  enerqetic  prepolymers  such 
as  g ly ci dy la z idopol ymer s  (GAP)  or  polydin itropr opyl  acrylates 
(PDNPA)  is  that  the  fraction  of  the  polymer  chains  bearing 
the  load  applied  to  the  binder  network  is  far  less  than 
that  of  conventional  prepolymers  such  as  polyethylene  qlycols 
(PFG).  For  example,  only  40  weight  %  of  GAP  will  participate 
in  the  strength  and  strain  capabilities  during  deformation 
in  contrast  to  91  weight  %  of  PEG.  Hence,  these  energetic 
prepolymers  usually  give  poorer  mechanical  properties  than 
PFG  when  formulated  into  elastomeric  networks. 

The  objective  of  our  project  was  to  evaluate  potential 
improvements  in  the  str ess / stra in  behavior  of  energetic 
propellant  binders  by  bimodal  network  formation;  i.e.  by 
blending  these  enerqetic  short  chains  with  long  chain  pre¬ 
polymers  prior  to  cure  into  rubbers.  The  extent  of  improve¬ 
ments  of  the  ultimate  properties  were  examined  in  terms  of 
the  chemical  type,  the  length  and  the  relative  amount  of 
the  long  chain  polymers  present,  and  also  of  the  visco¬ 
elastic  properties.  The  program  objective  and  the  results 
are  summarized  in  Figure  1. 

The  extent  of  improvement  appears  to  be  strongly  depen¬ 
dent  on  the  chemical  type  of  polymer  systems.  The  PEG/PEG 
system  shows  negligible  improvement  in  comparison  with 
PDNPA/PEG  systems.  For  example,  both  PEG  600/PEG  8000  and 
PEG  1000/PEG  8000  show  little  improvement  at  30  weight  per¬ 
cent  PFG  8000,  whereas  PDNPA(M-1 3 ) /PEG  8000  and  GAP /PEC  8000 
systems  with  short  chain  lengths  comparable  to  those  of 
short  chain  PEG,  show  two-  to  threefold  increase  in  strain 
(Figure  2).  Analysis  of  data  indicates  that  the  differences 
of  improvement  of  these  systems  appear  to  be  mainly  due  to 
the  degree  of  nonaffine  deformation  (lack  of  extension  of 
short  chains  is  compensated  by  the  neighboring  long  chains), 
rather  than  due  to  the  long  chains  actinq  as  inhibitor  of 
the  qrowth  of  rupture  nuclei.  A  more  flexible  crosslinker 
would  allow  more  extensive  nonaffine  deformation  at  a  given 
strain  rate  and  thus  further  improve  the  ultimate  properties. 

Other  effects  such  as  strain  induced  crystallization, 
filler  content,  arid  strain  rate  are  also  being  investigated. 
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Figure  1:  OBJECTIVE  AND  RESULTS  OF  THE  PROGRAM 


Objective ;  Imrovement  of  mechanical  properties  of 
propellant  binder  through  analysis  of 
bimodal  and  multimodal  elastomeric  network 
concept . 


Results : 

I.  GAP/PEG  8000  bimodal  binder  networks  give  signifi¬ 
cantly  improved  mechanical  proper  ites  over  monomod.il 
GAP. 

II.  PDNPA/PEG  8000  and  R-45M/HTPB  systems  also  show 
significant  improvements. 

III.  PEG/PEG  bimodal  systems  show  nj>  advantage. 

IV.  All  the  bimodal  systems  give  results  which  are 
consistent  with  "weakest  link"  theory. 

V.  Unfilled  bimodal  systems  show  appreciable  improve¬ 
ments  in  both  stress  and  strain  properties  over 
short  chain  monomodal  system.  As  filler  content 
increases,  the  favorable  effect  on  stress  tends  t.o 
diminish  while  that  on  strain  is  retained. 
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Figure  2  .  Improvement  of  Ultimate  Strain  Properties  of  Various 
Polymer/PEG  8000  Bimodal  Systems  vs  Wt  %  PEG  8000  at 
66°C;  Strain  Rate  4.2  x  10-4  m/s;  OPEG  600  (40  sk. 
atoms),  QpEG  1000  (70  sk.  atoms),  A  GAP  (70  sk.  atoms), 
X  PDNPA  (50  sk.  atoms) 
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INITIAL  THERMOCHEMICAL  DECOMPOSITION  MECHANISMS 
OF  ENERGETIC  INGREDIENTS:  DEUTERIUM  ISOTOPE  EFFECTS  AND 
ISOTHERMAL  REACTION  STUDIES 

Scott  A.  Shackelford,  Michael  B.  Coolidge, 

Stephen  L.  Rodgers,  Berge  B.  Goshgarian 

Air  Force  Rocket  Propulsion  Laboratory 
Edwards  AFB,  CA  93523 

The  thermochernical  decomposition  mechanism  and  kinetics  of  selected  energetic 
ingredients  are  being  investigated.  Special  attention  is  being  given  to  the  very  early  stages 
of  the  thermochemical  decomposition  reactions  since  earlier  investigations  with  TNT 
revealed  that  key  chemical  steps  can  occur  prior  to  detectable  exothermic  activity.  We  are 
using  deuterium  isotope  effects  in  conjunction  with  isothermal  thermochernical  analyses  to 
identify  in-si tu  and  in  real  time,  the  key  chemical  mechanisms  and  associated  bond  ruptures 
which  control  complex  thermochernical  decomposition  and  deflagration  processes.  Three 
tasks  compromise  this  effort,  and  two  are  in  progress. 

Task  I  The  thermochernical  decomposition  mechanism  and  kinetics  of  HMX  and  HMX-d8 
were  Investigated  collaboratively  with  R.N.  Rogers,  J.L.  Janney,  and  M.H.  Ebinqer  of  the 
Lps  Alamos  National  Laboratory.  Isothermal  DSC  (IDSC)  investigations  were  conducted  2- 
7°C  below  the  melting  point  of  HMX  (553-548°C),  and  three  distinct  physical  states  resulted 
during  the  overall  decomposition:  a  solid,  a  mixed  melt,  and  a  liquid  phase.  Three  different 
deuterium  isotope  effects  (DIE)  were  found  depending  upon  which  physical  state 
predominated  during  decomposition.  The  solid  phase  produced  a  primary  DIE  indicating  that 
C-H  bond  rupture  was  the  rate-limiting  step;  the  mixed  melt  phase  produced  an  apparent 
inverse  isotope  effect  indicating  that  C-H  bond  contractions  occur  as  the  intermolecular 
hydrogen  bond  forces  are  overcome  during  the  initial  melting  process;  and,  the  liquified 
stage  afforded  a  value  consistent  with  a  secondary  DIE  caused  by  rate-controlling  C-N  bond 
ring  rupture.  While  slow  thermochernical  decomposition  in  an  IDSC  cell  produced  three  DIE 
results,  isothermal  pyroprobe  deflagration  at  similar  temperatures  produced  only  a  primary 
DIE;  but,  very  highly  confined  time  to  explosion  (Henkin)  tests  provided  evidence  that  an 
inverse  DIE  operates  in  this  thermal  explosion.  Preliminary  evidence  suggests  that  the  rate¬ 
determining  reactions  in  thermochernical  HMX  responses  may  depend  upon  the  predominant 
physical  state  that  results  from  a  given  type  of  thermochernical  event  (i.e.,  slow 
decomposition,  pyrolysis,  thermal  explosion).  Similar  investigations  will  be  conducted  with 
RDX/RDX-dg  and  DNNC/DNNC~2,  d4,  d6  compounds. 

Task  2  Isothermal  DSC  (IDSC)  investigations  of  deuterated  liquid  TNT,  HNBB,  and  HNS  are 
being  investigated  to  determine  systematic  relationships  between  ingredient  chemical 
structure  and  thermochernical  decomposition  mechanisms.  Earlier  IDSC/deuterium  isotope 
effect  investigations  at  another  AF  laboratory  (FJSRL)  revealed  the  thermochernical 
decomposition  of  TNT  and  its  "dimeric"  analogue  (HNBB)  both  are  controlled  by  C-H  bond 
rupture.  HNS  is  an  unsaturated  analogue  of  HNBB  and  is  significantly  more  thermally 
stable;  but,  its  decomposition  mechanism  is  unknown.  A  mechanistic  elucidation  with  HNS 
should  reveal  how  its  vinyl  C-H  bond  structure  improves  thermochernical  stability.  By 
employing  a  newly  reported  synthesis  procedure,  we  successfully  synthesized  the  deuterated 
HNS  (HNS-d2)  needed  for  this  study;  HNS-d2  previously  was  an  unknown  compound. 

Task  3  Investigations  to  use  elucidated  thermochernical  decomposition  mechanisms  in  the 
modification  of  thermochernical  rate  processes  are  scheduled  for  a  later  start. 

These  three  tasks  address  the  problem  of  systematic  thermochernical  properties  control  vja 
either  initial  molecular  structural  design  and/or  chemical  additive  modification  of  the 
inherent  thermochernical  decomposition  mechanisms. 


INITIAL  THERM0C11EMICAL  DECOMPOSITION  MECHANISM/KINETICS  ELUCIDATION  - ►-  •  RATE-CONTROLLING  STEP 


TASK  1:  MECHANISM/KINETICS  ELUCIDATION 


HNS—  d ^  HAS  BEEN  SYNTHESIZED: 
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NEW  SYNTHETIC  TECHNIQUES 

FOR  ADVANCED  PROPELLANT  INGREDIENTS: 

SELECTIVE  TRANSFORMATIONS  AND  NEW  STRUCTURES 

Scott  A.  Shackelford,  Robert  D.  Chapman, 

Stephen  P.  Herrlinger,  John  L.  Andreshak 

Air  Force  Rocket  Propulsion  Laboratory 
Edwards  AFQ,  CA  93523 


New  chemical  transformations  and  methods  are  being  investigated  which  potentially 
lead  to  energetic  target  compounds  with  new  and  novel  molecular  structures,  or  which 
improve  synthetic  routes  to  known  compounds  that  are  currently  expensive  propellant 
luxuries.  Three  of  four  scheduled  research  tasks  are  in  progress. 

Task  I  Nucleophilic  substitutions  of  alkyl  bromides  with  silver  triflate  were  conducted  to 
form  precursor  triflate  ester  compounds  for  conversion  into  key  energetic  ingredients. 
Triflate  ester  intermediates  provide  a  facile  and  selective  pathway  to  energetic  compounds 
not  available  by  direct  synthetic  routes.  The  scope  and  limitations  of  this  triflate  ester 
synthesis  have  not  been  investigated  systematically,  especially  with  dihaloalkanes.  TheOf,uT- 
dibromoalkanes  react  much  more  rapidly  with  silver  triflate  than  do  corresponding 
monobromoalkanes,  and  unlike  the  monobromides,  dibrornoalkanes  resist  unwanted  iso¬ 
merization  to  2  triflate  products  in  CCI^  solvent.  The  behavior  of  this  is  explained  by  an 
SNj-type  reaction  where  intramolecular  or  anchimeric  assistance  occurs  in  bromoalkyl 
monotriflate  product  formation,  followed  by  intramolecular  anchimeric  stabilization  by  the 
triflate  group  during  alky Iditrif late  product  formation.  There  is  also  a  greater  tendency  to 
form  the  ditriflute  product  as  chain  length  increases  in  theOQD -dibrornoalkanes  up  through 
1 ,5-dibroinopentane  where  this  effect  plateaus.  Higher  reaction  temperatures  render  the 
intermediate  triflate  product  unstable;  decomposition  with  polymerization  then  occurs  in 
CCI^  solvent  and  aromatic  alkylation  results  in  benzene  solvent. 

Task  2  The  reaction  of  alkenes,  alkyl  carbonyls,  or  organometallic  compounds  with 
energetic  triflate  esters  could  permit  selective,  direct  syntheses  of  energetic  alkanes  or 
alkenes  in  which  an  overall  anti-Markovnikov  introduction  of  energetic  groups  is  desired. 
Preliminary  reactions  between  fluorodinitroethyl  triflate  and  carbonyl  compounds  imply 
low-yield  formation  of  energetic  -unsaturated  ethers.  Reaction  between  nitronium 
triflate  and  organolithium  compounds  showed  preliminary  promise  as  a  selective  non-acidie 
nitration  technique. 

T ask  3  Collaborative  AFRPL/Pt  Loina  College  (D.  F.  Shellhamer)  studies  are  being 
c onduc ted  with  olkene  methoxyf luorinat ion  reactions  using  methanol ic  XeF^  suspensions. 
Elucidation  of  this  mechanism  could  permit  an  exclusive  regioselectlve  control  which 
produces  the  anti-Markovnikov  ether  adduct.  Substitution  of  energetic  alcohols  for 
methanol  could  directly  produce  stable,  partially  fluorinated  energetic  ethers  rather  than 
undesireable  acetal-type  products.  Reactions  with  norbornene  produced  different  fluoro- 
methoxynorbornane  isomers  identified  as  those  when*  the  methoxy  moiety  exclusively  added 
first  when  catalyzed  by  I3T  New  isomers  occurred  with  HF  catalysis  in  which  the  F  atom 
added  first.  This  is  consistent  with  our  initially  proposed  reaction  mechanism  in  which 
effecting  methoxy  addition  first  produces  the  desired  anti-Markownikoff  product.  Future 
studies  will  be  conducted  with  the  energetic  fluorodinitroethanol  substrate. 

1  ask  IV  The  potential  syntheses  and  characterization  of  energetic  bora/ine  compounds  is 
scheduled  for  a  later  start. 


We  are  addressing  the  necessary  elements  for  an  effective  synthesis  effort:  (a)  novel 
transformations,  (b)  improved  synthetic  routes,  and  (c)  new  target  compound  classes. 


APPROACH 

Task  I:  Reactions  of  triflate  salts  (selective  intermediate  syntheses) 
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MECHANISTIC  COMBUSTION  INTERACTIONS  OF 
ISOTOPJCALLY  LABELED  PROPELLANT  INGREDIENTS 


Scott  A.  Shackelford 

Air  Force  Rocket  Propulsion  Laboratory 
Edwards  AFB,  CA  93523 

We  are  elucidating  in-si tu  and  in  real  time  the  key  fundamental  chemical  reaction 
mechanisms  and  ingredient  interactions  which  control  the  global  combustion  process. 
Isotopically  labeled  propellant  ingredients  are  being  used  to  detect  potential  deuterium 
isotope  effects  and  to  determine  the  randomness  of  isotopic  scrambling  experienced  during 
combustion.  These  approaches  should  provide  the  mechanistic  chemical  reaction  data 
needed  to  address  current  technological  barriers  in  nitramine  burn-rate  limitations,  and 
scientific  propellant  performance  tailorability.  Early  investigations  are  being  conducted  in 
a  pressurized  window  bomb  with  pure  pressed  pellet  ingredients.  Mixed  strand  burn 
formulations  and  later  micromotor  firings  will  follow  to  elucidate  fundamental  chemical 
reaction  mechanisms  and  ingredient  interactions  under  the  ballistic,  pressure,  and 
temperature  conditions  encountered  in  a  rocket  combustor.  Two  of  three  tasks  are  in 
progress. 

Task  I  integrates  the  deuterium  isotope  effect  concept  directly  to  the  combustion  process. 
In  a  collaborative  venture  with  Thiokol  (D.  A.  Flanigan  and  R.  E.  Askins)  and  LANL  (R.  N. 


Rogers),  HMX  and  HMX-dg  pressed  pellet  bum  rates  were  determined  in  a  pressurized 
window  bomb.  Burn  rate  comparison  of  identically  synthesized  HMX  and  HMX-dg  samples 
revealed  a  primary  deuterium  isotope  effect  at  500  and  1000  psai  pressures.  This  strongly 
suggests,  under  these  static  combustion  conditions,  that  C-H  bond  rupture  during  solid  state 
thermochemical  decomposition  ultimately  controls  the  burn  rate  of  HMX.  By  applying 
deuterium  isotope  effects  to  the  actual  combustion  processes,  we  also  provided  the 
combustion  research  community  a  new  in-si  tu  scientific  method  for  investigating  flame 
reaction  mechanisms.  Further  window  bomb  pressed  pellet  experiments  will  be  conducted 
with  deuterium  and  non-deuterium  labeled  RDX  cmd  DNNC  samples. 

Task  2  addresses  intramolecular  diffusion  in  .a  rocket  combustor.  Individual  propellant 
ingredients  appropriately  labeled  with  JC,  N,  and/or  °0  will  be  fired  in  a  rocket 
micromotor,  and  the  kinetically  frozen  products  will  be  collected  for  analysis.  The  degree 
of  randomness  found  in  scrambled  isotopic  product  labels  compared  to  that  expected  if 
complete  thermodynamic  equilibrium  were  achieved  prior  to  nozzle  kinetic  freeze  out  will 
reveal  the  amount  of  intramolecular  diffusion  or  thermodynamic  versus  kinetic  control 
exhibited  by  a  given  ingredient.  HMX,  RDX,  and  DNNC  will  be  studied. 

Task  3  will  study  ingredient  -  ingredient  chemical  interaction  end  diffusion  control  in  mixed 
ingredients  and  formulated  propellants.  Initially,  window  bomb  strand  burns  with  mixed 
ingredients  will  be  studied  in  which  only  one  ingredient  is  deuterium  labeled,  and  burn  rates 
will  be  analyzed  to  detect  any  change  in  key_  rate^controJUng  processes  from  pure 
ingredients.  Micromotor  firings  will  follow  with  1 JC,  ^N,  or  aO  labeled  ingredients,  and 
collected  combusiton  products  will  be  analyzed.  The  degree  of  random  isotopic  product 
scrambling  will  be  compared  to  that  expected  if  a  complete  ingredient-ingredient 
equilibration  is  established  prior  to  nozzle  kinetic  freeze-out.  Bi-ingredient  formulations 
will  be  built  one  ingredient  at  a  time  until  formulated  propellant  samples  are  fired.  Vv  ;  are 
initiating  stand  burn  experiments  this  year  with  HMX  and  HMX-dg  sample  mixtures. 

These  three  tasks  address  a  mechanistic  need  to  (a)  identify  rate-limiting  molecular  bond 
ruptures,  (b)  determine  the  degree  of  molecular  kinetic  vs.  thermodynamic  control,  and  (c) 
measure  the  amount  of  ingredient-ingredient  diffusion  encountered  in  propellant 
combustion. 


•  DEUTERIUM  ISOTOPE  EFFECTS  (DIE)  IN  CHEMICAL  COMBUSTION  MECHANISMS  -  RESULTS 
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NEW  ENERGETIC  POLYMER  SYNTHESIS 

AND  BINDER  CONCEPTS 
Stanley  D.  Morse 

University  of  Dayton  Research  Institute 
Dayton,  OH 

New  polymeric  syntheses  and  chain  modification  techniques  are  being  investigated  to 
obtain  improved  mechanical  properties  of  known  energetic  binders,  new  energetic 
prepolymer  and  binder  systems,  and  correlations  of  polymeric  chemical  structure  to 
physical/mechanical  properties.  Its  approach  centers  on  new  energetic  polymer  syntheses 
via  cationic  polymerization  or  by  synthetic  modifications  of  known  polymers.  Two  tasks  are 
active. 

Task  I  seeks  to  improve  the  physical  and  mechanical  properties  of  known  energetic 
polyether  polymers  derived  from  glycidyl-substituted  monomers.  The  modified  copolym¬ 
erization  technique,  successfully  demonstrated  with  the  PECH  polymer  in  this  project's 
predecessor  work  unit  and  completed  during  the  initial  stages  of  this  work  unit  (J.  Polym. 
Sci.  -  in  press)  will  be  used  with  GAP,  PGN,  and  poly-G  polymer  systems  to  maximize 
stress,  strain,  and  modulus  properties.  Epibromohydrin  polymerizations  were  conducted  to 
assess  the  PEBH  modified  polymer  as  a  viable  precursor  for  the  aforementioned  energetic 
binder  materials;  limited  success  resulted. 

Task  2  will  seek  to  obtain  a  hydroxy-terminated  poly-dinitropropyl  vinyl  ether  polymer 
suitable  for  propellant  urethane  curing  procedures.  The  dinitropropyl  vinyl  ether  monomer 
is  obtained  in  a  high-yield  one-step  synthesis  and  has  successfully  been  polymerized  into  an 
energetic  non-hydroxyterminated  binder  for  pressed  HMX  explosive  fills.  Its  inherent 
superior  thermochemical  stability,  coupled  with  its  excellent  hydrolytic  stability  and 
chemical  compatibility  at  elevated  temperatures  with  HMX,  makes  it  a  very  attractive 
energetic  binder  candidate  for  minimum  smoke  propellants. 

Task  3  addresses  the  polymerization  of  any  new  energetic  monomer  structures. 
Polymerization  reactions  will  be  conducted  to  characterize  new  energetic  polymer 
structure,  functionality,  and  physical  properties.  Urethane  curing  at  various  cure  ratios  will 
be  conducted  with  promising  new  energetic  prepolymers;  then,  mini-dogbones  will  be 
fabricated  and  evaluated  to  determine  resultant  mechanical  properties. 

Task  4  employs  polymerization  syntheses  and  techniques  which  provide  structurally 
controlled  polymer  backbones  for  correlation  investigations  of  structure  and 
physical/mechanical  properties.  Unlike  Task  I  which  seeks  to  optimize  the  physical  and 
mechanical  properties  of  GAP,  PGN,  and  poly-G  via  individual  synthetic  modification 
procedures,  Task  4  seeks  to  modify  only  the  energetic  pendant  substituent  on  an  identically 
synthesized  polymer  backbone  structure.  Any  differences  in  mechanical  or  physical 
properties  of  GAP,  PGN,  or  poly-G  samples  will  be  caused  solely  by  the  -N3,  -0N02,  or  - 
OCH^CfNCMof7  groups  themselves.  Because  of  the  ease  with  which  a  pendanf  triflate^ group 
coula  be  displaced  by  energetic  but  weak  nucleophiles  (e.g.  -N~,  -0N02,  -OFDNE),  a  poly- 
glycidyl  triflate  ether  polymer  was  synthesized.  The  isolate  glycidyl  triflate  polymer 
afforded  a  number  average  molecular  weight  of  6IOOg/mole.  Its  synthesis  so  far  is  limited 
by  low  monomer  yields.  Reaction  between  glycidyl  bromide  and  silvertriflate  produced  this 
monomer  which  was  isolated  and  characterized  by  NMR,  IR,  and  Mass  Spectrometry.  While 
yields  are  low  to  date,  reaction  between  glycidyl  bromide  and  silver  triflate  demonstrated 
tlie  first  triflate  synthesis  in  which  the  reactant's  epoxide  bond  was  left  intact.  Other 
synthetic  methods  produce  epoxide  ring  opening  and  destroy  the  glycidyl  structure  required 
for  subsequent  polymerization  reactions. 

The  overall  goal  of  this  project  is  two-fold:  new  improved  energetic  polymer /binder  systems 
and  systematic  polymer  structure/mechanical  property  correlation. 
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APPENDIX  1 


PROPULSION  RESEARCH  GOALS 

GENERAL  GOALS:  Research  In  support  of  rocket  propulsion  Is  directed  at 
providing  a  fundamental  basis  for  a  new  generation  of  propulsion  concepts  and 
Improving  the  scientific  understanding  of  phenomena  associated  with  the 
generation  of  propulsive  power.  Research  Is  needed  from  molecular  to 
macroscopic  scales  In  areas  which  Include  plasma  acceleration,  dynamics  of 
rocket  combustion,  the  behavior  and  synthesis  of  advanced  energetic  materials, 
characteristics  of  exhaust  plume  formation  and  radiation,  and  the  dynamics  of 
advanced  propulsion  concepts.  Further  knowledge  Is  needed  In  these  areas  to 
meet  mid-term  and  long-term  Air  Force  technology  goals  related  to  Improving 
performance,  reliability,  penetrability,  and  durability  for  Air  Force  mission 
applications  involving  space  systems,  ballistic  missiles,  and  air  launched 
missiles.  However,  specific  provisions  are  also  Included  for  research 
directed  toward  scientific  opportunities  which  cannot  be  related  to  particular 
mission  areas  at  this  time. 

Future  Air  Force  space  missions  for  communications,  surveillance,  and  weapons 
systems  will  require  substantially  higher  power  density  and  greatly  Improved 
propulsion  for  orbit  raising.  Research  Is  expected  to  lead  to  Improvements  of 
20%  In  propulsion  efficiency  with  Improved  liquid  and  solid  propellants  and 
100  to  400%  with  non-conventlonal  systems  such  as  electromagnetic  plasma 
acceleration  and  energy  beaming,  magnetoplasmadynamlcs  and  nuclear  power. 


SPECIFIC  GOALS: 

Ihe  following  topics  summarize  many  of  the  current  goals  which  are  being 
addressed  by  the  ongoing  research  programs. 

Solid  Propellant  Mechanics  To  Improve  understanding  of  damage  tolerant  design 
for  new  materials.  To  Investigate  fatigue,  damage  propagation  and  detection, 
fracture,  and  effects  of  spectrum  loading  and  sequencing.  To  extend  the 
existing  theories  to  three  dimensions,  Including  the  Interactions  between 
viscoelastic  behavior  and  damage,  the  mlcromechanlcs  of  propellant  structure, 
and  effects  of  temperature  to  obtain  a  approaches  that  can  predict  propellant 
stress-strain  behavior  from  relatively  simple  tests.  To  establish  an  age-life 
model  based  on  chemical  reactions.  To  establish  the  relationship  between 
carbon/carbon  material  processing  variables  and  the  kinetics  of  reaction  of 
these  materials  with  gases  representative  of  rocket  exhausts. 

Nonlnterfer Inq  Diagnostic  Techniques  To  achieve  sensors  which  will  enable  the 
adaptive  control  and  autonomous  operation  of  high  performance  propulsion 
systems.  To  achieve  diagnostic  techniques  for  measuring  the  gas  and 
gas-particle  flow  properties  representative  of  rocket  systems  Including 
plumes.  To  devise  diagnostic  techniques  to  obtain  rapid  and  quantitative 
spatially  and  temporally  resolved  measurements  of  temperature,  pressure, 
velocity,  species  concentrations  and  densities,  and  particle/droplet  size  and 
distribution  In  optically  dense,  high-temperature  reacting,  unsteady, 
turbulent  flowing  media.  To  establish  techniques  for  obtaining  rate  date  for 
energetic  material  being  heated  In  excess  of  105  K/s. 
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Exhaust  Plumes  To  provide  mechanistic  basis  for  analyzing.  Interpreting,  and 
altering  plume  radiance.  Establish  the  Influence  of  composition  and  flow  en¬ 
vironment  on  particle  formation  and  growth  mechanisms.  To  verify  and  Improve 
predictive  models  for  nucleatlon  and  condensation  In  multi -component,  multi¬ 
phase  systems.  To  determine  reaction  pathways  and  to  obtain  chemical  rate 
data  for  combustion  chamber,  nozzle,  and  exhaust  plume  flows.  To  Improve  and 
devise  diagnostic  methodology  for  determining  the  spatial  and  temporal  distri¬ 
bution  of  Important  species  so  that  rate  processes  can  be  measured.  To  deter¬ 
mine  the  mechanism  of  exhaust  plume  Igntlon  In  the  presence  of  flame  sup¬ 
pression  additives. 

Solid  Propellants  and  Energetic  Materials  To  control  the  behavior  of 
energetic  materials  during  combustion.  To  establish  relationships  between 
molecular  structure,  decomposition  and  combustion  processes.  To  synthesize 
energetic  Ingredients  leading  to  tough,  high  energy  propellants  with 
deslreable  properties  which  Include  low  glass-transition  temperatures,  high 
thermochemical  stability  temperatures,  high  solid  loadings,  and 
smokelessness.  To  Improve  techniques  to  determine  polymer  properties. 

Combustion  and  Reacting  flows  To  understand  Important  bond  rupture  and  heat 
release  mechanisms  so  that  research  to  modify  burning  rate  and  combustion 
efficiency  can  be  directed  at  specific  reaction  sites.  To  explore  methods  for 
obtaining  the  key  reaction  parameters  for  energetic  materials  being  heated  at 
rates  In  excess  of  10^K/s.  To  determine  the  effects  of  acceleration-fields 
on  propellant  ballistics.  To  characterize  the  mechanisms  leading  to  changes 
In  solid  propellant  properties  and  burning  modes  that  produce  transition  from 
normal  deflagration  to  detonation  of  rocket  motor  grains.  To  increase  the 
scope  and  usefulness  of  the  JANAF  Thermochemical  Property  Data. 

Combustion  Instability  To  Improve  analytical  methods  for  interpreting, 
predicting,  and  avoiding  Instability  behavior  of  solid  motors  through 
mechanistic  understanding  of  contributions  such  as  nozzle  damping,  acoustic 
eroslvlty,  pressure  coupling,  vortex  streaming,  vortex  sheading,  distributed 
combustion,  particle  and  structural  damping,  and  high  velocity  effects.  To 
acount  for  the  more  significant  effects  of  2-0  and  3-D  geometries.  To  defeat 
combustion  Instablltly  by  purposely  dissipating  acoustic  energy  in  prescribed 
modes.  To  examine  more  direct  means  of  measuring  acoustic  admittance.  To 
conduct  verification  of  new  combustion  response  measurement  techniques,  lo 
Investigate  methodologies  for  measuring  the  unsteady  velocity,  pressure,  and 
temperature  components  In  multi-dimensional  unsteady  reacting  flows.  To 
visualize  nonsteady,  multi-phase  flow  and  condensed-phase  breakup  for  the 
purpose  of  understanding  the  role  of  transient  processes  or.  acoustic  energy 
gains  or  losses.  To  determine  and  characterize  acoustic  energy  gains  and 
losses  In  terms  of  flow  field  parameters.  To  understand  contributions  of 
heterogeneous  structures  In  composite  propellants. 

Hetal  Combustion  To  characterize  the  basic  mechanisms  and  chemistry  Involved 
In  the  formation  of  metal  oxide  particulates.  To  experimentally  determine 
physical  and  chemical  processes  leading  to  the  acceleration  of  metal  droplet 
Ignition  and  combustion.  To  explore  methods  of  greatly  reducing  the  size  of 
metal  oxide  agglomerates  entering  nozzle  convergent  sections.  To  Improve 
methods  for  simultaneously  measuring  temperature  and  concentration  of 
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turbulent  flame  zones  of  multi -phase  media  so  that  conditions  leading  to  metal 
com-  bustlon  under  solid  propellant  ram-rocket  conditions  (particularly  high 
altl-  tude)  can  be  characterized.  To  assess  the  role  of  fluorine  In  the 
combustion  of  fluorlnated  metallic  propellants  and  to  determine  mechanisms 
that  affect  oxide  particle  size. 

Beamed  Energy  To  quantify  mechanisms  of  laser  energy  transfer  to  working 
fluids,  to  assess  the  barriers  to  optical  access  for  energy  transmission,  and 
to  establish  approaches  to  plasma  confinement.  To  enable  consideration  of 
candidate  short  wavelength  energy  beams  with  hydrogen;  to  determine  If  other 
working  fluids  and  fuels  offer  promise. 

Solar  Thermal  Propulsion  Establish  the  mechanisms  of  solar  radiation  energy 
absorption  under  the  desired  thruster  conditions. 

Electromagnetic  Acceleration  To  understand  magnetoplasmadynamlc  processes  for 
the  purposes  of  establishing  the  conditions  of  sustained  high  energy  density 
operation  and  to  project  upper  limits  of  power  density.  To  establish  the 
mechanism  leading  to  extended  electrode  life.  To  extend  operational  time  of 
electrode  and  contact  surfaces.  To  produce  stable  magnetoplasmadynamlc  flows 
for  sustained  high  energy  density  operation.  To  produce  conditions  that 
vaporize  and  disperse  metallic  fuels  being  considered  for  use  In  low  Earth 
orbits.  To  overcome  cooling  limitations  on  electrical  propulsion  systems 
operating  continuously  at  megawatts  power  levels.  To  realize  more  weight 
efficient  and  less  vulnerable  heat  rejection  concepts  for  electrodes, 
chambers,  and  nozzles. 

Ultra-High  Energy  density  To  decrease  the  losses  In  thermionic  flows.  To 
explore  novel  approaches  to  high  energy  density  storage.  To  devise  means  of 
accommodating  plasma  temperatures  greatly  In  excess  of  material  limits. 
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